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THERMO-MECHANICAL  AND  THERMAL 
BEHAVIOR  OF  HIGH-TEMPERATURE  STRUCTURAL  MATERIALS 


PREFACE 


Technical  ceramics  because  of  their  chemical  inertness,  high  melting 
point,  good  wear  resistance,  excellent  mechanical  stability  at  high  temp¬ 
erature  and  other  unique  -properties,-  represent  a  class  of  materials  eminently 
suited  for  many  critical  engineering  applications.  Unfortunately,  because 
of  their  brittleness  and  unfavorable  combination  of  pertinent  material 
properties,  technical  ceramics  generally  are  highly  susceptible  to  catas¬ 
trophic  failure  in  non-uniform  thermal  environments,  which  give  rise  to 
thermal  stresses  of  high  magnitude. _ 

Thermal  stress  failure  analysis  of  structural  materials  represents 
a  multi-disciplinary  problem  which  involves  the  principles  of  heat  trans¬ 
fer,  mechanics  and  materials  engineering.  Over  the  last  few  decades  much 
general  understanding  of  the  nature  of  thermal  stress  failure  of  brittle 
materials  has  been  generated.  However,  due  to  the  multi-disciplinary 
nature  of  the  problem,  the  ability  to  predict  thermal  stress  failure  quan¬ 
titatively  for  design  or  other  purposes  has  lagged  behind  the  progress  made 
in  other  engineering  fields.  >Jhe  objective  of  the  present  program  is  to 
improve  the  qualitative  and  quantitative  understanding  of  the  nature  of 
thermal  stress  failure  of  brittle  structural  materials,  including  the  ex¬ 
perimental  as  well  as  theoretical  variables.  In  order  to  achieve  this 
objective,  the  participating  investigators  and  scope  of  the  program  are 
organized  such  that  full  advantage  is  taken  of  the  combined  inputs  from  a 
number  of  engineering  disciplines.  In  a  similar  spirit,  a  number  of  studies 
were  conducted  in  cooperation  with  investigators  at  other  institutions. 

The  effort  of  this  program  consists  of  four  main  themes,  including: 
experimental  thermal  shock  testing  with  supporting  analyses,  measurement 
of  thermophysical  properties  relevant  to  thermal  stress  failure,  the 
analysis  of  mechanisms  of  thermal  stress  failure  and  the  dissemination 
of  information  on  thermal  stresses  in  the  form  of  review  articles,  con¬ 
ferences,  etc. 
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A  total  of  sixty-four  (64)  published,  accepted  or  submitted  technical 
publications  have  resulted  from  this  research  program  since  its  initiation. 
The  status  of  these  publications  is  summarized  in  Table  1. 

Studies  completed  within  the  period  covered  by  this  report  are  pre¬ 
sented  as  individual  chapters  in  the  main  body  of  this  report.  The  title 
of  these  chapters  with  a  brief  comment  are  as  follows: 


Chapter  I:  H.  Tawil,  L.  D.  Bentsen  and  D.  P.  H.  Hasselman,  "Effect 

of  Heat  Treatment  on  the  Thermal  Diffusivity  of  an  Alumino¬ 
silicate  Fiber-Reinforced  Carbon  Matrix" 

This  study  showed  that  composites  of  a  carbon  matrix  reinforced  by  alu¬ 
minosilicate  fibers  exhibited  an  irreversible  increase  in  thermal  diffusivity 
when  thermally  cycled  to  temperatures  in  excess  of  600°  C.  X-ray  analysis 
showed  that  this  observation  is  attributable  to  crystallization  of  both  the 
carbon  matrix  and  the  aluminosilicate  fibers. 


Chapter  II:  L.  D.  Bentsen,  0.  P.  H.  Hasselman,  "Role  of  Porosity 

in  the  Effect  of  Microcracking  on  the  Thermal  Conductivity 
of  Brittle  Ceramic  Composites" 

This  study  revealed  that  the  formation  of  microcracks  and  their  effect 
on  thermal  conductivity  is  a  very  sensitive  function  of  the  pore  content. 
This  effect  arises  because  for  pore  sizes  below  a  minimum  value,  the  stress 
intensity  factor  due  to  internal  stresses  resulting  from  thermal  expansion 
mismatches  is  too  low  to  result  in  microcrack  formation. 


Chapter  III:  L.  D.  Bentsen,  D.  P.  H.  Hasselman,  J.  J.  Brennan,  "Radia¬ 
tive  Contribution  to  the  Thermal  Diffusivity  and  Conducti¬ 
vity  of  a  Silicon  Carbide  Fiber  Reinforced  Glass-Ceramic" 

It  was  found  experimentally  that  the  thermal  diffusivity  and  conductivity 
of  a  lithium  aluminosilicate  glass-ceramic  with  or  without  SiC-fiber  rein¬ 
forcements  increased  with  increasing  specimen  thickness.  These  observations 
indicate  the  existence  of  a  radiative  contribution  to  the  total  heat  transfer. 


Chapter  IV:  R.  Ruh,  L.  D.  Bentsen,  D.  P.  H.  Hasselmen,  "Thermal  Diffu¬ 
sivity  of  SiC-BN  Composites" 

Composites  of  hot-pressed  Si C  with  a  dispersed  phase  of  BN  showed  signi¬ 
ficant  anisotropy  in  thermal  diffusivity  relative  to  the  hot-pressing  direc¬ 
tion.  This  effect  is  due  to  the  preferred  orientation  of  the  BN  particles 
promoted  during  hot-pressing.  The  anisotropy  in  thermal  diffusivity  was 
maintained  to  temperatures  as  high  as  1400°  C. 


Chapter  V:  L.  D.  Bentsen,  T.-Y.  Tien,  D.  P.  H.  Hasselmen,  "Effect  of 
Crystallization  of  the  Grain  Boundary  Phase  on  the  Thermal 
Diffusivity  of  a  Sialon  Ceramic" 

The  results  of  this  study  showed  that  the  degree  of  crystallinity  of 
the  grain  boundary  phase  has  a  significatn  effect  on  the  thermal  diffusivity 
of  a  sialon  ceramic. 


Chapter  VI:  J.  P.  Singh,  G.  Ziegler,  D.  P.  H.  Hasselman,  "Effect  of 
Drop-Height  on  the  Critical  Temperature  Difference  (aT  ) 
for  Brittle  Ceramics  Subjected  to  Thermal  Shock  by  Quench¬ 
ing  into  Water" 

The  critical  temperature  difference  (AT)  required  to  initiate  thermal 
stress  fracture  in  a  number  of  brittle  ceramics  subjected  to  thermal  shock 
by  quenching  into  water  was  found  to  be  independent  of  drop-height.  These 
findings  were  interpreted  in  terms  of  the  relative  contributions  to  the  heat 
transfer  during  the  quench. 


Chapter  VII:  T.  D.  Nguyen,  J.  R.  Thomas,  Jr.,  D.  P.  H.  Hasselman, 
"Effect  of  Thermal  Diffusivity  on  Magnitude  of  Thermal 
Stress  in  a  Brittle  Ceramic  Subjected  to  Rapid  Thermal 
Cycling" 

An  analysis  was  conducted  of  the  magnitude  of  thermal  stresses  in  a 
ceramic  component  subjected  to  thermal  cycling  of  such  high  frequency  that 
thermal  equilibrium  between  thermal  cycles  is  not  achieved.  It  was  found 
that  under  such  conditions  the  magnitude  of  thermal  stress,  in  addition  to 
the  usual  properties  such  as  coefficient  of  thermal  expansion.  Young's  modu¬ 
lus  and  Poisson's  ratio,  also  depends  simultaneously  on  the  magnitude  of 
the  thermal  diffusivity  and  thermal  conductivity.  It  was  also  found  that 
the  magnitude  of  the  thermal  stresses  decreased  as  the  frequency  increased. 


Chapter  VIII:  H.  Hencke,  J.  R.  Thomas,  Jr.,  D.  P.  H.  Hasselmen,  "Role 
of  Material  Properties  in  the  Thermal  Stress  Fracture  of 
Brittle  Ceramics  Subjected  to  Conductive  Heat  Transfer" 

The  magnitude  of  thermal  stresses  in  a  ceramic  component  subjected  to 
heating  or  cooling  by  direct  thermal  contact,  in  addition  to  being  a  function 
of  the  coefficient  of  thermal  expansion.  Young's  modulus  and  Poisson's  ratio, 
also  was  found  to  be  a  function  of  the  thermal  conductivity,  thermal  diffu¬ 
sivity  and  specific  heat  of  both  the  component  as  well  as  the  surrounding 
medium.  Sample  calculations  indicated  that  thermal  stress  fracture  of  cer¬ 
amic  specimens  subjected  to  thermal  shock  by  quenching  in  a  fluid  is  the 
result  more  of  conductive  rather  than  convective  heat  transfer. 


Chapter  IX:  D.  P.  H.  Hasselman,  J.  P.  Singh,  "Criteria  for  the  Thermal 
Stress  Failure  of  Brittle  Structural  Materials" 

This  review  article  presents  all  known  information  on  the  material 
variables  which  affect  the  thermal  stress  fracture  of  brittle  ceramics. 
Appropriate  experimental  data  were  presented. 
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ABSTRACT 

The  effect  of  heat-treatment  on  the  thermal  diffusivity  c c 
an  aluminosilicate  fiber-reinforced  carbon  matrix  was  measured  by 
the  laser-flash  method  while  the  specimen  was  being  heat  treated. 
A  permanent  increase  in  thermal  diffusivity  was  observed,  which 
was  attributed  to  the  combined  effects  of  structural  modifica¬ 
tions  in  the  carbon  at  the  atomic  or  molecular  level,  a  change  in 
the  morphology  of  the  pore  phase  within  the  carbon  matrix,  and 
possibly  the  crystallization  of  the  fibers. 


*0n  leave  from  Societe  Europeenne  de  Propulsion,  Bordeaux, 
France . 
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Ceramic  fiber-reinforced  carbon  composites  are  of  interest 
for  many  industrial  and  aerospace  applications.  The  ceramic 
fibers  serve  to  enhance  fracture  toughness  and  thermal  insulating 
ability.  Heat- treatment  of  carbon  can  lead  to  irreversible 
changes  (increases)  in  the  thermal  conductivity  and  diffusivity 
due  to  atomic,  molecular  or  other  structural  rearrangement  (1,2). 
Such  effects  should  also  be  present  in  ceramic  fiber-reinforced 
carbon,  as  reported  in  the  present  paper. 
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2.  EXPERIMENTAL 

The  fiber-reinforced  carbon  material  used  for  this  study  was 
obtained  from  a  commercial  source*  in  the  form  of  a  single  block 
measuring  approximately  20  by  20  by  2  cm.  The  fibers  consisted 
of  an  aluminosilicate  wool+.  In  the  manufacture  of  the  composite 
the  aluminosilicate  fibers  were  first  densified  with  a  phenolic 
resin  followed  by  carbonization  which  transforms  the  resin  to 
coke.  The  composite  was  then  impregnated  with  a  petroleum  pitch 
and  once  more  carbonized  at  600  to  900°C  for  10  to  20  hours  to 
yield  the  carbon  matrix.  In  the  final  composite  the  volume  con¬ 
tent  of  the  aluminosilicate  fibers  was  approximately  45%.  The 
absence  of  x-ray  diffraction  peaks  for  the  composite  indicated 
that  both  the  carbon  matrix  as  well  as  the  aluminosilicate 
fibers  were  primarily  amorphous. 


♦CERASEP  A102-SEP  Industrie,  Bordeaux,  France. 

+Kaowool,  Morganite  Ceramic  Fibers,  S.A.  Liege,  Belgium. 
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The  bulk  density,  apparent  density,  and  open  porosity  of 
disc-shaped  samples  approximately  12  mm  in  diameter  by  2  mm  thick 
were  measured  by  Archimedes'  method  using  boiling  water  as  the 
fluid  medium.  The  data  obtained  for  the  densities  and  the  poros¬ 
ity  indicated  the  existence  of  considerable  heterogeneity  within 
the  block.  For  eight  immediately  adjacent  samples  cut  from  a 
single  cylinder  core-drilled  from  the  block,  the  bulk  density 
ranged  from  1.25  to  1.55  g/cmJ  with  a  mean  value  of  1.41  g/cm  . 
Similarly,  the  apparent  density  ranged  from  1.81  to  1.88  g/cm 
with  a  mean  value  of  1.85  g/cm 3  .  The  open  porosity  ranged  from 
16  to  33%  with  a  mean  value  of  24%.  However,  for  all  cores 
drilled  from  widely  different  areas  of  the  original  block,  the 
densities  and  porosity  showed  much  greater  variation,  ranging 
from  1.05  to  1.55  g/cm3  for  the  bulk  density,  and  from  1.46  to 
2.16  g/cm 3  for  the  apparent  density.  The  open  porosity  ranged 
from  13.5  to  51%. 

Fig.  1  shows  a  typical  scanning  electron  fractograph  of  a 
composite  sample.  For  comparison  and  interpretation  of  the  data, 
a  few  samples  of  the  carbon  matrix  without  the  ceramic  fibers 
were  included  in  the  study  as  well. 

The  effect  of  heat- treatment  on  the  heat  conduction  behavior 
was  determined  by  measurements  of  the  thermal  diffusivity  by  the 
laser-flash  method  (3)  on  the  same  samples  used  for  the  density 
measurements.  In  principle,  this  method  consists  of  subjecting 
one  side  of  a  thin  disc-shaped  specimen  to  a  single  uniformly 
distributed  laser-flash.  The  time  dependence  of  the  temperature 
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response  of  the  opposite  face  is  then  monitored.  The  thermal 
diffusivity  (k)  is  calculated  from: 


<  =  AL  /t 


where  L  is  the  specimen  thickness  and  ^1/2  is  the  Period  of 
time  following  the  laser-flash  for  the  rear  surface  to  reach 
one-half  of  the  final  temperature  rise  when  the  heat  from  the 
flash  is  uniformly  distributed  within  the  specimen.  A  is  a  con¬ 
stant,  the  value  of  which  is  determined  by  the  heat  loss  from  the 
specimen  during  the  laser-flash  as  described  by  Heckman  (4). 

The  laser-flash  technique  was  particularly  convenient  for 
the  present  study  as  the  thermal  diffusivity  could  be  measured 
while  the  specimens  were  being  heat  treated.  The  heat  treatment 
and  measurement  of  the  thermal  diffusivity  were  carried  out  in  a 
nitrogen  atmosphere  within  a  resistively  heated  carbon  furnace*, 
modified  for  the  laser-flash  diffusivity  method.  Because  the 
changes  in  specimen  thickness  could  not  be  monitored  during  heat 
treatment,  the  thermal  diffusivity  was  calculated  on  the  basis  of 
the  initial  specimen  thickness.  Two  heat  -  treatment  schedules 
were  followed.  The  first  consisted  of  heating  to  and  cooling 
from  a  pre-selected  temperature  over  a  total  period  of  some  5  to 
6  hours.  The  second  heating  schedule  consisted  of  heating  at  a 
rate  of  80°C/min  to  the  desired  temperature  followed  by  holding 
the  specimens  under  isothermal  conditions  for  periods  ranging 
from  2  to  3  hours. 


♦Model  1000A,  Astro  Industries,  Santa  Barbara,  CA. 
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RESULTS 


3  . 

The  thermal  diffusivlty  at  room  temperature  prior  to  heat 
treatment  of  adjacent  samples  cut  from  a  single  core  showed  no 
significant  variation,  despite  the  aforementioned  wide  variations 
in  bulk  and  apparent  density  and  porosity.  However,  at  a  given 
value  of  density  and  porosity,  specimens  cut  from  different  areas 
in  the  original  block  showed  considerable  variation  in  thermal 

diffusivity.  For  instance,  at  a  value  of  bulk  density  of  1.4 

3  -  3  ? 

g/cm  ,  the  thermal  diffusivity  ranged  from  2.0  x  10  cm  /s 

3  9 

3.3  x  10  cm  Vs.  These  variations  must  be  attributable  to  local 
variations  in  composition  and  microstructure.  For  this  reason, 
no  clear-cut  dependence  of  the  thermal  diffusivity  on  density 
could  be  established. 

Fig.  2  shows  the  thermal  diffusivity  for  two  samples  of  the 
fiber-carbon  composite  cycled  to  1000  and  1400°C  and  for  a  sample 
of  the  carbon  matrix  without  fibers  cycled  to  1400°C.  All  three 
samples  show  a  major  increase  in  the  thermal  diffusivity  on 
return  to  room  temperature.  Such  an  increase  at  least  in  part 
must  be  attributed  to  the  permanent  structural  modification  of 
the  carbon  matrix,  also  responsible  for  the  observations  of  Tay¬ 
lor  (1)  and  Wagner  (2)  for  the  effect  of  heat-  treatment  on  the 
thermal  conductivity  of  carbons.  Of  interest  to  note  was  that 
this  structural  modification  of  the  present  samples  did  not 
result  in  a  major  change  in  the  x-ray  diffraction  patterns  of  the 
carbon  matrix  following  heat- treatment .  Comparison  of  the  magni¬ 
tudes  of  the  thermal  diffusivity  of  the  composite  and  carbon  sam- 


pies  given  in  Fig.  2  indicates  that  the  effect  of  the  fibers  on 
the  thermal  insulating  ability  is  retained  following  thermal 
cycling . 


Fig.  3  shows  the  time  dependence  of  the  thermal  diffusivity 
of  fiber-carbon  samples  held  isothermally  at  600,  1000  and  1400°C 
and  for  a  carbon  sample  held  at  1400°C  following  heating  to  these 
temperatures  at  a  rate  of  80°C/min.  These  data  show  that  at  any 
temperature  the  thermal  diffusivity  is  invariant  with  time.  This 
suggests  that  the  structural  modifications  in  the  carbon  phase 
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responsible  for  the  increase  in  the  thermal  diffusivity  appear  to 
be  a  function  only  of  the  maximum  temperature  reached.  Such 
structural  changes  appear  not  to  be  the  result  of  a  thermally 
activated  process,  for  which  a  time  dependence  of  the  thermal 
diffusivity  at  any  given  temperature  would  be  expected.  In  this 
respect,  the  findings  of  this  study  confirm  those  of  Rappeneau 
and  co-workers  (5)  that  the  thermal  as  well  as  the  mechanical 
properties  are  not  influenced  by  the  rate  of  heating  during  car¬ 
bonization. 

The  observed  changes  in  the  thermal  diffusivity  with  heat - 
treatment  also  were  accompanied  by  weight  loss,  and  by  micros- 
tructural  and  dimensional  changes.  For  a  number  of  specimens 
held  at  1000°C  for  approximately  1  1/2  hours,  the  weight  loss 
ranged  from  3  to  12%.  The  corresponding  weight  loss  at  1400°C 
ranged  from  9  to  25%.  For  temperatures  below  1400°C,  the  weight 
loss  of  the  composite  samples  can  be  attributed  primarily  to  the 
weight  loss  of  the  carbon,  also  observed  by  Wallouch  and  Fair  (6) 
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and  Markovic  and  Marinkovic  (7).  For  the  composite  samples  at 
1400°C,  decomposition  of  the  aluminosilicate  fibers  probably  con¬ 
tributed  to  the  weight  loss  as  well.  Evidence  for  this  latter 
effect  is  given  in  Fig.  4  which  compares  the  morphology  of  fiber 
surfaces  before  and  after  a  2-hour  isothermal  heat- treatment  at 
1400°C.  For  the  fibers  of  Fig.  4,  Fig.  5  compares  the  relative 
intensities  of  the  A1  and  Si  in  the  fibers  as  determined  by  EDAX 
analysis  before  and  after  heat- treatment  at  1400°C  for  2  hours. 
The  decrease  in  the  amount  of  Si  relative  to  the  A1  suggests  that 
the  weight  loss  of  the  composite  also  can  be  attributed  to  the 
evaporation  of  the  silicon  m  the  form  of  one  (or  more)  of  its 
oxides . 

Heat-treatment  of  the  composite  samples  also  caused  a  sig¬ 
nificant  increase  in  the  open  porosity  For  four  samples  held  at 
1000°C  for  2  hours,  the  open  porosity  increased  from  a  mean  valu<- 
of  approximately  24%  to  about  33%  with  corresponding  changes  in 
the  bulk  and  apparent  densities  of  +7.4%  and  +19.7%,  respec¬ 
tively.  For  four  samples  held  at  1400°C,  the  mean  value  of  open 
porosity  increased  from  25%  to  35%,  with  corresponding  increases 
in  the  bulk  and  apparent  densities  of  1.9%  and  23.9%,  respec¬ 
tively.  No  significance  should  be  attached  to  the  changes  in 
these  properties  at  1000  and  1400°C  in  view  of  the  limited  number 
of  samples. 

A  change  in  open  porosity  implies  a  rearrangement  of  the 
pore  morphology.  Such  a  rearrangement,  however,  was  not  readily 
ascertained  by  scanning  electron  fractography  possibly  because  of 
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the  highly  heterogeneous  nature  of  the  microstructure.  As  indi¬ 
cated  by  the  theory  of  the  thermal  conductivity  of  composites, 
the  nature  of  the  distribution  of  second  phase  inclusions, 
including  pores,  has  a  significant  effect  on  the  composite  ther¬ 
mal  conductivity  and  diffusivity  (8).  For  this  reason,  the 
observed  changes  in  the  thermal  diffusivity  on  heat-treatment  at 
least  in  part  also  may  have  been  due  to  changes  in  the  pore 
structure.  A  quantitative  assessment  of  this  effect,  however,  is 
not  feasible  at  this  time  because  of  the  difficulty  in  obtaining 
a  quantitative  description  of  the  microstructural  features  of-  the 
composite  samples,  including  the  pores. 

Dimensional  changes  of  the  specimens  during  heat  treatment 
were  also  observed.  At  1000  and  1400°C  for  2  hours,  the  dimen¬ 
sions  decreased  by  an  average  of  about  4  and  5%,  respectively. 
Because  the  thermal  diffusivity  was  calculated  on  the  basis  of 
the  original  thickness  of  the  specimens,  the  data  shown  in  Fig. 
2,  with  the  exception  of  the  initial  data  point  at  room  tempera¬ 
ture,  represent  overestimates  by  no  more  than  10%. 

The  evaporation  of  the  Si  (or  oxide)  from  the  fibers 
resulted  in  a  change  in  their  composition.  An  x-ray  analysis  of 
a  sample  heat-treated  at  1400°C  for  2  hours  revealed  diffraction 
peaks  which  corresponded  to  the  alumina- si lica  compound  known  as 
mullite.  Crystalline  materials  generally  have  higher  thermal 
conductivity  than  the  corresponding  amorphous  phase.  The 
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improved  thermal  conductivity  of  the  fibers  is  thought  to  account 
for  the  observations  as  shown  in  Fig.  2  that  the  relative  perma- 


nent  increase  in  the  thermal  diffusivity  following  heat  treatment 
at  1400°C  is  higher  for  the  fiber-carbon  composites  than  for  the 
carbon  matrix  itself.  The  absence  of  a  pronounced  time  depen¬ 
dence  of  the  thermal  diffusivity  as  shown  in  Fig.  3,  suggests 
that  the  crystallization  of  the  fibers  is  essentially  complete 
when  the  1400°C  temperature  is  attained. 
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ABSTRACT 

Fracture  mechanical  principles  were  used  to  predict  the  role  of 
porosity  in  the  formation  of  microcracks  in  brittle  ceramic  compo¬ 
sites  containing  a  dispersed  phase  with  a  coefficient  of  thermal 
expansion  less  than  that  of  the  matrix.  It  is  shown  that  micro¬ 
cracking  will  occur  only  if  the  amount  of  porosity  lies  within  a 
region  Pi<P<Pa,  and  that  all  the  raicrocracks  will  be  of  equal  size. 
Using  this  criterion,  the  effect  of  the  microcracks  and  pores  on 
the  thermal  conductivity  is  examined.  The  expected  behavior  in¬ 
cludes  a  sharp  drop  in  the  relative  thermal  conductivity  and 
diffusivity  at  P^,  followed  by  a  porosity-independent  value  between 
and  Pa.  These  conclusions  were  confirmed  by  measurements  of 
the  thermal  diffusivity  of  composites  containing  silicon  carbide 
dispersed  in  a  matrix  of  magnesium  oxide,  beryllium  oxide,  and 
aluminum  oxide,  each  with  a  range  of  compositions  and  amounts  of 
porosity. 


INTRODUCTION 

Many  brittle  materials  such  as  geological  formations  and 
structural  ceramics  for  high  temperature  purposes  can  undergo 
extensive  microcracking.  The  principle  causes  for  the  formation 
of  microcracks  include  phase  transformations,  thermal  expansion 
anisotropy  of  the  grains  of  polycrystalline  aggregates,  or  mis¬ 
matches  in  the  thermal  expansion  of  the  individual  components 
within  a  composite.  These  latter  two  effects  result  in  the 
generation  of  internal  stresses  upon  changes  in  temperature  of  the 


polycrystal  or  composite  from  the  manufacturing  temperature.  The 
magnitude  of  these  stresses  is  a  function  of  the  range  of  tempera¬ 
ture  change,  the  degree  of  thermal  expansion  anisotropy  or  mismatch, 
and  the  elastic  properties  of  the  grains  or  individual  components. 
These  internal  stresses  can  be  large  enough  to  result  in  localized 
fracture  and  raicrocrack  formation.  Such  microcracks  can  have  a 
major  effect  on  the  continuum  behavior  of  brittle  materials,  af¬ 
fecting  the  elastic  moduli,  strength,  fracture  energy,  and  toughness. 
Microcracking  can  also  lead  to  significant  improvements  in  machin- 
ability  and  thermal  shock  resistance. 


Microcracks  represent  barriers  to  heat  flow  by  phonon  and 
electron  transport  and  have  a  significant  effect  on  the  thermal 
conductivity  and  diffusivity  at  temperatures  low  enough  that 
radiative  heat  transfer  across  the  cracks  does  not  make  a  major 
contribution.  The  continuum  heat  transport  properties  of  micro- 
cracked  materials  have  received  much  theoretical  and  experimental 
attention.  Microcracks  in  structural  ceramics  of  technical  inte¬ 
rest  were  found  to  decrease  the  thermal  conductivity  and  diffusi¬ 
vity  by  as  much  as  a  factor  of  three  (1-4).  A  number  of 
theoretical  expressions  for  the  effect  of  cracks  on  thermal 
conductivity  have  appeared  in  the  literature.  For  a  matrix  with 
randomly  oriented  penny-shaped  cracks  of  radius  b,  the  effective 
thermal  conductivity,  Keff  is  (5) 


Keff  "  K0<l+8Nb3/9)  '1 


(1) 


where  Kq  is  thermal  conductivity  of  the  crack-free  material  and 
N  is  the  number  of  cracks  per  unit  volume. 


Many  ceramic  materials  also  contain  a  pore  phase.  Porosity 
of  10  to  30%  is  desirable  for  improving  the  thermal  insulating 
ability  and  thermal  shock  resistance  of  refractory  ceramics  used 
in  furnace  linings  and  related  applications.  A  residual  pore 
phase  of  only  a  few  percent  in  structural  ceramics  for  high 
performance  applications,  such  as  components  of  high  temperature 
turbine  engines,  generally  results  from  incomplete  densif ication 
during  sintering  or  hot-pressing  of  compacts  made  by  powder 
metallurgical  techniques. 


For  temperatures  at  which  radiation  across  the  pores  is 
negligible,  a  porous  material  represents  a  special  case  of  a 
composite  consisting  of  a  continuous  matrix  with  a  dispersed  phase 
of  zero  thermal  conductivity.  Theoretical  solutions  for  the 
effective  thermal  conductivity  of  composites  indicate  that  the 
pore  phase  should  lead  to  a  decrease  in  thermal  conductivity, 
with  the  relative  decrease  determined  by  the  volume  fraction, 
shape,  and  orientation  of  the  pores  (6-7) .  Such  a  decrease  in 


thermal  conductivity  or  diffusivity  was  confirmed  by  experimental 
data  (7-8) .  For  the  morphology  and  orientation  of  the  residual 
pore  phase  found  in  sintered  powder  compacts,  the  relative 
decrease  in  thermal  conductivity  generally  was  found  to  be  of  the 
order  of  three  times  the  pore  volume  fraction.  Because  of  the 
accompanying  decrease  in  the  specific  heat  per  unit  volume,  the 
corresponding  relative  decrease  in  thermal  diffusivity  is  expected 
to  be  about  two  times  the  fractional  porosity.  For  spherical  pores 
the  relative  decrease  in  thermal  diffusivity  is  only  one-half 
the  pore  volume  fraction. 

Microcracked  brittle  materials  usually  also  contain  a  residual 
pore  phase.  For  this  reason,  the  thermal  transport  behavior  of 
ceramics  prone  to  microcracking  should  be  affected  by  the  combined 
presence  of  the  microcracks  and  the  pore  phase.  In  general,  at 
concentrations  sufficiently  low  that  the  local  temperature  fields 
around  second  phase  inclusions  or  other  thermal  discontinuities 
are  not  affected  by  neighboring  inclusions  or  discontinuities,  the 
relative  effects  of  microcracks  and  porosity  on  thermal  conductivity 
are  additive.  At  first  sight,  then,  it  appears  that  the  combined 
effects  of  the  cracks  and  the  pore  phase  on  the  heat  conduction  be¬ 
havior  of  microcracked  materials  could  be  ascertained  from  information 
on  the  separate  effects  for  the  cracks  and  pores,  as  determined 
either  from  theory  or  experiment. 

For  microcracked  materials,  however,  this  approach  is  compli¬ 
cated  by  the  very  role  of  the  pore  phase  in  the  formation  of 
microcracks.  For  any  residual  stress  distribution,  the  magnitude 
of  peak  stress  alone  is  not  sufficient  to  describe  the  generation 
of  microcracks.  The  appropriate  criterion  for  microcrack  formation 
is  that  the  magnitude  of  the  stress  intensity  factor,  which  is 
governed  by  the  magnitude  of  stress  and  size  of  the  precursor  crack, 
equals  or  exceeds  the  critical  stress  intensity  factor  (Kjc) 
appropriate  to  the  preferred  plane  of  crack  propagation,  which  may 
be  grain  boundaries  in  polycrystalline  materials  or  phase  boundaries 
in  compositess.  For  materials  prone  to  microcracking,  pores  repre¬ 
sent  the  very  source  of  microcrack  precursors,  especially  the  crack¬ 
like  pores  located  at  the  triple  points,  grain  boundaries,  and 
interfaces. 

The  fracture  mechanics  of  microcrack  stability  and  propagation 
has  received  considerable  theoretical  attention  (9-11) .  Figure  1 
illustrates  schematically  the  stability  of  a  circumferential  micro¬ 
crack  (i.e.,  pore)  around  a  spherical  inclusion  contained  within  a 
brittle  matrix.  The  thermal  expansion  coefficient  of  the  matrix 
is  greater  than  that  of  the  inclusion,  resulting  in  tensile 
tangential  stress  in  the  matrix  near  the  inclusion  when  the  composite 
is  cooled  from  the  fabrication  temperature,  at  which  the  internal 


crack  length  for  a  radial  crack  produced  by 


stresses  are  zero.  Figure  1  indicates  that  the  stress  intensity 
factor  (Kj)  is  a  function  of  the  ratio  of  the  crack  size  to  the 
inclusion  size.  Crack  instability  and  propagation  will  occur  only 
between  crack  sizes  c^  and  ca.  For  c<c^  no  microcracking  will 
occur  since  the  crack  size  is  so  small  that  Kj<KIc  for  the  maxi¬ 
mum  stress  involved.  The  absence  of  crack  propagation  for 
c<ca,  the  crack  length  for  crack  arrest,  is  due  to  the  combined 
effects  of  the  decrease  in  stress  with  increasing  distance  from  the 
inclusion,  and  the  stress  relief  due  to  the  formation  of  the  crack, 
which  permits  accomodation  of  the  mismatch  in  thermal  expansions 
of  the  matrix  and  inclusion.  The  conclusion  obtained  from  Figure 
1  is  critical  to  the  objectives  of  this  study  and  the  experimental 
results  to  be  presented.  The  figure  shows  that  the  crack  length 
for  crack  arrest,  ca,  is  independent  of  the  initial  crack  length 
c>ci.  Furthermore,  the  value  of  the  initial  crack  length  is  in¬ 
cluded  in  the  value  of  crack  length  at  crack  arrest.  It  is  for 
this  reason  that  the  combined  effect  of  porosity  and  microcracks 
are  expected  to  be  coupled  rather  than  simply  additive. 

A  direct  relationship  generally  exists  between  the  pore 
content  and  the  number  of  pores  per  unit  volume,  their  geometry 
and  their  dimensions.  For  simplicity,  it  will  be  assumed  that  the 
dimensions  and  geometry  of  all  pores  are  identical.  The  existence 
of  a  minimum  crack  size,  c^  (i.e.,  pore  size),  implies  that  no 
microcracking  will  occur  below  a  minimum  value  of  porosity,  Pj_. 

Thus  for  values  of  porosity  0<P<Pa,  the  thermal  conductivity  will 
be  governed  by  the  effect  of  the  pore  phase  only.  The  same  effect 
will  prevail  for  P>Pa  for  which  no  microcracks  can  form  either. 


For  Pa<P<Pa,  however,  microcrack  formation  will  occur,  with 


all  cracks  being  of  equal  size.  Because  the  size  of  the  micro¬ 
crack  precursor  (i.e.,  the  pore)  is  included  within  the  dimensions 
of  the  microcrack  at  arrest,  it  is  expected  that  for  P^<P<Pa,  the 
thermal  conductivity  will  be  governed  only  by  the  microcracks  and 
should  be  independent  of  pore  content. 


The  implications  of  the  above  conclusions  with  regard  to  the 
effect  of  porosity  on  the  thermal  conductivity  of  brittle  materials 
prone  to  microcracking  is  shown  schematically  in  figure  2.  For 
P<Pi  and  P>Pa,  the  thermal  conductivity  decreases  monotonically 


with  increasing  porosity.  At  P^,  the  thermal  conductivity  decreases 


discontinuous ly  to  a  value  which  remains  constant  for  P^<P<Pa. 


Clearly,  the  dependence  of  thermal  conductivity  on  pore  content 
shown  in  figure  2  is  valid  for  the  idealized  situation  of  all  pores 
having  equal  size  and  geometry.  In  practice  however,  the  pores  will 
exhibit  a  range  of  dimensions  and  geometries .  Depending  on  the 
size  and  shape  distributions,  the  dependence  of  thermal  conductivity 
on  pore  content  shown  in  figure  2  will  show  a  more  continuous 


POROSITY 

Predicted  effect  of  porosity  on  the  rela¬ 
tive  thermal  conductivity  of  brittle 
composites  prone  to  microcracking  due  to 
thermal  expansion  mismatch. 


variation  preceeded  by  an  initial  rapid  decrease  in  thermal  conduc 
tivity  at  P-P^.  Furthermore,  for  Pi<P<Pa,  a  slight  decrease  in 
the  thermal  conductivity  is  expected  for  those  small  pores  which 
did  not  contribute  to  the  formation  of  microcracks.  It  should  be 
noted  that  microcrack  formation  itself  also  contributes  to  the 
creation  of  porosity.  For  this  reason,  the  measured  pore  content 
represents  the  sum  of  the  original  pore  phase  plus  the  extra 
porosity  created  by  microcracks.  For  two-component  composite  ma¬ 
terials,  an  upper  bound  on  the  pore  content  created  by  the  micro¬ 
cracks,  (Pmc) ,  can  be  estimated  assuming  that  the  total  thermal 
expansion  mismatch  between  the  two  components  can  be  accounted  for 
by  the  volume  of  the  cracks.  This  yields: 


P  -  3AciAt 
me 

where  Act  is  the  mismatch  in  the  coefficients  of  thermal  expansion 
of  the  two  components  and  AT  represents  the  temperature  range  of 
cooling  from  the  processing  temperature. 


The  above  conclusions  regarding  the  effect  of  porosity  on  the 
thermal  conductivity  of  microcracked  materials  are  supported  by 
experimental  data  reported  herein. 


EXPERIMENTAL 

The  microcracked  brittle  composites  consisted  of  aluminum 
oxide,  beryllium  oxide,  or  magnesium  oxide  with  a  dispersed  phase 
of  silicon  carbide.  Table  I  lists  the  coefficients  of  thermal 
expansion  (12)  for  the  three  matrix  materials  and  the  silicon 
carbide,  together  with  the  values  of  thermal  conductivity  (13), 
density,  specific  heat  (14) ,  and  thermal  diffusivity  at  room 
temperature.  The  latter  four  values  are  required  for  the  calcula¬ 
tion  of  the  composite  thermal  conductivity  without  microcracks. 
Because  the  coefficients  of  thermal  expansion  for  the  three  c xide 
materials  exceed  the  value  for  the  silicon  carbide,  microcrack 
formation  will  occur  in  the  matrix  phase  when  the  composites  are 
cooled  from  the  manufacturing  temperature. 

The  specific  process  by  which  these  composites  were  made 
consisted  of  pressure  sintering  mixtures  of  the  appropriate  powders 
in  graphite  dies.  The  amount  of  residual  pore  phase  could  not  be 
closely  controlled  during  the  pressure  sintering  operation;  there¬ 
fore,  a  number  of  samples  were  obtained  to  assure  a  sufficient 
range  of  pore  content.  Typical  microstructures  for  composites 
consisting  of  magnesium  oxide  with  a  range  of  silicon  carbide 
contents  are  shown  in  figure  3.  The  mean  grain  size  of  the  magnes¬ 
ium  oxide  matrix  was  of  the  order  of  1  ym,  with  similar  values  for 
the  composites  based  on  aluminum  oxide  and  beryllium  oxide.  The 
silicon  carbide  particles  were  somewhat  disc-shaped  with  a  diameter 
of  about  6  urn  and  a  thickness  of  about  1-2  ym. 

The  porosity  of  the  samples  was  determined  from  the  measured 
density  of  the  specimens  and  the  theoretical  density  of  the  pore- 
free  composites  as  calculated  from  the  values  of  density  listed  in 
table  I.  The  effect  of  the  pores  and  microcracks  on  the  heat 
conduction  behavior  was  determined  by  measurements  of  the  thermal 
diffusivity  by  the  flash  method  (15)  using  a  glass-Nd  laser  as  the 
flash  source.  All  measurements  were  made  at  room  temperature  using 
a  liquid  nitrogen-cooled  infrared  detector  to  monitor  the  transient 
temperature  of  the  specimen.  The  specimens  were  either  circular 
disks  about  12  mm  in  diameter  or  square  plates  about  1  cm  on  a 
side,  each  about  2  mm  thick. 

The  specimens  for  all  three  series  of  composites  varied  in 
pore  content  as  well  as  silicon  carbide  content.  To  facilitate 
reporting  the  data,  the  actual  experimentally  determined  diffusivity 
data  were  multiplied  by  the  density  and  specific  heat  to  obtain 
the  thermal  conductivity.  These  values  of  thermal  conductivity 


were  then  dividied  by  the  value  of  the  conductivity  without  cracks. 
This  latter  value  was  calculated  using  the  Bruggeman  mixture 
equation  and  the  appropriate  values  of  the  thermal  conductivity 
of  the  individual  components  as  listed  in  table  I . 


RESULTS  AND  DISCUSSION 

Figure  4  compares  typical  experimental  data  for  the  thermal 
diffusivity  of  the  alumina-silicon  carbide  composites  as  a  function 
of  silicon  carbide  content  to  the  values  for  the  crackfree  composites 
calculated  from  composite  theory.  The  experimental  data  clearly  fall 
far  below  the  calculated  values.  The  difference  is  attributable  to 
the  formation  of  the  microcracks  due  to  the  thermal  expansion  mis¬ 
match.  Dividing  the  experimental  value  by  the  calculated  value  for 
any  given  silicon  carbide  content  yields  the  relative  thermal 
conductivity,  which  reflects  the  effect  of  the  microcracks  alone, 
without  the  accompanying  change  in  the  thermal  conductivity  due 
to  the  silicon  carbide  phase.  These  relative  values  can  then  be 
plotted  as  a  function  of  pore  content  in  order  to  verify  the  role 
of  porosity  in  the  heat  conduction  behavior  of  microcracked  com¬ 
posites. 

Figures  5,  6,  and  7  show  the  data  for  the  relative  thermal 
conductivity  of  the  composites  of  silicon  carbide  dispersed  in 
alumina,  beryllia,  and  magnesia,  respectively,  as  a  function  of 
the  pore  content.  All  three  sets  of  data,  especially  the  extensive 
set  of  data  for  the  MgO-SiC  system,  are  in  general  agreement  with 
the  behavior  shown  in  figure  2,  predicted  from  fracture-mechanical 
principles.  The  very  rapid  decrease  in  thermal  conductivity  with 
increasing  porosity  at  the  lower  values  of  porosity  is  clearly  evi¬ 
dent.  The  invariance  of  the  thermal  conductivity  with  porosity  at 
the  higher  pore  contents  is  compatible  with  the  conclusion  that  the 
final  crack  size  is  independent  of  the  initial  pore  size  (i.e.,  ini¬ 
tial  pore  content  at  a  given  crack  density) .  The  data  scatter  for 
the  thermal  conductivity  at  the  lowest  values  of  porosity  is  too 
large  to  make  a  reliable  extrapolation;  nevertheless,  the  general 
trend  of  the  data  suggests  the  existence  of  a  minimum  value  of  poro¬ 
sity  below  which  no  microcracking  occurs. 

The  amount  of  porosity  created  by  the  formation  of  the 
microcracks  cannot  be  ascertained  from  the  present  data.  However, 
the  differences  in  porosity  at  a  given  value  of  silicon  carbide 
content  must  represent  the  amount  of  original  pore  phase  prior  to 
the  formation  of  the  microcracks. 

Comparison  of  the  data  shown  in  figures  5,  6,  and  7  indicates 
that  the  relative  thermal  diffusivity  at  the  higher  pore  contents  is 
lowest  for  the  magnesium  oxide  matrix,  followed  by  the  beryllium 
oxide  and  then  the  aluminum  oxide  matrix.  This  is  expected  from 
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the  differences  in  the  mismatch  in  the  coefficients  of  thermal  ex¬ 
pansion  between  the  silicon  carbide  and  the  three  matrices.  The 
extent  of  crack  propagation  and  resulting  decrease  in  thermal 
diffusivity  is  expected  to  increase  with  increasing  thermal  expan¬ 
sion  mismatch. 

In  summary,  on  the  basis  of  the  stability  of  microcracks  in 
internal  stress  fields  produced  by  thermal  expansion  heterogeneity, 
a  prediction  was  made  concerning  the  relative  effect  of  the  pore 
content  on  the  thermal  diffusivity  and  conductivity  of  brittle 
materials  subject  to  microcracking  and  was  confirmed  by  experiment. 
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ABSTRACT 

The  thermal  diffusivity  and  conductivity  of  a  silicon  carbide 
fiber-reinforced  lithium  aluminosilicate  glass-ceramic  was  measured 
using  the  laser-flash  method.  As  indicated  by  the  effect  of  speci¬ 
men  thickness  on  thermal  conductivity  and  its  positive  temperature 
dependence,  heat  transfer  by  radiation  makes  a  significant  contri¬ 
bution  to  the  total  thermal  conductivity.  The  observed  dependence 
of  the  effect  of  specimen  thickness  on  thermal  conductivity  sug¬ 
gests  that  radiative  heat  transfer  between  the  carbon  coated 
surfaces  of  dielectric  materials  for  the  laser-flash  technique  may 
be  partially  governed  by  the  view  factor  between  the  specimen  sur¬ 
faces. 


INTRODUCTION 

Electrons,  phonons,  and  photons  represent  the  primary  mechan¬ 
isms  of  heat  transport  in  solids.  Each  of  these  mechanisms  dis¬ 
plays  a  unique  temperature  dependence  (1).  In  dielectric  materials 
the  thermal  conductivity  due  to  electron  transport  is  generally 
not  significant,  and  above  room  temperature  the  conductivity  is 
approximately  inversely  proportional  to  the  temperature  as  the 
result  cf  increased  phonon  collisions. 

The  thermal  conductivity  due  to  photon  transport  is  a  function 
of  the  optical  absorption  and  emission  behavior.  For  highly  trans¬ 
parent  materials,  heat  transfer  can  occur  by  direct  radiation  be- 
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tween  the  surfaces  of  the  material.  In  this  case  the  effective 
radiative  thermal  conductivity  for  surfaces  with  equal  emissivity 
is  (2) : 

K  =  (r2-)  4an2T3d  (1) 

R  Z”6 

where  e  is  the  emissivity  of  the  specimen  surface,  o  is  the 
Stefan-Boltzman  constant,  n  is  the  index  of  refraction,  T  is  the 
absolute  temperature,  and  d  is  the  thickness  of  the  material  in 
the  direction  of  heat  flow. 


Heat  transport  by  radiation  can  also  occur  by  emission, 
absorption,  and  re-emission.  When  the  distance  between  absorption 
and  re-emission  of  a  photon  is  small  compared  to  the  specimen 
dimensions,  the  radiative  conductivity  becomes  (2) : 

KR=^On2T31  (2) 

where  1  is  the  photon  mean-free-path  which  is  equal  to  the  recipro¬ 
cal  of  the  extinction  coefficient. 

Equations  1  and  2  indicate  that  the  thermal  conductivity  due 
to  radiation  between  the  specimen's  surfaces,  and  the  thermal 
conductivity  due  to  absorption  and  re-emission  of  radiation  are 
proportional  to  T3.  This  implies  that  the  radiative  thermal  con¬ 
ductivity  should  exhibit  a  strong  positive  temperature  dependence. 
Equations  1  and  2  also  show  that  the  thermal  conductivity  due  to 
radiation  directly  across  the  material  is  directly  proportional  to 
the  specimen  thickness,  whereas  radiative  heat  transfer  by  photon 
absorption  and  re-emission  is  independent  of  the  specimen  size. 

In  general,  the  conduction  of  heat  in  dielectric  materials 
occurs  by  the  combined  contributions  of  phonon  and  photon  transport. 
Above  room  temperature,  the  negative  temperature  dependence  of  the 
thermal  conductivity  due  to  phonon  transport  will  be  offset  by  the 
positive  dependence  of  the  radiative  contribution.  In  crystalline 
dielectric  materials  which  have  a  large  phonon  mean-free-path, 
such  as  aluminum  oxide,  beryllium  oxide,  and  diamond,  the  negative 
temperature  dependence  of  phonon  transport  generally  overrides 
the  positive  temperature  dependence  of  radiative  heat  transfer  up 
to  temperatures  of  some  1000°C.  In  highly  amorphous  dielectrics, 
phonon  transport  tends  to  be  suppressed.  For  these  materials, 
the  relative  contribution  of  radiative  heat  transfer  can  be  such 
that  the  combined  thermal  conductivity  due  to  phonons  and  photons 
produces  a  significant  positive  temperature  dependence. 

Such  a  positive  temperature  dependence  was  observed  in  a 
feasibility  study  of  the  composite  method  for  determining  the 
thermal  conductivity  of  amorphous  silicon  carbide  fibers  contained 
within  a  continuous  matrix  of  a  lithium  aluminosilicate  (LAS) 


glass -ceramic  (3) .  The  experimental  data  indicated  that  radiative 
heat  transfer  within  the  fibers  as  well  as  the  glass-ceramic 
matrix  made  significant  contributions  to  the  total  thermal  conduc¬ 
tivity  of  the  composite.  A  few  exploratory  measurements  indicated 
the  existence  of  the  effect  of  specimen  thickness  on  the  thermal 
conductivity,  suggestive  of  radiative  heat  transfer  between  the 
specimen  surfaces  as  described  by  equation  1.  The  purpose  of  this 
study  was  to  investigate  the  effect  of  specimen  size  on  the  thermal 
conductivity  of  the  SiC-LAS  composites  in  order  to  determine  in 
further  detail  the  nature  of  the  radiative  contribution  to  the 
total  thermal  conductivity. 


EXPERIMENTAL  PROCEDURES 
Materials 

The  silicon  carbide  fibers  from  a  commercial  source*  were 
made  from  an  organometallic  polymer  by  the  method  of  Yajima  (4) . 
Chemical  analysis  indicated  that  the  fibers  consisted  of  about 
65%  SiC,  25%  Si02 ,  and  10%  C.  The  fibers  were  in  the  form  of  tows 
of  yarn  containing  about  500  fibers  per  tow.  The  average  fiber 
diameter  was  approximately  10  ym  and  the  density  was  about  2.55 
g*cm~3.  The  degree  of  peak  broadening  in  x-ray  analysis  showed 
that  the  fibers  were  nearly  amorphous  with  a  crystallite  size  of 
2.5  to  3.0  nm. 

The  matrix  material  consisted  of  a  lithium  aluminosilicate 
glass-ceramic  with  a  composition  nearly  identical  to  that  of  a 
commercial  glass-ceramic**;  however  the  3%  Ti02  nucleating  agent 
was  replaced  by  about  2%  Zr02  in  order  to  eliminate  the  reactivity 
between  the  TiC>2  and  the  SiC.  The  density  of  the  LAS  was  2.52 
g*cm~3. 

Composites  were  fabricated  by  passing  the  SiC  yarn  through  a 
slurry  of  the  uncrystallized  glass  powder  dissolved  in  isopropyl 
alcohol.  After  drying,  the  coated  yarn  was  cut  to  fit  a  graphite 
die,  and  was  hot-pressed  in  vacuum  for  5  to  60  minutes  at  1400  to 
1500°C  and  14  MPa.  Following  hot-pressing,  the  glass  was  crystal¬ 
lized  by  heat  treating  the  composite  for  1  to  2  hours  at  tempera¬ 
tures  of  880  to  1100°C.  Both  uniaxial  and  0/90  cross-ply  composites 
were  fabricated  with  SiC  contents  up  to  50  vol.%.  Figure  1  shows 
a  photomicrograph  of  a  composite  with  uniaxial  fibers,  and  figure 
2  shows  a  scanning  electron  fractograph  of  a  0/90  cross-ply  com¬ 
posite. 

*N ippo n  Carbon  Company,  Japan. 

**C-9608,  Corning  Glass  Works,  Corning,  New  York. 


Fig.  1.  Photomicrograph  of  a  composite  of  LAS  glass 
ceramic  containing  48  vol.%  uniaxially 
aligned  SiC  fibers. 


Measurement  and  Evaluation  of  Thermal  Properties 


The  laser-flash  method  (5)  was  used  to  measure  the  thermal 
diffusivity  of  the  composites  and  the  matrix.  The  specimens  of 
the  composites  cut  from  the  hot-pressed  blanks  were  in  the  form 
of  nearly  square  plates  about  9.5  mm  on  a  side  with  thicknesses 
ranging  from  0.5  to  4.5  mm.  Samples  were  cut  to  permit  measurement 
of  the  thermal  diffusivity  both  parallel  and  perpendicular  to  the 
fibers.  Specimens  of  the  LAS  matrix  were  circular  disks  12.8  mm 
in  diameter  with  thicknesses  of  0.3  to  4.0  mm.  For  all  samples, 
direct  transmission  of  the  laser  beam  was  prevented  by  coating  the 
specimen  surfaces  with  carbon.  The  transient  temperature  response 
of  the  specimen  was  monitored  with  remote  optical  sensors.  For 
measurements  above  room  temperature,  the  specimens  were  supported 
within  a  carbon  resistance  furnace  containing  a  nitrogen  atmosphere 

In  the  analysis  of  the  data,  corrections  for  heat  loss  were 
taken  into  account  using  the  method  of  Heckman  (6) .  Changes  in 
the  thickness  and  density  of  the  specimens  due  to  thermal  expansion 
were  also  accounted  for.  The  specific  heat  of  the  fibers  and  the 
glass-ceramic  matrix  was  determined  with  differential  scanning 
calorimetry. 
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Fig.  2.  Scanning  electron  fractograph  of  a  composite 
of  LAS  glass-ceramic  containing  45  vol.%  of 
crossplied  SiC  fibers. 


From  the  measured  value  of  the  thermal  diffusivity  of  either 
the  fiber-free  matrix  or  the  composite,  the  corresponding  value  of 
the  thermal  conductivity  can  be  calculated  from: 

K  =  Kpc  (3) 

where  K  is  the  thermal  diffusivity,  p  is  the  density,  and  c  is  the 
specific  heat.  The  specific  heat  of  the  composite  was  calculated 
using  the  rule  of  mixtures  and  the  measured  specific  heats  of  the 
fibers  and  the  matrix  shown  in  Figure  3 . 

The  thermal  conductivity  and  diffusivity  of  the  fibers  were 
evaluated  from  the  corresponding  data  for  the  matrix  and  the  com¬ 
posites  using  the  theory  for  the  thermal  conductivity  of  mixtures 
with  a  dispersed  phase  in  the  shape  of  cylinders  (7).  For  heat 
flow  parallel  to  uniaxially  aligned  fibers,  the  thermal  conductivity 
of  the  composite  is: 


K  =  K  V  +  KV 
c  mm  f  £ 


where  K  is  the  thermal  conductivity,  V  is  the  volume  fraction,  and 
the  subscripts  c,  m,  and  f  refer  to  the  composite,  matrix  and 
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Fig.  3.  Specific  heat  of  LAS  glass-ceramic  and  of 
amorphous  SiC  fibers. 


fibers,  respectively.  For  heat  flow  perpendicular  to  the  fibers, 
the  thermal  conductivity  is: 


After  determining  the  thermal  conductivity  of  the  fibers  using 
equations  4  and  5,  the  thermal  diffusivity  can  be  calculated  by 
means  of  equation  3. 


RESULTS  AND  DISCUSSION 

Figure  4  shows  the  data  for  the  effect  of  specimen  thickness 
on  the  thermal  conductivity  and  diffusivity  of  the  LAS  glass- 
ceramic  matrix  and  two  composites  with  heat  flow  parallel  and  per 
pendicular  to  the  fiber  axes.  The  thermal  conductivity  was  cal¬ 
culated  using  Eq.  5  and  the  experimental  data  for  the  thermal 
diffusivity,  density,  and  specific  heat.  The  data  indicate  a 
significant  specimen  size  effect  for  these  three  materials,  sug¬ 
gestive  of  radiative  heat  transfer  between  the  carbon-coated 
surfaces  of  the  specimens. 

A  similar  size  effect  is  observed  for  the  values  of  thermal 


Fig.  4.  Thickness  dependence  of  the  thermal  diffusivity 
and  conductivity  at  25°C  for  LAS  glass-ceramic 
and  two  composites  of  LAS  with  SiC  fibers  oriented 
parallel  (middle)  and  perpendicular  (bottom)  to  the 
direction  of  heat  flow. 


THERMAL  DIFFUSIVITY 
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conductivity  and  diffusivity  of  the  SiC  fibers  at  room  temperature 
as  shown  in  Figure  5.  These  values  were  calculated  from  the  ex¬ 
perimental  data  for  the  composite  and  matrix  samples  using  equations 
4  and  5.  The  values  are  essentially  independent  of  the  direction 
of  heat  flow  relative  to  the  fibers.  This  suggests  that  any 
existing  interfacial  effects  do  not  influence  heat  flow  perpendi¬ 
cular  to  the  fibers. 

For  a  quantitative  interpretation  of  the  size  effect,  it  is 
important  to  note  that  the  thermal  conductivity  does  not  rise 
linearly  with  specimen  thickness,  which  is  expected  from  Eq.  1  if 
all  the  heat  transfer  through  the  specimen  occurred  by  radiation 
between  the  surfaces.  This  indicates  that  additional  mechanisms 
of  heat  transfer  independent  of  the  specimen  thickness  also  con¬ 
tribute  to  the  thermal  conductivity. 

Such  mechanisms  include  heat  transfer  by  phonons  and  by 
continuous  emission  and  absorption  of  photons  with  short  mean- 
free-path.  Both  processes  are  independent  of  specimen  thickness. 
However,  an  additional  factor  to  explain  the  observed  size  effect 
should  be  considered  as  well.  First  it  should  be  noted  that 
equation  1  for  the  effective  thermal  conductivity  due  to  radiative 
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SAMPLE  THICKNESS  (mm) 

Thermal  diffusivity  and  conductivity  of  SiC 
fibers  at  room  temperature  determined  from 
experimental  values  for  the  LAS  matrix  and 
the  composites. 


heat  transfer  between  sample  surfaces  is  valid  when  the  ratio  of 
the  length  of  the  surface  to  the  distance  between  the  surfaces 
approaches  infinity.  Under  this  condition,  all  the  radiation 
emitted  by  one  surface  is  absorbed  by  the  second  surface.  This 
condition  is  not  met  for  the  specimens  used  in  this  study,  where 
in  the  case  of  the  thicker  specimens,  the  ratio  is  only  two  to 
three.  Thus  an  appreciable  fraction  of  radiation  emitted  hemis- 
pherically  from  the  hotter  surface  will  not  be  intercepted  and 
absorbed  by  the  colder  surface.  This  condition  is  governed  by  the 
view  factor  from  one  surface  to  another,  which  is  defined  as  the 
fraction  of  the  total  radiation  emitted  by  the  one  surface  that  is 
directly  incident  on  the  other  (8) .  For  parallel  circular  disks 
of  equal  diameter  positioned  directly  opposite  each  other,  the 
view  factor  (F)  can  be  written  (9) : 


where  r  is  the  radius  of  the  disks  and  d  is  the  distance  between 
them. 

The  effect  of  specimen  geometry  on  the  radiative  thermal  con¬ 
ductivity  of  typical  carbon-coated  laser-flash  specimens  can  be 
taken  into  account  by  incorpora ting  the  view  factor  of  equation  6 
into  equation  1,  which  yields: 

K  =  (~~)  4F0n2T3d  (7) 

R  2-e 

Figure  6  shows  the  normalized  radiative  conductivity,  (2-e)/ 
4eOn2T3 ,  as  a  function  of  thickness  for  a  range  of  specimen  dia¬ 
meters,  which  for  the  LAS  matrix  in  this  study  is  12.8  mm.  Com¬ 
parison  of  the  data  for  the  thermal  conductivity  as  a  function  of 
specimen  thickness  shown  in  Figure  4  with  the  curve  for  specimen 
diameter  of  12.8  mm  shown  in  Figure  6  indicates  that  the  variation 
in  the  observed  thermal  conductivity  can  be  partially  accounted 
for  by  the  change  in  specimen  shape  with  specimen  thickness.  Such 
an  effect  should  be  taken  into  account  in  the  interpretation  of 
the  data  for  the  thermal  diffusivity  of  transparent  materials 
measured  by  the  laser-flash  method. 

Figure  7  shows  the  temperature  dependence  of  the  thermal 
conductivity  for  two  values  of  thickness  of  specimens  of  the  glass- 
ceramic  matrix  and  two  composites  with  heat  flow  parallel  and  per¬ 
pendicular  to  the  fiber  axes.  The  increase  in  thermal  conductivity 
with  temperature  for  all  three  materials  is  indicative  of  the 
contribution  of  radiation  to  the  total  heat  transfer.  The  relative 
increase  in  thermal  conductivity  with  temperature,  however,  is  well 
below  the  T3  dependence  expected  if  radiation  were  the  only  con¬ 
tributing  mechanism  of  heat  conduction.  This  indicates  that  in 


THICKNESS  (mm) 


Fig.  6.  Normalized  thermal  conductivity  due  to  radia¬ 
tion  between  the  surfaces  of  samples  of 
various  diameters  as  a  function  of  the 
sample  thickness. 


spite  of  the  observed  specimen  size  effect  and  the  positive  tempera¬ 
ture  dependence,  phonon  transport  still  makes  the  primary  contri¬ 
bution  to  the  thermal  conductivity.  Of  course,  the  temperature 
and  spectral  dependence  of  the  absorption  coefficient  may  also  play 
a  role.  However,  at  this  time  no  data  for  these  effects  appear  to 
be  available  for  these  composites. 


The  data  presented  in  Figure  7  indicate  that  the  curves  for 
the  temperature  dependence  of  the  thermal  conductivity  for  the 
two  values  of  specimen  thickness  are  nearly  parallel  for  all  three 
materials.  This  suggests  that  the  relative  effect  of  specimen  thick¬ 
ness  on  thermal  conductivity  decreases  with  increasing  temperature. 
This  in  turn  suggests  that  heat  transfer  by  photon  emission  and  ab¬ 
sorption  makes  a  relatively  larger  contribution  to  the  total  thermal 
conductivity  with  increasing  temperature. 

It  is  interesting  to  note  that  compared  to  the  present  materials 
the  contribution  of  radiation  to  the  thermal  conductivity  of  such 
highly  transparent  materials  as  quartz  and  sapphire  appears  to  be 


m 


m 


I 


much  higher,  as  observed  by  Howlett  (10)  .  This  suggests  that  the 
suppression  of  radiative  heat  transfer  by  photon  scattering  at  op¬ 
tical  discontinuities  could  play  a  significant  role  in  the  conduction 
of  heat  in  the  materials  investigated  in  this  study. 

In  summary,  it  was  found  that  radiative  heat  transfer  between 
the  specimen  surfaces  contributes  to  the  overall  thermal  conducti¬ 
vity  of  a  silicon  carbide  fiber-reinforced  lithium  aluminosilicate 
glass-ceramic.  The  dependence  of  thermal  conductivity  on  specimen 
thickness  suggests  that  the  radiative  heat  transfer  is  partially 
governed  by  the  view  factor  between  the  carbon  coated  laser-flash 
specimens . 
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ABSTRACT 

The  thermal  diffusivity  of  hot  pressed  silicon  carbide-boron  nitride 
composites  is  shown  to  exhibit  a  high  degree  of  anisotropy  with  respect  to  the 
hot  pressing  direction  due  to  preferred  orientation  of  the  boron  nitride 
phase.  This  effect  is  particularly  pronounced  at  the  higher  boron  nitride 
contents  at  both  room  and  elevated  temperature. 

INTRODUCTION 

Single  crystals  of  most  materials  exhibit  anisotropy  i  many  of  their 
tensor  properties,  including  the  electrical  and  thermal  conductivity,  and  the 
elasticity.^  While  such  anisotropy  is  most  pronounced  with  the  crystal 
structures  of  materials  such  as  mica,  graphite  and  boron  nitride,  it  is  also 
exhibited  by  polycrystalline  aggregates  of  anisotropic  single  crystals  if  the 
grains  have  a  preferred  orientation  (i.e.,  texture)  introduced  during  or  after 
manufacture.  In  composites,  the  individual  components  frequently  are 
oriented  preferentially.  This,  in  addition  to  other  variables  such  as  the 
distribution  of  the  phases,  also  contributes  to  the  anisotropy  of  the  continuum 
properties  of  the  composites.  This  was  recently  demonstrated  by  experimental 
data”  for  the  room  temperature  thermal  diffusivity  of  hot-pressed  composites 


of  silicon  nitride  with  a  dispersed  phase  of  boron  nitride.  At  the  highest 


concentration  of  BN  investigated  (30  wt.%)  these  composites  exhibited  a  factor 
of  two  difference  in  the  thermal  diffusivity  parallel  and  perpendicular  to  the 
hot-pressing  direction.  X-ray  analysis  confirmed  that  this  anisotropy 
resulted  primarily  from  the  preferred  orientation  of  the  boron  nitride  phase, 
probably  introduced  during  compaction  of  the  powders  and  facilitated  by  the 
plate-like  shape  of  the  individual  BN  particles.  The  purpose  of  this 
communication  is  to  present  further  confirming  data  for  the  effect  of  the 
preferred  orientation  of  boron  nitride  on  the  room  temperature  thermal 
diffusivity  of  hot-pressed  composites  of  silicon  carbide  and  boron  nitride  over 
the  total  range  of  BN  content. 

EXPERIMENTAL 

A  range  of  SiC-BN  compositions  were  prepared  using  a  sinterable  f$-SiC 

*  4 

powder  and  a  BN  powder  with  a  particle  size  of  15  ym.  The  powders 
were  mixed  in  a  plastic  ball  mill  with  silicon  nitride  balls  for  4  hours  in 
isopropyl  alcohol.  After  drying,  the  powders  were  pelletized  and  then 

A  A  0 

uniaxially  hot-pressed  in  vacuum  at  2100  C  for  1  hour  at  a  pressure  of 

♦  + 

35  MPa  in  graphite  dies  lined  with  graphite  foil.  After  the  specimens  were 
removed  from  the  die,  the  surfaces  were  ground  and  the  density  was 
determined  from  the  dimensions  (=2.5  cm  diameter  by  1  cm  thick)  and  the 
mass . 

X-ray  diffraction  analysis  of  the  hot-pressed  blanks  revealed  the 
presence  of  several  SiC  polytypes,  including  21R,  15R,  6H,  4H,  and  3C. 


*  Lot  S  1626,  beta  silicon  carbide,  Herman  C.  Starck,  Goslar,  Federal 
Republic  of  Germany. 

♦  B-27,  Research  Organic/Inorganic  Chemical  Corp.,  Belleville,  NJ. 

**  Model  300,  Vacuum  Industries  Div.,  GCA  Corp.,  Chicago,  IL. 

*♦  Carbon  Products  Division,  Union  Carbide  Corp.,  Chicago,  IL. 


With  this  multiphase  assemblage  it  was  not  possible  to  quantify  the  amount  of 
each  phase  present,  but  the  SiC  patterns  were  similar  regardless  of  the  BN 
content.  Boron  nitride  was  always  present  as  the  hexagonal  graphite 
structure,  and  no  evidence  of  reaction  or  solid  solution  between  the  SiC  and 
BN  was  observed. 

The  effect  of  the  concentration  and  the  orientation  of  the  BN  on  the  heat 

transfer  behavior  of  the  composites  was  determined  by  measurements  of  the 

3 

thermal  diffusivity  by  the  flash  method  using  a  glass-Nd  laser  as  the  flash 
source  and  remote  optical  sensing  for  monitoring  the  transient  temperature 
response  of  the  specimen.  At  room  temperature,  the  measurements  were 
conducted  in  air,  whereas  at  elevated  temperatures,  the  specimen  was 
contained  within  a  carbon  resistance  furnace  with  a  protective  nitrogen 
atmosphere.  The  specimens  used  for  these  measurements  were  cut  from  the 
hot-pressed  blanks  and  were  in  the  form  of  square  plates  =2  mm  thick  with  a 
diagonal  of  =12  mm.  A  number  of  specimens  were  cut  from  each  blank  in 
order  to  obtain  measurements  of  the  thermal  diffusivity  in  directions  pa.aiiel 
and  perpendicular  to  the  hot-pressing  direction. 

RESULTS  AND  DISCUSSION 

Table  I  lists  the  composition,  density,  and  porosity  of  each  composite. 

The  amount  of  porosity  was  determined  from  the  measured  density  and  the 

3 

theoretical  density  based  on  values  of  2.28  g/cm  for  the  density  of  boron 

3 

nitride  and  3.21  g/cm  for  that  of  silicon  carbide.  As  indicated  by  the 
table,  the  pore  content  generally  increases  with  increasing  BN  content. 

Figure  1  shows  the  experimental  data  for  the  thermal  diffusivity  at  room 
temperature  both  parallel  as  well  as  perpendicular  to  the  hot-pressing 
direction.  The  scatter  in  the  data  at  a  given  BN  content  and  direction  of 


heat  flow  indicates  minor  variations  from  specimen  to  specimen  cut  from  the 

hot-pressed  blank.  For  heat  flow  perpendicular  to  the  hot-pressing 

direction,  the  thermal  diffusivity  is  nearly  independent  of  BN  content.  This 

is  expected  since  the  measured  thermal  diffusivity  of  the  BN  perpendicular  to 

the  hot-pressing  direction  is  nearly  the  same  as  that  of  the  silicon  carbide. 

9 

Note  however  that  this  value  for  the  SiC  (=0.4  cm  /s)  is  well  below  the 

4 

value  of  thermal  diffusivity  observed  for  other  samples  of  SiC,  which  can 

2 

be  as  high  as  0.8  cm  /s.  The  lower  value  is  most  likely  a  reflection  of  the 

nature  and  amount  of  impurities,  to  which  the  diffusivity  of  SiC  is  especially 
4 

sensitive. 

The  thermal  diffusivity  of  the  composites  parallel  to  the  hot-pressing 

direction  decreases  rapidly  with  increasing  BN  content,  resulting  in  a 

maximum  anisotropy  ratio  of  about  four  between  the  thermal  diffusivity  in  the 

two  directions.  Since  the  thermal  diffusivity  along  the  basal  plane  of  boron 

nitride  is  greater  than  that  perpendicular  to  the  basal  plane,  the  anisotropy 

observed  in  the  composites  suggests  that  the  BN  is  preferentially  oriented 

with  the  basal  plane  perpendicular  to  the  direction  of  hot-pressing.  For  the 

2 

silicon  nitride-boron  nitride  composites  investigated  in  an  earlier  study, 

this  preferred  orientation  was  verified  by  x-ray  analysis.  In  the  present 

samples  the  orientation  of  the  boron  nitride  can  be  seen  in  the  scanning 

electron  fractographs  of  the  specimen  containing  30  at.o  BN  shown  in  figures 

2a  and  2b  for  directions  parallel  and  perpendicular  to  the  hot-pressing 

5  6 

direction.  In  two  separate  investigations  '  of  boron  nitride  single 
crystals,  the  ratio  of  the  thermal  conductivity  at  300  K  between  the  direction 
perpendicular  to  the  basal  plane  and  that  parallel  to  the  basal  plane  was 
about  forty  and  one  hundred.  This  indicates  that  the  orientation  of  the  BN 
in  the  present  composites,  while  partly  preferential,  is  mostly  random. 


The  values  of  thermal  diffusivity  for  the  silicon  carbide  also  indicate  a 

slight  anisotropy.  Since  other  high  conductivity  ceramics  with  the  wurtzite 

crystal  structure  exhibit  very  little  anisotropy  in  the  thermal 
7  8 

conductivity,  '  the  observed  differences  in  the  silicon  carbide  data  may  be 
the  result  of  a  preferred  orientation  of  certain  particle  morphologies, 
especially  the  hexagonal  polytypes. 

The  thermal  diffusivity  will  also  be  affected  by  the  amount  of  porosity 

9  10 

present.  Theoretical  estimates  by  Raleigh  and  Maxwell  show  that  at 
low  porosity  contents,  the  relative  decrease  in  thermal  conductivity  is  1.5 
times  the  relative  pore  content.  Due  to  the  associated  change  in  density,  the 
corresponding  relative  change  in  thermal  diffusivity  should  be  0.5  times  the 
relative  pore  content.  Assuming  that  this  value  is  applicable  to  the  present 
composites,  the  thermal  diffusivity  was  corrected  for  porosity  and  the  values 
are  listed  in  table  l.  Because  the  correction  is  relatively  small,  the 
conclusions  based  on  the  original  data  remain  unchanged. 

A  comparison  of  the  temperature  dependence  of  the  thermal  diffusivity 
perpendicular  to  the  hot-pressing  direction  with  that  parallel  to  the  hot- 
pressing  direction  is  shown  in  figures  3a  and  3b  for  composites  with  30  and 
80  at.o  BN,  respectively.  These  results  indicate  that  for  both  compositions, 
the  relative  temperature  dependence  of  the  thermal  diffusivity  in  one  direction 
of  heat  flow  is  not  significantly  different  from  that  in  the  other  direction. 
This  also  implies  that  the  orientation  of  the  BN  is  mostly  random. 

The  observed  anisotropy  of  the  thermal  diffusivity,  especially  at  the 
higher  BN  contents,  could  be  of  practical  interest  for  engineering 
applications .  This  is  particularly  so  for  structures  or  components  which,  for 
improved  energy  efficiency,  require  high  thermal  insulating  ability  in  one 
direction  of  heat  flow  and  high  thermal  conduction  in  a  perpendicular 


direction,  thus  reducing  the  possibility  of  thermal  stress  failure  due  to 
spatially  non-uniform  heating. 
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ABSTRACT 

Crystallization  of  the  glassy  grain  boundary  phase  in  a  series  of  Sialon 
ceramics  fabricated  using  a  range  of  hot  pressing  schedules  increased  the 
thermal  diffusivity  at  room  temperature  by  an  average  of  about  10o.  For 
samples  made  by  a  given  hot  pressing  schedule,  the  relative  difference  in  the 
thermal  diffusivity  between  composites  containing  either  a  glassy  or  a 
crystallized  grain  boundary  phase  decreased  with  increasing  temperature. 
This  was  attributed  to  enhanced  phonon  scattering  in  the  crystalline  grain 
boundary  phase. 

INTRODUCTION 

The  densification  during  sintering  of  polycrystalline  ceramics  for  high- 
performance  applications  can  be  enhanced  significantly  by  the  addition  of  a 
sintering  aid.  This  is  particularly  so  for  silicon  nitride-based  ceramics, 
which  do  not  density  readily  without  the  assistance  of  additives.  For  silicon 
nitride,  such  additives  commonly  consist  of  the  oxides  of  magnesium,  yttrium 
or  aluminum.  At  the  fabrication  temperatures,  the  presence  of  these 
additives  results  in  the  formation  of  a  liquid  phase  at  the  grain  boundaries. 
This  liquid  phase  can  act  either  as  a  lubricant  to  promote  particle  transport 


or  as  a  medium  for  material  transport  by  solution  and  re-precipitation  of  the 
matrix  phase.  Upon  cooling,  the  liquid  at  the  grain  boundaries  and  triple 
points  becomes  a  glass. 

The  presence  of  this  grain  boundary  phase,  which  can  vary  from  just  a 
few  percent  to  over  10»,  is  expected  to  have  an  effect  on  the  continuum 
properties  of  the  final  densified  product.  This  effect  is  well  known  for 
mechanical  behavior,  especially  at  temperatures  sufficiently  high  that  the 
grain  boundary  phase  can  undergo  significant  viscous  deformation. 

The  effect  of  the  grain  boundary  phase  on  the  continuum  properties  of 
polycrystalline  aggregates  also  is  expected  to  be  a  function  of  the  degree  of 
crystallinity.  In  terms  of  the  objectives  of  this  study,  it  should  be  noted 
that  due  to  greatly  enhanced  phonon  scattering,  glassy(amorphous)  materials 
generally  exhibit  a  much  lower  thermal  conductivity  and  diffusivity  than  their 
crystalline  counterparts  of  identical  composition!  1 -2) .  A  polycrystalline 
ceramic  with  a  grain  boundary  phase  can  be  considered  to  be  a  composite. 
Theory(3)  indicates  that  the  thermal  conductivity  of  composites  is  a  function 
of  the  thermal  conductivity,  volume  fraction,  and  distribution  of  each 
component.  Continuous  phases  are  particularly  effective  in  affecting  the 
thermal  conductivity  of  composites.  For  this  reason,  the  presence  of  a  grain 
boundary  phase  and  its  degree  of  crystallinity  are  expected  to  have  a 
significant  effect  on  the  thermal  transport  properties  of  polycrystalline 
aggregates.  Experimental  data  are  presented  in  support  of  this  hypothesis. 

EXPERIMENTAL 

The  materials  for  this  study  consisted  of  a  $-Sialon  (ft-SigN^  solid 
solution  containing  5  eq.°  Al,  i.e.,  B-SL  Al  O  N„  where  x=0.4) 

O-  X  X  x  o~  X 

with  a  sintering  aid  corresponding  to  the  composition  of  the  mineral  garnet. 


1.6.  , 


Y3AI5°12  (4).  The  samples  were  prepared  by  hot-pressing 

* *  **  + 
appropriate  mixtures  of  powders  of  Si^N^  ,  AIN  ,  A^O^  ,  and 

Y7Og  in  boron  nitride  coated  graphite  dies  at  a  pressure  of  -20  MPa 

in  a  nitrogen  atmosphere  at  temperatures  ranging  from  1550  to  1690°C  for 

time  periods  ranging  from  15  to  120  minutes.  The  densities  achieved  were  at 

least  99.5°  of  theoretical.  The  amount  of  densification  aid  added  was 

approximately  10  wt.°.  Following  hot-pressing,  the  samples  were  annealed  at 

1350°C  for  20  hours,  which  was  found  to  result  in  nearly  complete 

crystallization  of  the  garnet  grain  boundary  phase.  Details  of  the  x-ray  and 

TEM  studies  of  the  crystallization  process  were  presented  in  an  interim 

technical  report(4)  intended  for  formal  publication  at  a  later  date.  In  the  "as 

hot-pressed"  samples,  B-Sialon  was  the  only  crystalline  phase,  whereas  garnet 

was  the  only  second  crystalline  phase  in  the  annealed  specimens. 

The  effect  of  the  crystallization  of  the  grain  boundary  phase  on  thermal 
transport  behavior  was  determined  by  measurement  of  the  thermal  diffusivit/ 
by  the  flash  method(5)  using  a  glass-Nd  laser  as  the  flash  source.  The 
specimens  for  these  measurements  consisted  of  square  plates  about  6  mm  on  a 
side  by  about  1.5  mm  thick.  Measurement  of  the  thermal  diffusivity  at  room 
temperature  was  conducted  in  air;  at  elevated  temperatures  the  specimens 
were  contained  in  a  carbon  furnace  with  a  protective  nitrogen  atmosphere. 

RESULTS  AND  DISCUSSION 

Figure  la  compares  the  data  for  the  thermal  diffusivity  at  room 
temperature  prior  to  and  following  the  annealing  of  samples  hot-pressed  for  30 
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min  at  various  temperatures.  Similarly,  Figure  1b  compares  the 
corresponding  data  for  specimens  hot-pressed  at  1730°C  for  various  periods 
of  time.  Both  sets  of  data  show  an  increase  in  the  thermal  diffusivity  of 
about  10°  after  annealing,  regardless  of  the  hot-pressing  conditions. 

Figure  2  compares  the  temperature  dependence  of  the  thermal  diffusivity 
for  Sialon  samples  with  a  glassy  and  a  crystalline  grain  boundary  phase. 
Both  samples  were  hot  pressed  for  30  minutes  at  1690°C.  The  relative 
difference  in  the  thermal  diffusivity  decreases  with  increasing  temperature. 
This  is  expected  since  enhanced  phonon  scattering  generally  produces  a 
greater  relative  decrease  in  the  thermal  conductivity  and  diffusivity  of 
crystalline  materials  than  for  glassy  materials.  The  data  in  figure  2  were 
obtained  up  to  1000°C  in  order  to  avoid  crystallization  of  the  grain 
boundary  phase  at  higher  temperatures.  Such  crystallization  would  lead  to  a 
permanent  increase  in  the  thermal  diffusivity  of  the  as  hot  pressed"  sample 
if  measurements  had  been  conducted  to  temperatures  as  high  as  1200  or 
1300°C.  Although  the  temperature  dependence  of  the  thermal  diffusivity 
was  measured  for  only  one  hot  pressing  condition,  it  is  anticipated  that  the 
samples  for  the  other  hot  pressing  conditions  would  exhibit  a  behavior  similar 
to  that  shown  in  figure  2. 
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TEMPERATURE  (*C) 


Temperature  dependence  of  the  thermal  dif- 
fusivity  of  B-Sialon/garnet  hot-pressed  at 
1690°C  for  30  min  compared  to  B-Sialon/garnet 
hot-pressed  at  the  same  conditions  and  then 
annealed  at  1350°C  for  20  hours.  The  solid 
lines  represent  linear  least  squares  fits  of 
the  measured  thermal  diffusivity  to  the 
reciprocal  absolute  temperature. 
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EFFECT  OF  DROP-HEIGHT  ON  THE  CRITICAL  TEMPERATURE  DIFFERENCE 
(ATC)  FOR  BRITTLE  CERAMICS  SUBJECTED  TO  THERMAL  SHOCK  BY  QUENCHING 

Jitandra  P.  Singh,  Gunther  Ziegler,  D.  P.  H.  Hasselman 

ABSTRACT 

Values  for  the  critical  quenching  temperature  difference  (AT  )  required 
for  the  initiation  of  thermal  stress  failure  of  soda-lime-silica  glass, 
polycrystall ine  alumina,  mullite  and  reaction-sintered  silicon  nitride 
subjected  to  thermal  shock  by  quenching  into  water  were  found  to  be  indepen¬ 
dent  of  the  height  from  which  the  specimens  were  dropped.  These  results 
suggest  that  within  the  range  of  usual  laboratory  practice  no  close  control 
over  the  drop  height  needs  to  be  exercised.  The  relevance  of  the  present 
observations  in  terms  of  the  mechanism  of  heat  transfer  responsible  for 
thermal  stress  failure  in  a  water  quench  is  briefly  discussed. 

The  quench  test  is  a  favorite  method  for  the  measurement  of  the  thermal 
stress  resistance  of  brittle  ceramics1-^.  Usually,  the  method  consists  of 
dropping  appropriate  specimens  under  free-fall  conditions  into  a  fluid  at  or 
near  room  temperature  from  an  oven  or  furnace  at  elevated  temperature.  The 
critical  quenching  temperature  difference  ( ATC )  required  to  initiate  fracture 
as  evidenced  by  extensive  crack  formation  is  a  measure  of  the  thermal  stress 
resistance  of  the  particular  ceramic  being  tested.  Unfortunately,  the 
quantitative  interpretation  of  the  data  for  ATC  is  handicapped  by  uncertainti 


in  the  relative  contributions  of  heat  transfer  by  convective  processes, 
nucleate  boiling  and  steam-film  formation. 

It  is  expected  that  these  heat  transfer  mechanisms  could  be  affected 
by  the  velocity  at  which  the  specimens  hit  the  fluid  surface,  which,  in 
turn  is  a  function  of  the  height  over  which  the  specimens  are  dropped. 

Very  few,  if  any,  reports  on  the  results  for  ATC  obtained  by  water 
quenching  give  values  for  the  drop-height.  For  this  reason  the  purpose 
of  this  study  was  to  obtain  experimental  data  for  this  effect  which  should 
also  throw  a  light  on  the  dominant  mechanism  of  heat  transfer. 

Four  different  ceramic  materials  were  selected,  as  listed  in  Table  1, 
together  with  the  specimen  size,  geometry  and  values  for  the  relevant 
material  properties,  corresponding  to  the  approximate  mean  value  at  the 
initial  temperature  of  the  specimens  at  AT  =  ATc.  Five  specimens  were 
quenched  simultaneously  from  a  range  of  values  of  height  and  temperature 
difference  (aT)  from  an  electrically  heated  laboratory  oven  into  a  water 
bath  held  at  33°C.  Each  specimen  was  tested  only  once  in  order  to 
eliminate  the  possible  effect  of  fatigue  on  aTc.  The  circular  specimens  hit 
the  surface  in  a  horizontal  position,  whereas  the  square  specimens  of 
silicon  nitride,  tested  in  another  laboratory,  hit  the  surface  in  a  vertical 
position.  Following  quenching,  the  strength  retained  by  the  specimens  was 
measured  in  4-point  bending  at  a  cross-head  speed  of  0.05  cm/min.  with  an 
inner  and  outer  span  for  the  circular  specimens  of  1.6  and  5.0  cm,  respectively 
For  the  silicon  nitride,  the  corresponding  values  were  2  and  4  cm,  respectively 
The  value  of  ATC  was  defined  by  the  range  and  mean  of  the  values  of  quenching 
temperature  difference  (AT)  at  which  1  or  5  of  the  five  specimens  tested 
exhibited  a  significant  strength  loss.  At  values  of  AT  near  ATC,  AT  was  varied 
over  10°C  intervals  in  order  to  define  ATC  as  closely  as  possible. 


The  experimental  results  are  shown  in  Fig.  1.  For  all  four  materials 

no  significant  difference  exists  between  the  value  of  ATC  over  the  total 
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range  of  values  of  height  used.  In  view  of  the  square  root  dependence 
of  specimen  velocity  on  drop-height,  the  range  of  velocities  involved  a 
factor  in  excess  of  six.  During  their  transit  from  the  oven  to  the  water 
bath,  the  specimens  should  lose  some  heat  to  the  air  by  convection  or  by 
radiation.  This  effect,  however,  would  serve  to  increase  ATC  with  in¬ 
creasing  height,  other  effects  being  absent.  Due  to  the  very  low  heat 
transfer  coefficient,  the  loss  of  heat  to  the  air  by  convection  (8)  should 
have  a  negligible  effect  on  ATC.  The  slight  apparent  increase  in  ATc 
with  increasing  height  for  the  mullite  and  silicon  nitride  specimens  which 
require  the  higher  values  for  the  initial  specimen  temperature  possibly 
may  be  attributed  to  radiative  heat  loss  during  the  transit  to  water  bath. 
However,  from  the  dependence  of  radiative  heat  transfer  on  absolute  tempera¬ 
ture  it  is  easily  calculated  that  significant  changes  in  ATC  due  to  radia¬ 
tive  heat  losses  would  occur  only  for  initial  specimen  temperatures  well 
in  excess  of  1000°C  for  the  range  of  values  of  drop-height  used  in  this 
study.  Regardless  of  the  details,  these  experimental  results  must  be 
considered  as  rather  gratifying  as  they  suggest  that  for  the  range  of 
quenching  conditions  for  a  water  bath  generally  encountered  in  practice, 
close  control  over  the  drop-height  need  not  be  exercised. 

The  implications  of  the  present  results  for  the  mechanisms  of  heat 
transfer  can  be  speculated  on  briefly.  The  independence  of  ATC  from  drop- 

★ 

It  is  not  the  purpose  of  this  study  to  compare  the  values  of  ATC  for  the  four 
materials  because  of  the  differences  in  geometry  involved.  It  should  be 
noted,  however,  that  the  relative  order  of  ATC  is  the  same  as  the  order  of 
the  parameters  (8,9):  af(l-v)/otE  and  o^(l-v)K/aE. 


height  suggests  that  the  mechanism  of  heat  transfer  responsible  for  the 
thermal  fracture  is  not  governed  by  the  velocity  at  which  the  specimens 
hit  the  surface  and  the  corresponding  forced  convective  heat  transfer 
coefficient.  From  heat  transfer  theory16,  this  latter  quantity  can  be 
calculated  to  vary  by  a  factor  of  about  two  for  the  range  of  drop-heights 
used.  This,  in  turn  suggests  that  the  transfer  of  heat  is  controlled  by 
either  forced  convection  corresponding  to  the  "terminal"  velocity  at  which 
the  specimen  falls  through  the  water  (which  is  independent  of  the  drop- 
velocity)  or  by  nucleate  boiling  or  steam-film  formation.  In  fact,  for  the 
circular  samples  of  aluminum  oxide  for  which  ATC  is  not  strongly  affected 
by  sub-critical  crack  growth  as  demonstrated  by  computer  simulation  ,  the 
value  of  AT„  calculated  from  the  heat  transfer  coefficient  for  the  free-fall 
velocity6’14  and  the  appropriate  thermo-elastic  equation,  agrees  very  well 
with  the  observed  value.  Such  agreement,  however,  must  be  concerned  fortuitous 
as  for  the  initial  specimen  temperature  nucleate  boiling  and  steam-film  forma¬ 
tion  are  expected  to  play  a  significant  role.  In  fact,  the  effective  heat 
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transfer  coefficient  in  the  presence  of  steam- film  formation  ’  is  of 
approximately  the  same  magnitude  as  the  one  calculated  for  free-fall  conditions. 
As  indicated  by  experimental  evidence  ’  ,  the  rate  of  heat  transfer  by 

nucleate  boiling  and  steam- film  also  is  relatively  independent  of  the  relative 
velocity  between  the  water  and  the  hot  surface.  For  this  reason,  ATC  for 
heat  transfer  by  these  latter  mechanisms  would  also  be  independent  of  specimen 
velocity  (i.e.,  drop-height). 
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Fig.  1.  Effect  of  drop-height  on  the  critical  quenching 
temperature  difference  ATC  required  to  initiate 
thermal  stress  fracture  in  four  ceramic  materials 


TABLE  1 .  CERAMIC  MATERIALS 


imina,  McDanel  Porcelain,  Beaver  Fal 
McOanel  Porcelain,  Beaver  Falls,  P/ 
;  DFVLR,  Cologne,  FRG. 
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ABSTRACT 


A  study  was  conducted  of  the  effect  of  the  magnitude  of  the  thermal 
diffusivity  on  the  thermal  stresses  in  a  ceramic  plate  subjected  to 
thermal  cycling  by  changes  in  ambient  temperature  involving  convective 
heat  transfer,  at  frequency  sufficiently  high  that  thermal  equilibrium 
(i.e.,  temperature  uniformity)  is  not  achieved  within  each  cycle.  The 
results  obtained  indicate  that  at  any  specific  value  of  the  Biot  number 
the  magnitude  of  maximum  thermal  stress  within  any  cycle  decreases  with 
decreasing  thermal  diffusivity  and  increases  with  increasing  plate 
thickness.  The  approximate  expression  for  the  thermal  stress  resistance 
parameter  appropriate  for  very  high  cycling  frequencies  and  low  values 
of  the  Biot  number  indicates  that  good  thermal  stress  resistance  under 
these  conditions  requires  a  high  value  for  the  specific  heat  per  unit 
volume. 
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I.  INTRODUCTION 
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Ceramic  materials  for  structural  purposes  frequently  exhibit  high 
values  of  Young’s  modulus  and  coefficient  of  thermal  expansion  and 
relatively  low  values  of  tensile  strength.  This  combination  of  properties, 
together  with  a  high  degree  of  bri ttleness,  make  structural  ceramics 
highly  susceptible  to  catastrophic  failure  under  the  influence  of 
steady-state  or  transient  thermal  stresses  usually  encountered  during 
normal  use  at  elevated  temperatures.  For  this  reason,  it  is  imperative 
for  design  purposes  that  the  variables  which  influence  the  magnitude  of 
thermal  stress  and  the  resulting  failure  are  well  understood.  For 
this  purpose  a  number  of  analytical  studies  have  appeared  in  the  ceramics 


I i terature. 


These  studies  have  permitted  the  derivation  of  so-called 


thermal  stress  "resistance  parameters"  or  "f igures-of-meri t"  which 
permit  the  selection  of  the  ceramic  from  among  a  number  of  candidate 
materials  with  optimum  thermal  stress  resistance  for  a  specific  thermal 
environment  and  failure  mode. 

The  derivations  of  the  thermal  stress  resistance  parameters  for 
ceramic  materials  presented  to  date,without  exception,  have  considered 
steady-state  heat  flow,  or  in  the  case  of  transient  heat  transfer, 
situations  in  which  initial  and/or  final  thermal  equilibrium  is  attained. 
In  practice,  however,  situations  involving  thermal  shock  are  encountered 
in  which  the  ceramic  is  not  in  a  state  of  thermal  equilibrium.  For 
instance,  this  is  the  case  for  ceramic  components  of  internal  combustion 
or  turbine  engines,  which  involve  thermal  cycling  of  such  frequency  that 
no  thermal  equilibrium  is  attained  within  the  time-duration  of  each  cycle. 
It  is  anticipated  that  for  these  conditions,  the  thermal  stress 


WOT 
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resistance  parameters  will  differ  from  those  appropriate  to  steady-state 
heat  flow  or  those  which  involve  thermal  equilibrium  following  or 
preceding  transient  heat  flow. 

It  is  the  purpose  of  this  study  to  develop  an  analysis  of  the 
material  properties  which  affect  thermal  stress  failure  of  ceramic 
components  subjected  to  rapid  thermal  cycling  and  to  establish  the 
appropriate  thermal  stress  resistance  parameters. 

I  I .  THEORY 


The  geometry  selected  for  this  study  consists  of  an  infinite  flat 
plate  of  thickness  2a  located  in  the  y-z  plane  with  -a  <  x  <  a.  The 
plate  is  exposed  to  an  ambient  environment,  with  heat  transfer  to  and 
from  the  plate  taking  place  by  Newtonian  convection.  The  plate  is 
subjected  to  rapid  thermal  cycling  by  changes  in  the  temperature  of 
the  ambient  alternating  instantaneously  between  upper  and  lower 
temperatures  T^  and  TQ,  respectively,  for  equal  duration  with  a  total 

period  of  At  illustrated  schematically  in  Fig.  1. 
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A  solution  of  the  heat  conduction  equation; 

k32T/9x2  =  9T/3t  da) 


is  required  for  each  time  interval; 


At  <  t  <  |  At 


(1b) 


where  k  is  the  thermal  diffusivity  and  m  refers  to  the  m"l’h  half-cycle 
with  m  =  I,  2,  3  ....  », subject  to  the  boundary  condition  of  symmetry: 

9T/9x  =  0  at  x  =  0  (2) 
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with  Newtonian  convection  at  the  surfaces  x  =  -  a, 


*(!;]_'  "£<».+>  -  V+>] 

X—  3 


(3) 


where  k  is  the  thermal  conductivity  and  h  is  the  heat  transfer  coefficient 
(assumed  to  be  the  same  for  heating  and  cooling)  and  Tm  =  or  T  ,  for 
the  heating  or  cooling  phase  of  the  cycle,  respectively. 

The  initial  conditions  are  different  for  each  time  interval, 
but,  in  general,  the  initial  temperature  distribution  for  each  time 


interval  will  be  identical  to  the  distribution  at  the  end  of  the  pre¬ 


ceding  interval.  The  solution  is  found  by  standard  techniques 
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th 


the  m  half-cycle,  we  find 


°°  -KA^t 

T ( x,t)  =  T  +  l  C.  e  1  cos(A„x) 
m  dm  £, 


(4) 


where 


T  = 
m 


To,  m  even; 


T^ ,  m  odd 


The  eigenvalues  A?  are  solutions  of  the  transcendental  equation: 


A^a  tan  (A^a)  =  3 


(5) 


with  Biot  number 


8  -  — 
6  ~  k 


(6) 


th 


The  expansion  coeff'cients  Cf  for  the  m  time  interval  are 


2(T,-T  )$  sin  A„a 


f  'o 


Urn 


A  a(6  +  sin  A  a) 


m=  1 


J.  n  c-n'-J’e 

LJ=0 


. ,  -jrA  "At/2 

it  m 


m  >  1 


(7) 


The  thermal  stress  is  computed  from  the  temperature  distribution 
accord ing  to: ^ 


a  (x,t)  = 
yz  1-v 


.  a  ,  a 

-T(x,t)  +  2^  /  T(x',t)dx'  +  ~  /  T(x',t)x'dx' 
-a  2aJ  -a 


,  ( 


th 


to  yield,  for  the  m  half-cycle, 

-KAtfsin  A, 


o  (x,t)  =  7^—  y  C 
yz  1-v  k  , 


yz 


t- 1 


Urn 


V 


cos  A^a 


where  C  is  given  by  Eq.  (7)  and  a  is  the  thermal  expansion  coefficient 
E  is  Young’s  modulus,  and  v  is  Poisson's  ratio. 

The  temperature  was  calculated  in  terms  of  the  non-dimensional 
temperature,  T*: 

T*  =  (T  -  T  )/(T,  -  T  )  (1( 

o  f  o 

and  the  stress  in  terms  of  the  non-dimensional  stress,  a*: 


a*  =  a( 1-v)/aE(T,  -  T  ) 
f  o 

the  non-dimensional  time,  t*: 


(1 


t*  =  xt/a^ 


(i: 


and  the  non-dimensional  period  for  each  cycle: 


For  convenience,  these  non-dimensional  quantities  will  be  referred 


to  simply  as  temperature,  stress,  time  and  period,  respectively. 

II.  NUMERICAL  RESULTS  AND  DISCUSSION 

Figures  2a  and  2b,  for  6  =  2.5,  show  the  distribution  of  the 
temperatures  and  stresses  encountered  during  the  first  cycle  for  a 
range  of  values  of  time,  t*,  for  a  value  of  the  cycling  period 
At*  =  50.  The  correspond i ng  value  of  frequency  is  sufficiently  low 
that  during  the  heating  and  cooling  phases  of  the  cycle  thermal  equi¬ 
librium  (i.e.,  temperature  uniformity  at  T^  and  T  ,  respectively)  is 

achieved.  The  numerical  values  for  the  temperature  and  stresses  agree 
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with  other  literature  data  ,  which  supports  the  validity  of  the 
theory  presented  above.  For  purposes  of  further  discussion  it  should  be 
noted  that  the  data  shown  in  Fig.  2,  because  of  the  large  value  of  the  perioi 
At*  (i.e.,  low  frequency),  will  not  vary  from  cycle-to-cycle,  regardless 
of  the  number  of  cycles  to  which  the  plate  is  subjected. 

Figures  3a  and  3b,  for  a  value  of  B  =  2.5,  show  the  temperature 
and  stress  distribution  for  a  range  of  values  of  t*  during  the  first 
cycle,  with  a  period  At*  =  1.0.  This  value  of  period  is  sufficiently 
brief  that  within  the  cycle  thermal  equilibrium  is  not  achieved.  This 
is  clearly  evident  from  the  data  presented  because  the  value  of  tempera¬ 
ture  (T*)  does  not  attain  the  value  of  unity. 

In  general,  after  a  sufficient  number  of  cycles,  a  condition  will 
be  achieved  to  which,  for  a  specific  value  of  time  for  any  cycle,  no 
cycle-to-cycle  variation  in  temperature  or  thermal  stress  is  found. 

For  convenience,  this  condition  will  be  referred  to  as  cycle-to-cycle 
equi I ibrium. 


Figures  4a  and  4b  show  the  distribution  of  the  temperature  and 
stress,  respectively,  during  a  cycle  for  which  such  cycle-to-cycle 
equilibrium  was  attained  for  a  vlaue  of  3  =  2.5  and  At*  =  1.0,  identical 
to  the  values  of  Fig.  3.  Comparison  of  the  data  given  in  Fig.  4  with 
those  in  Fig.  3  shows  that  the  magnitude  of  the  maximum  compressive 
stresses  during  the  equilibrium  cycle  are  below  those  encountered  during 
the  first  cycle,  and  conversely  for  the  tensile  stress.  This  effect 
arises  because  during  heating  in  the  first  cycle  the  value  of  tensile 
stress  at  the  center  of  the  plate  rises  slower  than  the  compressive 
stress  in  the  surface.  For  this  particular  value  of  3  if  the  value  of 
At*  is  smaller  than  0.2,  the  cooling  phase  will  be  initiated  well  before 
tensile  stress  could  reach  a  peak  value,  if  heating  had  been  continued. 

It  should  be  noted  that  for  the  thermal  cycling  conditions  for  the 
data  in  Fig.  4,  the  stress  at  any  instant  of  time  may  exhibit  a  maximum 
at  a  position  in  the  plate  not  necessarily  at  the  center  or  surface. 
However,  it  was  found  that  the  maximum  value  of  stress  encountered  during 
any  complete  cycle  occurs  in  the  surface,  except  for  tensile  stress  in 
the  first  complete  cycle  which  for  values  of  At*  between  0.1  and  1.0  may 
occur  in  the  center. 

Figure  5  shows  the  envelope  of  the  absolute  value  of  the  maximum 

values  of  the  tensile  and  compressive  stresses  from  the  first  cycle 

.  # 

until  cycle-to-cycle  equilibrium  is  approached  for  values  of  At 
and  3  identical  to  those  for  Figs.  3  and  4.  As  can  be  seen  for  the 
latter  condition  the  magnitudes  of  the  maximum  tensile  and  compressive 
stresses  are  identical.  This  is  expected  because  at  cycle-to-cycle 
equilibrium  the  spatial-  and  time-averaged  temperature  of  the  plate 
lies  mid-way  between  the  upper  and  lower  ambient  temperatures  between 
which  the  plate  is  being  cycled. 
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Figures  6a  and  6b  give  the  maximum  values  of  the  compressive  and 
tensile  thermal  stresses  at  the  surface  and  center  of  the  plate,  respec 
tively,  during  the  heating  phase  of  the  first  cycle  as  a  function  of 
the  cycling  period.  At*,  and  a  range  of  value  of  the  Biot  number,  3. 

In  general,  as  found  earlier  for  a  specific  case  (Fig.  3),  the  absolute 
values  of  the  compressive  stresses  exceed  the  corresponding  values  of 
the  tensile  stresses. 

Figures  7a  and  7b  show  the  tensile  and  compressive  stresses, 
respectively,  in  the  surface  of  the  plate  during  the  cooling  phase  of 
the  first  cycle  as  a  function  of  the  cycling  period.  At*,  and  a  range 
of  values  ofthe  Biot  number.  Of  interest  to  note  is  that  at  the 
lowest  values  of  the  cycling  period  (i.e.,  high  cycling  frequency)  the 
stresses  in  the  surface  are  compressive  throughout  the  cycle. 

The  opposite  effect  can  be  seen  in  Figs.  8a  and  8b,  in  which  we 
give  the  values  for  the  maximum  compressive  and  tensile  stresses, 
respectively,  at  the  center  of  the  plate,  during  the  cooling  phase  of 
the  first  cycle,  as  a  function  of  the  cycling  period  At*  and  for  a 
range  of  values  of  Biot  number.  In  this  case,  at  the  lower  values  of 
cycling  period  (high  frequency)  the  center  of  the  plate  remains  in  a 
state  of  tension  throughout  the  first  cycle. 

The  values  for  the  maximum  tensile  (or  compressive)  stress  in  the 
surface  and  center  of  the  plate,  at  the  condition  of  rye  I p-to-cyr l r 
equilibrium  as  a  function  of  the  period  At*  and  a  range  of  values  fur 
the  Biot  number  are  shown  in  Figs.  9a  and  9b  ,  respectively.  The 
stresses  in  the  surface  are  significantly  higher  than  those  at  the 


center. 


In  terms  of  the  objectives  of  this  study,  it  is  of  interest  to 
note  that,  as  indicated  by  the  data  presented  in  Figs.  6,  7,  8  and  9, 
at  least  qualitatively  the  magnitude  of  maximum  stress  within  any  cycle 
decreases  with  decreasing  cycling  period  At*,  or  in  other  words  the 
maximum  stress  decreases  with  increasing  cycling  frequency.  Since  At* 
is  defined  as  <At/a  ,  where  At  is  the  period  in  real  time,  this  con¬ 
clusion  implies  that  for  any  given  value  of  At  and  a,  low  thermal 
stress  can  be  achieved  by  choosing  materials  with  a  value  of  thermal 
diffusivity  as  low  as  possible,  other  variables  being  constant.  This 
conclusion  is  in  direct  contrast  with  the  general  concept  that  for  good 
thermal  stress  resistance,  high  thermal  diffusivity  is  desirable.  Of 
course,  high  thermal  diffusivity  is  appropriate  for  heating  or  cooling 
conditions  involving  initial  and/or  final  thermal  equilibrium.  The 
requirement  of  low  thermal  diffusivity  as  concluded  from  the  present 
analysis  is  appropriate  for  thermal  cycling  at  frequencies  for  which 
thermal  equilibrium  is  not  achieved  within  each  cycle. 

Of  interest  to  note  also,  as  indicated  by  the  dependence  of  the 
magnitude  of  stress  on  At*  at  high  frequencies,  is  that  at  any  given 
value  of  Biot  number  the  stresses  decrease  with  increasing  plate 
thickness.  This  effect  also  is  opposite  to  the  general  effect  of 
dimensions  on  the  magnitude  of  thermal  stress,  as  indicated  for  instance 
by  the  dependence  of  stress  on  Biot  number.  It  is  postulated  that  at 
least  qualitatively,  high  thermal  inertia  is  required  for  low  magnitude 
of  thermal  stress.  Of  course,  high  thermal  inertia  is  achieved  by 
large  plate  thickness,  as  well  as  by  low  thermal  diffusivity. 
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It  is  appropriate  to  establish  the  applicability  of  the  present 
results  to  ceramic  practice,  by  means  of  a  numerical  example.  It  will 
be  assumed  that  the  value  for  the  Biot  number  is  (3  =  1.  At  this  value 
the  magnitude  of  At*  for  which  the  value  of  stress  at  cycle-to-cycle 
equilibrium  is  approximately  one  half  the  value  attained  for  At*  -►  ® 
corresponds  to  At*  ~  0.3,  with  corresponding  value  of  frequency 
f  =  (At)  =  </Q.3a  .  For  silicon  nitride  with  k  :  0.1  cm.  s  and  a 
value  of  a  =  0.25  cm,  this  frequency  is  approximately  5  cps  or  300  cpm, 
well  within  the  range  encountered  during  normal  operation  of  an  internal 
combustion  engine.  For  the  above  reasons,  it  is  predicted  that  structural 
ceramics  with  very  low  thermal  diffusivity  such  as  zirconia  or  micro- 
cracked  multi -oxide  ceramics  should  show  very  good  thermal  stress 
resistance  under  conditions  of  rapid  thermal  cycling. 

As  indicated  in  Fig.  6,  the  tensile  stresses  increase  with  increasing 
number  of  cycles  to  which  the  ceramic  is  subjected.  This  implies  that  if 
tensile  failure  is  of  primary  concern,  as  it  usually  is,  failure  more 
likely  will  occur  when  conditions  of  cycle-to-cycle  equilibrium  are 
reached  rather  than  during  the  first  few  cycles. 

This  conclusion  needs  revision  if  the  plate  originally  is  at  the 
higher  temperature  and  is  subjected  to  cooling  during  the  initial  phase 
of  the  first  cycle.  If  so,  the  signs  of  all  the  stress  values  reported  are 
reversed,  and  the  magnitude  of  the  tensile  stresses  during  the  first  cycle 
exceeds  that  for  cycle-to-cycle  thermal  equilibrium.  Under  this  condition, 
tensile  failure  is  more  likely  to  occur  during  the  first  few,  rather  than 
subsequent  cycles.  Obviously,  these  conclusions  would  need  to  be  re¬ 
examined  if  significant  thermal  fatigue  effects  are  present. 


For  purposes  of  material  selection,  it  is  convenient  to  derive  an 
appropriate  "thermal  stress  resistance  parameter"  on  the  basis  of  which, 
from  among  a  number  of  candidates,  the  material  can  be  selected  with 
optimum  thermal  stress  resistance  for  any  specific  failure  criterion 
and  mechanism  of  heat  transfer.  Because  of  the  specific  dependence  of 
o*  on  At*  found  in  the  present  study,  no  simple  linear  analytical 
expression  was  found  by  curve-fitting  or  other  method  from  which  the 
appropriate  thermal  stress  resistance  parameter  could  be  obtained. 
Nevertheless,  a  very  approximate  derivation  of  such  a  parameter  can  be 
presented.  This  derivation  will  concentrate  on  the  dependence  of  stress 
on  At*  at  the  condition  of  cycle-to-cycle  equilibrium. 

The  data  in  Fig.  9  illustrate  that  at  low  values  of  At*  (high 
frequencies)  the  magnitude  of  stress  increases  with  increasing  At*. 

I f  to  a  first  approximation  a  linear  dependence  can  be  assumed,  the 
stress  can  be  written: 

a  =  AaEATK(At)  (14) 

( 1-v)a^ 

The  data  in  Fig.  9  also  indicate  that  at  the  lower  values  of  the 
Biot  number  (8  <  1)  the  magnitude  of  maximum  stress  is  directly  proportional 
to  the  value  of  the  Biot  number,  8.  Incorporating  this  effect  in  eq.  14, 
yields: 

B«EATK(At)h  ,,  x 

o  *  TT-vVak (15) 

where  B  is  another  constant. 

Thermal  stress  failure  in  brittle  ceramics  usually  occurs  in  tension, 
Equating  the  stress,  o,  in  eq.  15,  to  the  tensile  strength,  S^,  an 
expression  can  be  obtained  for  the  maximum  range  of  temperature  difference 
ATmax  which  the  plate  can  be  subjected  for  a  given  value  of  AT: 


SA l-v)k 


B(At)h 


which  yields  the  thermal  stress  resistance  parameter: 


S+(1-v)k 


Because,  by  definition  <  =  k/pc  (where  p  is  the  density  and  c  is  the 


specific  heat),  this  parameter  can  be  written  also  as: 


S_j.(  I  -v)pc 


which  indicates  that  for  thermal  cycling  at  high  frequency,  materials 


should  have  a  high  specific  heat  per  unit  volume.  At  least  qualitatively. 


this  conclusion  seems  physically  realistic  since  a  high  specific  heat 


per  unit  volume  limits  the  range  of  temperature  variation  in  the  plate 


within  any  given  thermal  cycle.  Also,  high  specific  heat  per  unit  volume 


is  in  accordance  with  the  earlier  hypothesis  that  high  thermal  inertia 


is  required  to  minimize  the  stresses. 


In  summary,  the  present  results  indicate  under  conditions  of 


thermal  cycling  of  such  frequency  that  thermal  equilibrium  is  not 


achieved  within  each  cycle,  materials  with  low  values  of  thermal 


diffusivity  should  exhibit  better  thermal  stress  resistance  than  materials 


with  high  values  of  thermal  diffusivity. 


The  notation  Sf  for  tensile  strength  follows  the  original  notation  used 
by  Kingery  (footnote  1). 


-  A  _  *  tT m  —  •  .  —  - 
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ABSTRACT 


An  analysis  is  presented  of  the  thermal  stresses  in  brittle 
ceramics  subjected  to  thermal  shock  by  quenching  into  a  liquid 
medium  in  which  conduction  is  the  primary  mechanism  of  heat 
transfer.  The  magnitude  of  maximum  stress  was  shown  to  be  a 
function  of  the  coefficient  of  thermal  expansion,  Young's  modulus 
and  Poisson's  ratio  of  the  ceramic  as  well  as  the  values  of  the 
thermal  conductivity,  specific  heat  and  density  of  the  ceramic 
and  quenching  medium.  The  magnitude  of  thermal  stress  was  shown 
to  be  independent  of  specimen  size  and  on  cooling  to  occur  at  the 
instant  of  cooling  (t  =  0).  Appropriate  thermal  stress  resis¬ 
tance  parameters  were  derived. 


I.  INTRODUCTION 


Brittle  structural  ceramics  for  use  at  high  temperatures  are 
very  susceptible  to  catastrophic  failure  due  to  transient  or 
steady-state  thermal  stresses,  which  inevitably  arise  under  nor¬ 
mal  or  unusual  operating  conditions  involving  non-uniform  temper- 

1  2 

ature  contributions.  '  A  fundamental  understanding  of  the 
factors  which  affect  thermal  stress  fracture  of  brittle  ceramics 
is  imperative  for  reliable  performance  and  failure  analysis.  As 
an  alternative,  information  on  the  thermal  fracture  behavior  of 
ceramic  materials  can  be  obtained  by  experimental  means.  One 
test  for  this  purpose  which  is  highly  popular  because  of  its  sim¬ 
plicity,  consists  of  quenching  appropriate  specimens  from  an  oven 

3  4  5 

at  high  temperature  into  a  fluid  medium  such  as  water,  '  ' 

6  7  8  9 

oil,  fused  salts  or  fluidized  bed  '  maintained  at  a 

lower  temperature.  The  critical  quenching  temperature  difference 
(ATc)  required  to  initiate  thermal  stress  fracture  represents  a 
quantitative  measure  of  the  thermal  stress  resistance  of  the  par¬ 
ticular  ceramic. 

The  quantitative  interpretation  of  such  data,  without  excep¬ 
tion,  as  far  as  these  authors  are  aware,  is  based  on  thermo-elas¬ 
tic  theory  for  heat  transfer  by  convective  means.  The  assumption 
of  heat  transfer  by  convection  whether  by  natural  convection  for 
stationary  specimens  or  by  forced  convection  for  moving  speci¬ 
mens,10  implies  the  transport  of  heat  by  relative  motion 
(shear)  within  the  fluid.  It  should  be  noted  that  the  quantita¬ 
tive  understanding  of  convective  heat  transfer  is  based  on  meas- 


urement  during  steady-state  conditions.  In  this  respect  it 

should  be  noted  also  that  the  maximum  thermal  stresses  (i.e., 

fracture)  during  quenching  from  a  higher  to  a  lower  temperature 

occurs  a  relatively  brief  time  after  the  specimen  is  inserted 
3  7 

into  the  fluid  '  .  For  specimens  in  the  form  of  rods  with  cir¬ 
cular  or  square  cross-section  with  diameter  or  sides  of  a  few  mm 

and  values  of  thermal  diffusivity  ranging  from  0.02  to  0.2 
2  -1 

cm.  sec  ,  common  for  many  structural  ceramics,  the  time-to- 
maximum  stress  (i.e.,  time-to-f racture )  is  of  the  order  of  0.01 
seconds.  It  seems  unlikely  that  steady-state  conditions  for  con¬ 
vective  heat  transfer  could  be  established  in  a  time  span  of  this 

small  magnitude. 

In  order  to  provide  a  justification  for  the  objectives  of 
this  study,  it  should  also  be  recognized  that  a  thin  layer  of  the 
fluid  in  immediate  contact  with  the  specimen  is  stationary  with 
respect  to  the  hot  surface.  Any  heat  lost  from  the  specimen  must 
be  removed  by  conduction  through  this  stagnant  layer.  Such  con¬ 
duction  of  heat  requires  a  temperature  gradient  through  this 
layer  and  also  results  in  an  increase  in  temperature  of  this 
layer.  This  temperature  increase,  in  turn,  requires  an  amount  of 
heat  which  is  a  function  of  the  specific  heat  and  density  of  the 
fluid. 

As  indicated  by  the  values  of  peak  thermal  stress  and  asso- 

11  12  13 

dated  transient  temperature  distributions  '  '  on  cooling, 

tensile  thermal  stresses  in  the  surface  of  sufficient  magnitude 
to  cause  fracture  requires  the  removal  of  heat  from  a  relatively 


3 


thin  surface  layer  of  the  ceramic  specimen  which  is  a  small  frac¬ 


tion  of  the  total  thermal  energy  stored  within  the  specimen. 


Depending  on  the  relative  values  of  specific  heat  and  density  of 


the  specimen  and  fluid  medium  the  heat  removed  from  the  specimen 


surface  sufficient  to  cause  fracture  can  be  of  the  order  of  the 


amount  of  heat  absorbed  by  the  stagnant  layer.  If  this  hypothe¬ 


sis  is  correct  it  appears  conceivable  that  thermal  fracture  in  a 


quenching  experiment  may  not  require  the  transport  of  heat  by 


convection,  in  the  usual  sense.  Indeed,  if  this  statement  can  be 


proven  to  be  true,  quantitative  interpretation  of  experimental 


data  for  AT  should  be  based  on  conductive  heat  transfer  within 
c 


the  fluid,  rather  than  by  convective  heat  transfer  as  has  been 


common  practice. 


For  the  reasons  outlined  above  thermal  stress  fracture  in 


quenching  experiments  by  purely  conductive  processes  appears 


deserving  of  attention.  It  is  thought  that  during  a  quench. 


thermal  fracture  by  convection  and  conduction  are  competitive 


processes,  depending  on  the  relative  values  of  the  thermal  con¬ 


ductivity,  density  and  specific  heat  of  the  specimen  and  fluid 


medium  as  well  as  the  value  of  the  fluid  viscosity.  However,  in 


order  to  assess  the  dominant  mechanism  of  heat  transfer  responsi¬ 


ble  for  thermal  stress  fracture  in  a  quench,  an  analytical  solu¬ 


tion  for  the  thermal  stresses  generated  by  conductive  heat  trans¬ 


fer  is  required. 


The  purpose  of  this  study  was  to  derive  such  solutions  and 


to  assess  their  relevance  to  the  quenching  method  for  the  meas- 


urement  of  thermal  stress  resistance  of  brittle  structural  ceram¬ 
ics  . 

II.  ANALYSIS 

In  view  of  its  relevance  to  quenching  studies  a  circular 
cylinder  was  chosen  as  a  suitable  geometry  for  the  present  analy¬ 
sis.  The  cylinder,  at  initial  temperature  T  ,  was  assumed  to 
be  submerged  instantaneously  into  a  fluid  at  temperature  Tf. 
The  thermal  contact  between  the  cylinder  and  fluid  was  considered 
to  be  perfect  so  that  the  heat  transfer  was  not  impeded  by  a  con¬ 
tact  resistance.  The  transfer  of  heat  within  the  cylinder  and 
fluid  was  assumed  to  occur  in  a  radial  direction  only.  For  pur¬ 
poses  of  simplification  the  cylinder  was  assumed  to  be  infinitely 
long,  so  that  end  effects  could  be  ignored.  For  this  reason  the 
solutions  to  be  obtained  are  expected  to  describe  the  conditions 
near  the  center  of  quenching  specimens  with  a  length  large  com¬ 
pared  to  the  diameter.  The  solutions  for  the  transient  tempera¬ 
tures  were  obtained  first,  from  which  the  thermal  stresses  were 
then  readily  obtained. 


A.  Transient  Temperatures 

The  transient  temperatures  are  obtained  by  solution  of  the 
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heat  conduction  equations  : 


K  i  »_  (r  !!i,  .  i!i 

K1  r  9r  Sr  '  3t 

and 

K  i  i_  (r  !!ij  =  !!i 

k2  r  3r  'r  3r  '  3t 


(1) 


(2) 


where  8(r,t)  =  T  -  T^,  <  is  the  thermal  diffusivity,  r  is 

the  radial  coordinate  and  the  subscripts  1  and  2  refer  to  the 
cylinder  and  medium,  respectively. 

The  solution  of  eqs.  1  and  2  must  satisfy  the  boundary  con¬ 
ditions  : 


61  ( r, 0)  =  eQ  ;  02(r,O)  =  0 


T —  (0,t)  =  0  ;  Mm  6,(r,t)  =  0 


01(R,t)  =  82(R,t) 


oal 

ki  IT  <R't>  =k2^  <*'*> 


(3) 

(4) 

(5) 

(6) 


where  0O  =  T0  “  Tf,  R  is  the  radius  of  the  cylinder  and  k 
is  the  thermal  conductivity. 

A  general  solution  of  eqs.  1  and  2,  subject  to  the  boundary 
conditions  given  by  eqs.  3  thru  6,  was  presented  by  Carslaw  and 
Jaeger^  and  worked  out  in  detail  by  Hencke^.  The  final 
solution,  for  purposes  of  general  convenience,  was  expressed  in 
terms  of  the  thermal  conductivity  ratio,  K 


K  =  kl/k2  (7) 

the  thermal  diffusivity  ratio,  A 

a2  =  <^<2  (8) 

the  non-dimensional  time,  t* 

t*  =  k  jt/R2  <9> 

the  non-dimensional  radial  distance,  r* 

r*  =  r/R  (10) 

and  the  non-dimensional  temperature,  T* 


M  ta 


£<  « 


T*  =  (T-Tf)/(T  -Tf)  (M) 

For  convenience,  the  last  three  non-dimensional  quantities 
will  be  referred  to  as  time,  distance  and  temperature,  respec¬ 
tively. 

The  solutions  for  the  temperatures  T^*  and  T2*  in  the  cylin¬ 
der  and  surrounding  liquid,  respectively,  are: 


.  00  o  Jn  (ur*)  J,  (  m  )  dpi 

T,  (r* , t*)  =  —j  K  /  exp  (-  y  t*)  — = - j - 2 - 

1  7r  o  y^[$Z(y)  +  ^  (u)  ] 


(12a) 


T2  (r* , t* )  =  |  K  /  exp (- y2t*) 


J,  (y)  [J  (ayr*)  <J>  (y)  -Y  (ay  r*)  Q  (y)  ]  dy 


y  [ $2 (y )  +  n2(y)] 


(12b) 


where 


*(y)  -  K  J1(y)Yo(ay)  -  a  JQ(y)  Y^ay) 


fl(y)  =  K  J1(y)JQ(ay)  -  a  JQ(y)  ^  (ay) 


where  y  is  the  variable  of  integration  which  arises  from  invert- 

1  s 

ing  the  Laplace  transform  solutions  of  eqs.  1  and  2,  and 
J  ,  J^,  Yg,  and  Y^  are  the  Bessel  functions  of  the  first 
and  second  kind. 


a  J 

r. 


H  W 

d  a 
i 


The  temperatures  T^*  and  T2*  were  obtained  most  conven¬ 
iently  by  numerical  evaluation  of  the  integrals  in  eqs.  12a  and 

12b. 


B.  Thermal  Stresses 


General  expressions  for  the  thermal  stresses  in  long  circu¬ 
lar  cylinders  with  arbitrary  radial  temperature  distributions  are 
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given  by  Boley  and  Weiner  .  The  stress  distribution  consists 
of  radial,  tangential  and  axial  stresses.  From  the  point  of  view 
of  thermal  stress  failure  the  maximum  stresses  are  of  primary 
interest.  Examination  of  the  general  solutions  shows  that  the 

axial  stresses  exceed  the  radial  and  tangential  stresses  through¬ 
out  the  volume  of  the  cylinder  and  are  equal  to  the  tangential 
stresses  at  the  cylinder  surface  (r  =  R). 
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In  general, the  axial  stress  is  given  by  : 


°z  = 


-  TT^vT  I  f?  /  T*'ar'  - 

R  O 


where  a  is  the  coefficient  of  linear  thermal  expansion,  E  is 

Young's  modulus  and  v  is  Poisson's  ratio. 

For  convenience  the  results  will  be  reported  in  terms  of  the 

non-dimensional  stress  o  *,  defined  as: 

z 

_  oz(l-v) 

CTZ  “  aE (T  -T_)  (16) 

o  r 

In  this  notation,  substitution  of  eq .  (12a) into  eq.  (15)  yields: 


*  q  °°  o  J  -i  (y)dy 

az  (r* ' fc*)  *  h  K  /  **P<Vt*)  ~~ 3~~2~  - - ~57~7 

7i  o  w  [$  (u)  +  n  (y)  ] 


-  T.  (r*,t*) 


which  can  be  evaluated  numerically  as  eqs.  12. 
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III.  NUMERICAL  RESULTS  AND  DISCUSSION 


The  magnitude  of  the  transient  temperatures  and  stresses  are 
illustrated  most  conveniently  by  means  of  a  few  numerical  exam¬ 
ples.  For  two  of  these  examples  the  thermal  conductivity  of  the 
cylinder  will  be  taken  to  be  much  higher  or  much  lower  than  the 
thermal  conductivity  of  the  fluid.  For  the  third  example  the 
thermal  conductivity  of  the  cylinder  and  fluid  will  be  assumed  to 
be  equal.  The  stress  was  also  found  to  be  a  function  of  the 
volumetric  heat  capacity  ratio,  A  =  pici/p2c2  the  cyl“ 
inder  and  the  fluid.  However,  in  order  to  keep  the  number  of 
numerical  examples  to  a  reasonable  number,  the  volumetric  heat 
capacities  ratio  was  taken  as  unity. 

For  a  cylinder  being  cooled,  figs.  la,  lb,  and  lc  show  the 
transient  temperature  distributions  in  the  cylinder  and  fluid  for 
values  of  the  thermal  conductivity  ratio,  K  =  100,  1  and  0.01  , 

respectively,  for  a  range  of  values  of  time,  t* .  Comparison  of 
these  data  indicates  that  at  the  higher  value  of  K  the  tempera¬ 
tures  within  the  cylinder  at  any  instant  of  time  remain  fairly 
uniform.  This  effect  arises  because  the  heat  loss  from  the  cyl¬ 
inder  is  controlled  primarily  by  the  much  lower  value  of  the 
thermal  conductivity  of  the  fluid.  In  contrast,  at  the  value  of 
K  =  0.01,  because  of  the  much  higher  value  of  the  thermal  conduc¬ 
tivity  of  the  surrounding  liquid,  the  resulting  rapid  surface 
cooling  of  the  cylinder  results  in  a  relatively  high  degree  of 


non-uniformity  of  the  temperature  distribution  within  the  cylin¬ 
der.  For  K  =  1  the  degree  of  cooling  of  the  cylinder  and  heating 
of  the  fluid  at  the  surface  are  comparable. 


In  general,  the  magnitude  of  thermal  stresses  is  a  function 
of  the  degree  of  non-linearity  of  the  temperature  distribution. 

For  this  reason,  the  differences  in  the  temperature  distributions 
for  the  three  values  of  thermal  conductivity  ratio  shown  in  fig. 

L 

1  should  be  reflected  in  differences  in  the  magnitude  of  the 
thermal  stresses. 

S 

Figures  2a,  2b  and  2c  show  the  distributions  of  the  thermal 
stresses  within  the  cylinder  for  the  same  values  of  the  thermal  | 

conductivity  ratio  and  t*  as  those  for  the  temperature  distribu¬ 
tion  shown  in  figs,  la,  lb  and  lc.  As  expected  for  cooling,  the 

t 

maximum  tensile  stresses  are  found  at  the  surface  of  the  cylinder  ( 

(r  =  R),  whereas  the  maximum  compressive  stresses  are  found  at 
the  exact  center,  r  =  0.  For  heating,  the  signs  of  the  stresses 
are  reversed. 

The  maximum  value  of  thermal  stress  was  obtained  by  calcu¬ 
lating  the  stresses  as  a  function  of  time.  These  results  are 
shown  in  figs.  3a,  3b  and  3c,  with  the  stresses  in  the  surface 
being  tensile  or  compressive  for  cooling  and  heating,  respec¬ 
tively.  The  stresses  in  the  surface  reach  their  maximum  value  at 
t*  =  0.  In  contrast,  the  stresses  at  the  center  of  the  cylinder  1 

reach  their  maximum  value  after  some  time  has  elapsed. 

It  was  found  that  the  maximum  value  of  stress  in  the  surface 
could  be  expressed  by  a  simple  relation: 

o*(R)  =  (Nk1/2  +  l)"1 


where  N.  is  referred  to  as  the  conductive  heat  transfer  number 
k 

defined  by: 

=  klPlCl 

K  k2p2C2 

where  k,  p  and  c  are  the  thermal  conductivity,  density  and  spe¬ 
cific  heat,  respectively,  and  the  subscripts  1  and  2  refer  to  the 
cylinder  and  liquid  medium,  respectively.  Figure  4a  shows  the 
maximum  thermal  stress  in  the  surface  as  expressed  by  eq.  18,  as 
a  function  of  the  quantity  N^. 

The  "critical  quenching  temperature  difference,  ATc" 
required  for  the  initiation  of  thermal  stress  fracture  obtained 
in  a  quenching  experiment  is  a  measure  of  the  thermal  stress 
resistance  of  the  ceramic  being  tested.  At  the  value  of  AT  = 
ATc,  the  maximum  tensile  stress  equals  the  tensile  strength 
( St )  of  the  ceramic.  By  equating  this  value  of  thermal  stress 
to  the  tensile  strength,  as  is  common  practice  in  the  development 
of  thermal  stress  failure  theories  ,  with  the  aid  of  eq.  18  an 
expression  for  ATc  in  terms  of  the  relevant  material  properties 
can  be  derived: 


AT 

c 


St(l-v) 

rfWii 

1/2  -i 

i  1 

aE 

-'•k2P2C2- 

+  1 

(20) 


For  two  limiting  conditions,  eq.  20  permits  the  derivation 


of  thermal  stress  resistance  parameters  for  the  selection  of  the 
ceramic  v/ith  the  optimum  thermal  stress  resistance. 


the  appropriate  parameter 


For  k^P  ]_c]_/^2  P2C2 >>1 
in  terms  of  the  properties  of  the  cylinder  quenched  into  a  given 
fluid  is: 


St(l-v)(kpc)1/2 


aE 


(21) 


For  k^p  1c^/k2p2c2 '  the  appropriate  parameter 

i  s : 

St(l-v)/aE  (22) 

The  maximum  tensile  stress  at  the  center  of  the  cylinder 
encountered  during  heating  could  be  plotted  most  conveniently  as 
a  function  of  the  thermal  conductivity  ratio  (K)  and  the  volumet¬ 
ric  heat  capacity  ratio,  A,  as  shown  in  fig.  4b.  Because  of  the 
relatively  complex  dependence  of  the  stresses  on  K  and  A,  the 
derivation  of  thermal  stress  resistance  parameters  of  simple  form 
was  not  feasible. 

The  relevance  of  the  above  results  to  the  interpretation  of 
quenching  data  can  be  demonstrated  by  means  of  a  numerical  exam¬ 
ple  which  compares  the  magnitude  of  the  thermal  stresses  attained 
by  convective  and  conductive  heat  transfer.  For  this  example, 
four  different  liquid  media  were  selected  with  widely  different 
convective  and  conductive  heat  transfer  properties,  namely  sili¬ 
cone  oil,  water,  mercury  and  sodium.  Convective  heat  transfer 
can  occur  by  natural  (free)  or  forced  convection.  The  latter 


process  is  governed  by  the  relative  speed  between  the  cylinder 

and  the  fluid  medium.  Depending  on  the  choice  of  speed  almost 

any  value  for  the  forced  convective  heat  transfer  coefficient  can 

be  obtained.  The  value  of  speed  at  which  the  specimen  hits  the 

fluid  surface,  at  least  for  a  water  bath,  was  shown  in  a  recent 
17 

study  to  have  no  effect  on  ATc .  Perhaps  a  more  appropriate 
value  of  speed  is  the  terminal  velocity  at  which  the  specimen 
falls  through  the  fluid  medium  under  the  influence  of  grav¬ 
ity^.  Liquid  metal  media  such  as  lead  and  mercury,  however, 
have  a  density  higher  than  the  density  of  many  ceramic  specimens. 
This  difference  in  density  would  not  permit  quenching  by  the 
usual  free-fall  method,  as  the  specimen  would  float  only  par¬ 
tially  submerged  on  the  surface  of  the  fluid.  For  such  media 
with  high  density  it  appears  more  suitable  to  rapidly  immerse  the 
specimen,  attached  to  a  suitable  specimen  holder,  in  the  fluid 
after  which  it  is  held  stationary.  Under  this  condition  the 
specimen  is  subject  to  heat  flow  by  natural  convection,  which  was 
the  heat  transfer  conditions  selected  for  the  numerical  example 
to  be  presented  subsequently. 

The  heat  transfer  coefficient  under  conditions  of  free  con¬ 
vection  for  a  circular  cylinder  held  in  a  horizontal  stationary 

18 

position  was  calculated  from  equations  given  by  Holman 

For  the  silicone  oil  and  water  baths  the  heat  transfer  coef- 
1 8 

ficient  is  givey  by  : 

h  =  0.53(NGr*Npr)1/4(Kf/d) 
where  N„  is  the  Grashof  number: 


in  which  g  is  the  gravitational  constant,  3  is  the  coefficient 
of  volumetric  thermal  expansion,  d  is  the  diameter  of  the  cylinder 


and  v ^  is  the  viscosity  of  the  fluid. 

The  Prandtl  number,  Npr,  is  defined  by: 


NPr  =  VKf 


where  is  the  thermal  diffusivity  of  the  fluid. 
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For  the  liquid  metals,  the  heat  transfer  coefficient  is  : 
h  =  0.53(NGr  Npr2)1/4(kf/d) 

Because  the  values  of  h  given  by  eqs.  23  and  26  are  a  func¬ 
tion  of  the  temperature  difference  itself,  the  quantity 
(To-Tf)  was  taken  to  be  200°C  for  all  calculations.  The 
actual  value  of  (Tq-T^.)  can  be  obtained  experimentally  or  by 
the  method  of  successive  approximations.  For  the  water  bath  the 
heat  transfer  coefficient  was  assumed  to  result  from  free  convec¬ 
tion  only,  without  the  additional  contribution  from  nucleate 
boiling  or  steam-film  formation.  The  diameter  D  of  the  cylinder 
was  taken  to  be  6.35  and  12,7  mm  (  0.25  and  0.5  inch),  respec¬ 
tively,  which  represent  dimensions  typical  for  specimens  encoun¬ 
tered  in  laboratory  practice. 

The  thermal  stresses  for  conductive  heat  transfer  were  cal¬ 
culated  by  means  of  eqs.  18  and  19.  The  property  values  for  the 

water,  mercury  and  sodium  for  both  the  convective  and  conductive 

19 

heat  transfer  were  taken  from  values  tabulated  by  Chapman  for 


assumed  film  temperatures  of  140,  150  and  200°C,  respectively. 

-4  2- 

The  data  for  a  typical  silicone  oil  with  a  viscosity  of  10  ms 
were  taken  from  data  sheets  provided  by  a  manufacturer.* 

Table  I  lists  the  values  for  the  heat  transfer  parameters 
for  the  silicone  oil,  water,  mercury  and  sodium  required  for  the 
calculation  of  the  convective  heat  transfer  coefficient. 

Table  II  compares  the  values  of  maximum  thermal  stress  (a*) 
for  three  different  materials  with  property  values  typical  for  a 
glass,  polycrystalline  aluminum  oxide  and  silicon  carbide  at  a 
temperature  corresponding  approximately  to  the  initial  specimen 
temperature  during  an  actual  quench.  Comparison  of  these  stress 
values  indicates  that  for  all  quenching  media,  the  stresses  for 
conductive  heat  transfer  are  comparable  to  or  exceed  the  stress 
values  for  natural  convection  by  a  considerable  margin.  It 
should  be  emphasized  that  this  conclusion  is  valid  only  for  the 
specific  condition  of  natural  convection  chosen  for  the  mode  of 
convective  heat  transfer.  In  practice,  under  the  right  condi¬ 
tions  of  forced  convection,  almost  any  value  of  heat  transfer 
coefficient  and  corresponding  values  of  stress  (o*)  approaching 
unity  can  be  achieved. 

However,  if  the  quenching  test  is  conducted  by  holding  the 
specimen  stationary  in  a  non-stirred  fluid,  the  data  for  o*  in 
Table  II  indicate  that  the  possibility  of  thermal  fracture  by 
conductive  heat  tranfer  definitely  must  be  considered  in  the 
analysis  of  the  experimental  data.  Especially  for  liquid  metal 

*  Dow-Corning,  Midland,  Michigan 
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baths  with  their  much  higher  thermal  conductivity  used  for 

quenching  of  specimens  with  low  thermal  conductivity,  conductive 

heat  transfer  may  well  be  the  primary  mechanism  for  thermal 

stress  failure.  The  relative  contributions  of  convective  and 

conductive  heat  transfer  also  are  expected  to  be  affected  by 

specimen  size.  For  convective  heat  transfer  the  maximum  value  of 

thermal  stress  increases  with  specimen  size,  which  has  no  effect 

on  the  peak  stress  for  conductive  heat  transfer.  This  suggests 

that,  especially  for  small  specimens  commonly  used  for  laboratory 

studies,  the  contribution  of  conductive  heat  transfer  should  not 

be  overlooked.  If  indeed  conductive  heat  tranfer  is  the  primary 

cause  for  failure,  any  observed  effect  of  specimen  size  on  ATc 

is  attributable  to  the  statistical  nature  of  brittle  fracture, 

20 

such  as  expressed  by  the  Weibull  theory  .  Even  if  the  magni¬ 
tudes  of  peak  stress  due  to  convective  and  conductive  heat  trans¬ 
fer  are  comparable, failure  may  still  be  due  to  the  conductive 

heat  transfer  for  which  the  maximum  value  occurs  at  t  =  o; 

11  12  13 

whereas  for  convective  heat  transfer  '  '  ,  for  all  values  of 

the  Biot  number  6  <  00  ,  the  maximum  stress  occurs  at  t  >  o .  Ver¬ 

ification  of  the  analytical  results  of  this  study  is  handicapped 

by  the  absence  of  appropriate  data  in  the  open  literature.  The 

2 1 

only  available  data  for  specimens  held  stationary  were 
obtained  for  glass  thermally  cycled  in  water  as  the  quenching 
medium  for  which  the  interpretation  of  the  data  is  complicated  by 
extensive  sub-critical  crack  growth  and  quantitative  uncertain¬ 
ties  of  nucleate  boiling  and  steam-film  formation.  The  acquisi- 


tion  of  the  necessary  experimental  data  with  major  emphasis  on 
liquid  metals  is  intended  to  be  the  subject  of  a  future  study. 

As  a  final  comment  it  should  be  noted  that  eq.  20  for  T 
was  based  on  the  assumption  that  failure  would  occur  at  a  value 
of  maximum  tensile  stress  equal  to  a  well-defined  value  of  ten¬ 
sile  strength  without  regard  to  the  underlying  variables  which 
govern  or  contribute  to  tensile  failure.  It  should  be  noted  that 
making  this  assumption  has  been  common  practice  in  the  develop¬ 
ment  of  thermal  stress  failure  theories  of  brittle  cera- 
mics"*" '  ^ ^  >  25  ^  it  weil  known,  however,  that  tensile 

failure  of  brittle  materials  is  subject  to  statistical  and/or 
fracture-mechanical  variables.  For  this  reason  for  a  quantitative 
application  of  the  results  of  this  study,  a  separate  detailed 
analysis  of  the  tensile  failure  stress  will  be  required  for  the 
specific  ceramic  being  considered  and  for  the  specific  thermal 
conditions  to  which  it  is  being  subjected. 

In  summary,  an  analysis  was  conducted  of  the  variables  which 
affect  thermal  stress  failure  of  brittle  ceramics  subjected  to 
conductive  heat  transfer.  The  magnitude  of  the  peak  thermal 
stress  was  shown  to  be  a  function  of  the  coefficient  of  thermal 
expansion.  Young's  modulus  and  Poisson's  ratio  of  the  ceramic  as 
well  as  the  values  of  the  thermal  conductivity,  specific  heat  and 
density  of  the  ceramic  and  surrounding  medium.  Appropriate  ther¬ 
mal  stress  resistance  parameters  were  derived. 
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CRITERIA  FOR  THE  THERMAL  STRESS 
FAILURE  OF  BRITTLE  STRUCTURAL 
CERAMICS 


D.P.H.  Hasselman  and  J.P.  Singh 

Department  of  Materials  Engineering 
Virginia  Polytechnic  Institute  and 
State  University 
Blacksburg,  Virginia  24061  USA 

A  survey  is  presented  of  the  current  understanding  of  the 
response  of  brittle  structural  ceramic  materials  to  thermal 
stresses.  The  f igures-of-merit  available  at  this  time  for  the 
ceramic  with  optimum  thermal  stress  resistance  for  a  number 
of  thermal  boundary  conditions  and  failure  criteria  are 
reviewed.  These  f igures-of-merit  include  those  for  the 
initiation  of  thermal  stress  fracture  as  well  as  those 
appropriate  for  crack  propagation  and  arrest.  Experimental 
data  for  a  selection  of  brittle  ceramics  are  presented 
for  purposes  of  illustration  of  a  number  of  features  of  the 
nature  of  thermal  stress  fracture  of  these  materials. 


INTRODUCTION 

Structural  materials  for  numerous  applications  encountered 
in  industry,  aerospace  and  other  areas  of  technology  may  be 
required  to  provide  satisfactory  service  at  elevated  tempera¬ 
tures.  Such  applications  include  process  vessels  and  crucibles 
in  the  metal-working  and  chemical  process  industry,  heat  exchangers , 
heat-regenerators  or  recuperators,  components  for  energy-conversion 


devices  such  as  turbine  engines,  components  or  structure  for 
high-speed  aerospace  purposes,  and  many  others.  In  addition 
to  high  temperature,  the  materials  used  in  many  of  the  above 
applications  also  are  subjected  to  corrosion,  erosion  and 
mechanical  loads.  Obviously,  for  long-term  performance,  the 
materials  selected  for  these  applications  should  exhibit  high 
melting  point,  good  corrosion  and  erosion  resistance,  excellent 
retention  of  load-bearing  capability  as  well  as  resistance  to 
creep  at  high  temperature. 

The  only  class  of  materials  which  can  satisfy  the  above 
requirements  is  known  as  "structural  ceramics".  Such  ceramic 
materials  include  the  metal  oxides,  carbides,  nitrides,  silicides 
and  similar  compounds  as  well  as  carbon  and  graphite.  The  advan¬ 
tageous  properties  of  these  materials  at  high  temperature  basically 
derive  from  the  high  strength  of  the  bond  between  the  atoms. 

Such  bonding  usually  is  of  ionic  and/or  covalent  type.  Most 
unfortunately,  this  type  of  atomic  bonding  also  causes  these 
ceramic  materials  to  be  exceedingly  brittle.  This  brittleness 
introduces  a  great  deal  of  difficulty  in  the  manufacture  and  suc¬ 
cessful  use  of  structural  ceramics  in  those  applications  which 
require  close  design  tolerances  and  high  load-bearing  require¬ 
ments  . 

Structures  or  components  which  operate  at  elevated  temperatures 
invariably  will  be  subjected  to  transient  or  steady-state  heat 
flow  at  some  time  during  normal  use.  Such  heat  flow  usually  will 
involve  non-linear  temperature  distributions  which  inevitably  will 
give  rise  to  thermal  stresses.  The  combination  of  the  relevant 
properties  of  structural  ceramics  (i.e.,  low  thermal  conductivity. 


high  Young's  modulus)  is  such  that  these  thermal  stresses  are  of 
considerable  magnitude.  Because  of  the  high  degree  of  brittleness, 
these  thermal  stresses  are  not  relieved  by  non-linear  deformation. 
This  behavior  in  combination  with  a  low  fracture  stress  renders 
structural  ceramics  highly  susceptible  to  thermal  stress  failure. 
Furthermore,  the  extent  of  crack  propagation,  which  frequently 
occurs  in  a  highly  unstable  manner,  can  be  very  extensive.  In 
extreme  cases,  the  size  of  the  cracks  can  exceed  the  dimensions  of 
the  structure  or  component.  For  this  reason,  thermal  stress  fail  ire 
can  be  catastrophic,  rendering  the  material  completely  incapable  of 
continued  satisfactory  service.  Such  failure  is  demonstrated  in 
Fig.  1,  which  shows  a  series  of  thermally  shocked  spherical  test 
specimens  of  zirconium  carbide  with  varying  amounts  of  graphite 
dispersed  phase.  The  specimen  on  the  left  which  contains  no  graphite 
literally  exploded  on  fracture.  The  presence  of  the  graphite  phase 
inhibited  this  mode  of  failure,  although  the  initiation  of  fracture 
was  not  suppressed. 

Such  catastrophic  thermal  stress  failure  of  brittle  structural 
ceramics  constitutes  a  major  handicap  in  the  successful  design  of 
components  or  structures  intended  for  high-temperature  operation. 

It  is  ironic  that  the  very  class  of  structural  materials,  which  in 
view  of  a  combination  of  very  excellent  properties  is  most  eminently 
suited  for  high-temperature  use,  is  most  susceptible  to  catastrophic 
thermal  stress  failure. 

For  these  reasons,  the  development,  design  and  selection  of 
materials  for  high-temperature  application  requires  a  great  deal 
of  care.  The  role  of  the  pertinent  material  properties  and  oi her 


Fig.  1.  Catastrophic  thermal  stress  fracture  in  sphere 
of  zirconium  carbide  (left).  Other  spheres 
contain  a  dispersed  phase  of  graphite,  which 
inhibits  unstable  mode  of  crack  propagation 
(for  additional  details,  the  reader  is  referred 
to  reference  [l])  . 
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variables  which  affect  the  magnitude  of  thermal  stress  in  any  specific 
heat  flow  environment  must  be  well  understood,  both  qualitatively  as 
well  as  quantitatively.  All  possible  modes  of  failure  must  be  con¬ 
sidered.  Furthermore,  a  detailed  understanding  of  the  factors  which 
control  crack  propagation  following  the  initiation  of  fracture  is 
essential  for  predicting  service  performance  following  thermal  stress 
failure. 

The  purpose  of  this  paper  is  to  present  the  state-of-the-art  under¬ 
standing  of  the  thermal  stress  failure  of  brittle  structural  cerai lies. 
The  point  of  view  of  the  materials  engineer  or  scientist  responsible 
for  the  ceramics  with  optimum  thermal  stress  resistance  is  emphasized. 
The  underlying  thermo-elastic  theory  is  covered  by  other  studies  of 
this  series.  Much  of  the  experimental  work  in  the  area  of  thermal 
stress  fracture  of  brittle  structural  materials  was  done  under  poorl 
defined  heat  transfer  conditions.  For  this  reason,  no  exhaustive 
review  of  these  results  will  be  presented.  Instead,  for  the  purpose 
of  clarity  and  convenience  of  the  more  theoretically  oriented  reader 
of  this  series  only  those  experimental  data  which  will  illustrate  a 
particular  phenomenon  of  thermal  failure  will  be  presented. 
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CRITERIA  FOR  OPTIMUM  THERMAL  STRESS  RESISTANCE 
General 

At  this  time  the  selection  of  materials  with  optimum  thermal  stress 
resistance  is  based  on  two  distinctly  different  (as  will  be  shown) 
criteria . 

The  first  criterion  is  based  on  the  objective  that  the  magnitude 
of  thermal  stress  relative  to  the  failure  stress  be  reduced  at  much 
as  possible.  In  other  words,  the  basic  premise  of  this  criterion 
is  to  avoid  the  initiation  of  fracture.  On  the  basis  of  this 


6 

criterion,  appropriate  selection  rules  for  materials  with  optimum 
thermal  stress  resistance  can  be  developed.  For  this  purpose, 
simple  expressions  for  the  maximum  heat  transfer  conditions  based 
on  the  most  simplifying  assumptions  were  derived  to  indicate  the 
relative  role  of  the  pertinent  material  properties  in  governing 
thermal  stress  resistance.  From  these  expressions  appropriate 
f igures-of-merit  expressed  in  terms  of  these  material  properties 
then  form  the  basis  for  comparison  for  the  relative  thermal  stress 
resistance  of  candidate  materials  for  a  specific  application. 

Critical  to  note  is  that  the  derivation  of  such  f igures-of-merit 
should  not  be  regarded  as  design  calculations  which  is  so  only  for 
the  most  simple  case. 

In  view  of  the  high  susceptibility  of  brittle  structural  ceramics 
to  catastrophic  thermal  stress  failure,  the  derivation  of  such  f igures- 
of-merit  for  optimum  thermal  stress  resistance  appears  to  have  re¬ 
ceived  greater  attention  for  purpose  of  high-temperature  technology 
than  in  other  fields  of  science  or  engineering.  In  this  respect,  then, 
one  of  the  major  objectives  of  this  paper  is  to  review  the  thermal 
stress  f igures-of-merit  derived  to  date. 

The  above  approach  in  which  the  thermal  stresses  are  kept  as  small 
as  possible  relative  to  the  failure  stress  (or  alternatively:  the 
stress  intensity  factor  is  kept  below  the  critical  stress  intensity 
factor)  is  common  practice  in  the  area  of  structural  engineering. 

Most,  if  not  all,  readers  of  this  paper  should  have  no  reason  to  take 
issue  with  this  criterion.  In  practice,  however,  materials  can  be 
subjected  to  thermal  stresses  which  can  be  an  order-of-magnitude 
or  more  above  the  stress  required  for  fracture.  Clearly,  under  such 
conditions  the  initiation  of  failure  cannot  be  avoided.  Order-of- 
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magnitude  improvements  in  the  failure  stress  of  brittle  materials 
at  the  current  state  of  technology  is  not  expected  within  the  near 
or  even  distant  future.  At  first  sight,  then,  it  appears  that  in 
terms  of  conventional  engineering  practice  no  solutions  for  conditions 
of  excessive  thermal  stresses  appear  to  be  available. 

Nevertheless,  a  solution,  believed  to  be  unique  to  high  tempera¬ 
ture  technology,  is  based  on  the  criterion  that  since  the  initiation 
of  fracture  cannot  be  avoided,  the  resulting  crack  propagation 
should  be  limited  as  much  as  possible.  Basically  this  criterion 
relies  on  the  arrest  of  cracks  after  the  onset  of  propagation 
and  is  commonly  used  in  so-called  "refractory  practice".  Refractories 
are  the  type  of  materials  used  as  linings  of  blast-furnaces,  kilns, 
thermal  storage  units,  and  other  large  structures  which  must  operate 
at  elevated  temperature  for  long  duration  and  for  many  of  which 
heat  insulation  is  a  critical  aspect  of  their  performance. 

The  latter  criterion  also  can  be  used  to  derive  appropriate 
figures-of-merit .  These,  of  course,  because  of  their  basically 
different  foundation  are  expected  to  differ  from  those  derived 
from  the  criterion  that  the  initiation  of  thermal  fracture  must 
be  avoided.  In  this  respect,  then,  the  second  objective  of  this 
paper  is  to  present  the  state-of-the-art  understanding  of  the 
material  properties  and  other  variables  which  control  crack  arrest 
during  thermal  stress  failure.  In  the  view  of  these  writers,  this 
latter  problem  represents  the  most  promising  area  for  further 


theoretical  work. 


DERIVATION  OF  FIGURES— OF-MERIT 


General  Approach 

The  major  objective  of  the  derivation  of  f igures-of-meri t  is 
to  establish  the  role  of  the  pertinent  material  properties  which 
affect  thermal  stress  failure  and  to  develop  appropriate  rules 
for  the  selection  of  candidate  materials  with  optimum  thermal 
stress  resistance  for  a  particular  thermal  environment,  mode  of 
heat  transfer  and  failure.  Since  such  derivations  are  not  made 
for  design  purposes,  they  can  be  based  on  the  most  simplifying 
assumptions,  geometries  and  boundary  conditions.  All  material 
properties  will  be  assumed  to  be  independent  of  temperature  and 
position  within  the  material.  With  one  exception,  linear  elastic 
behavior  till  fracture  will  be  assumed  throughout.  Since  the 
geometry  does  not  affect  the  relative  role  of  the  material  prop¬ 
erties  for  any  given  thermal  environment,  the  appropriate  figure- 
of-merit  will  be  given  for  a  single  geometry  only,  in  order  to 
avoid  redundancy.  Poisson's  ratio  generally  plays  a  lesser  role 
in  establishing  the  magnitude  of  thermal  stress  than  other  mater¬ 
ial  properties  as  it  usually  occurs  in  factors  of  the  form  (l-v) 
or  (l-2v).  For  this  reason,  again,  to  avoid  redundancy,  a  single 
expression  involving  Poisson's  ratio  will  be  given  for  any  given 
thermal  environment  and  failure  criterion. 

Brittle  structural  materials  generally  exhibit  a  ratio  of 
tensi le-to-compressive  failure  stress  equal  to  an  order  of  magni¬ 
tude  or  so.  For  this  reason,  thermal  stress  failure  far  more 
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likely  will  occur  under  the  influence  of  a  tensile  thermal  stress 
than  by  the  compressive  thermal  stresses.  Such  tensile  failure 
will  be  assumed  throughout  Wj  th  the  exception  of  the  case  of 
external  constraints,  which  result  in  a  compressive  thermal 
stress  only.  For  the  f igures-of-meri t  for  limited  crack  propaga¬ 
tion,  the  derivation  will  be  based  on  either  mode  I  tensile  or 
mode  II  shear  crack-opening. 

The  general  approach  taken  in  the  development  of  a  specific 
f igure-of-merit  for  a  given  heat  transfer  condition  and  mode  of 
failure,  is  to  derive  or  obtain  from  the  literature  the  appropri¬ 
ate  expression  for  the  thermal  stresses  in  terms  of  the  relevant 
parameters  such  as  geometry,  size  and  boundary  conditions.  These 
latter  can  include  the  rate  of  heating,  magnitude  of  heat  flux 
and  many  others.  The  maximum  value  of  thermal  stress  should  not 
exceed  the  failure  stress  (such  as  the  tensile  strength).  By 
substitution  of  this  latter  value  for  the  maximum  stress,  fol¬ 
lowed  by  rearrangement,  an  expression  is  obtained  for  the  maximum 
thermal  conditions  ( i . e . ,  maximum  rate  of  heating)  to  which  the 
material  can  be  subjected  without  fracture.  Implicit  in  such  an 
approach  is  the  assumption  that  the  safety  factor  equals  unity. 
For  higher  values  of  the  safety  factor,  the  appropriate  expres¬ 
sion  can  be  obtained  by  a  simple  division. 

In  the  final  expressions  which  have  appealed  in  the  litera¬ 
ture,  the  relevant  material  properties  and  other  parameters  such 
as  geometry  and  dimensions  usually  are  uncoupled;  this  permits 
the  easy  extraction  of  the  appropriate  f igure-of-merit .  For  the 
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convenience  of  the  reader  and  to  keep  this  report  within  reason¬ 
able  size,  the  final  expression  for  the  maximum  thermal  condition 
will  be  given  together  with  the  appropriate  literature  reference 
if  available.  For  those  cases  which  involve  more  complex  rnateri- 
als-related  phenomena  such  as  sub-critical  crack  growth  and 
others,  with  which  some  of  the  readers  may  not  be  fully  familiar, 
a  brief  explanation  will  be  given  as  well. 

The  above  assumptions  and  approach  permit  expressing  the 
f igures-of-merit  in  simple  analytical  form.  For  those  cases 
where  the  stress  depends  non-linearly  on  some  relevant  parameter 
(as  in  the  case  for  the  dependence  of  the  thermal  stress  on 
optical  thickness  of  partially  transparent  materials  subjected  to 
radiation  heating)  limiting  conditions  are  selected  such  that  the 
f igure-of-merit  can  be  expressed  in  simple  algebraic  form.  For 
the  more  general  case,  the  reader  is  referred  to  the  original 
study.  For  the  thermal  conditions  and  failure  criteria,  figures- 
of-merit  derived  and  reported  in  the  literature  are  listed 
below. 

Figures-of-Merit  for  the  Resistance  to  the  Initiation  of 

Thermal  Stress  Fracture 

A.  Uniform  temperature  or  steady-state  heat  flow 

Case  1.  Flat  plate  at  uniform  temperature  constrained  from  in-plane 
expansion. 

The  flat  plate  initially  stress-free  is  at  uniform  tempera¬ 
ture  and  is  prevented  from  thermal  expansion  by  rigid  constraints 
at  the  edges  of  the  plate.  The  plate  will  be  considered  to  be 
sufficiently  thick  so  that  failure  due  to  elastic  instability 


(i.e.,  thermal  buckling)  is  avoided.  The  plate  undergoes  a  uni¬ 
form  temperature  change  AT.  For  a  positive  coefficient  of  ther¬ 
mal  expansion  (as  is  usually  the  case  for  most  materials)  the 
stresses  in  the  plate  will  be  tensile  or  compressive  depending  on 
whether  AT  is  negative  or  positive.  The  maximum  value  of  AT  to 
which  the  plate  can  be  subjected  is: 

AT  =  a.  (l-v)/aE  cooling  (la) 

in  3.x  z 

AT  =  a_(l-v)/aE  heating  (1L) 

m3  x  c 

iT^ax  °f  eqs.  la  and  lb  is  independent  of  the  size  of  the  plate 
and  has  a  geometric  constant  equal  to  unity. 

The  appropriate  f igures-of-meri t  are: 

ofc(l-v)/aE  and  0c(l-v)/aE  (lc) 

For  a  uniaxial  or  triaxial  constraint  the  term  (1-v)  in  eg. 
lc  is  replaced  by  unity  and  (l-2v),  respectively. 

Case  2.  Composite  structures. 

High  levels  of  thermal  stress  even  under  isothermal  condi¬ 
tions  can  be  generated  in  composites  and  structures  composed  of 
strongly  bonded  dissimilar  materials  with  mismatches  in  their 
coefficients  of  thermal  expansion.  Such  a  composite  or  structure 
will  be  stress-free  only  at  the  temperature  at  which  the  materi¬ 
als  which  constitute  the  composite  or  structure  were  fabricated 
or  joined.  A  change  from  this  value  of  temperature  creates  ther¬ 
mal  stress,  as  each  individual  material  is  prevented  from  free 
expansion  by  the  other  material(s).  Such  stresses  usual. ,  are 


undesirable,  with  the  possible  exception  of  the  bi-metallic  strip 
for  which  the  resulting  deformation  can  be  used  to  advantage  for 
the  operation  of  a  thermostat. 


A  simple  expression  can  be  given  for  the  stresses  which 
result  from  a  mismatch  in  the  coefficients  of  thermal  expansion 
for  a  composite  structure  consisting  of  a  thin  coating  on  an 
underlying  thick  substrate.  The  magnitude  of  the  stresses  in  the 
coating  equals  [2]: 

a  =  EAaAT/(l-v)  (2a) 


where  Aa  is  the  difference  in  the  coefficients  of  thermal  expansion 
of  the  substrate  and  the  coating  and  AT  is  the  difference  in 
temperature  at  which  the  coating  was  joined  to  the  substrate 
and  the  new  temperature  to  which  the  composite  is  heated 
or  cooled.  For  other  composite  structures  such  as  inclusions  [3,4] 
and  concentric  rings  or  cylinders,  the  expressions  for  the  stresses 
which  result  from  a  mismatch  in  the  coefficients  of  thermal  expansion 
are  more  complex  and  involve  the  elastic  properties  of  all 
individual  components  or  materials.  Regardless  of  the  geometry, 
however,  the  magnitude  of  these  stresses  is  proportional  to 
the  quantities  Aa  and  AT  and  can  be  negative  or  positive 
depending  on  the  sign  of  Aa  and  AT. 

For  failure  in  tension,  eq.  2a  can  be  rearranged  to  yield 
the  maximum  temperature  difference  to  which  the  coating  and  sub- 


Case  3 . 


Concentric  hollow  cylinder.  Radially  outward  steady-state 
heat  flow. 


The  concentric  hollow  cylinder  infinite  in  extent  with  inner 
and  outer  radius  of  a  and  b,  respectively,  is  subjected  to  radi¬ 
ally  outward  heat  flow.  This  results  in  a  radial  distribution  of 
temperature  with  a  temperature  difference  AT  between  the  inner 
and  outer  surfaces.  The  cylinder  is  not  constrained  externally. 
From  the  solutions  of  Timoshenko  and  Goodier  [5],  the  maximum 
value  of  temperature  difference  (ATmax)  across  the  cylinder  wall 


which  can  be  permitted  is: 


AT 


at(l-v) 


21n 


!)][- 


2a‘ 


b2-a~ 


In 


,-i-l 


(3a) 


max  aE 

By  noting  that  the  heat  flux,  q  per  unit  length  of  the  cyl¬ 
inder  can  be  related  to  the  temperature  difference  AT  across  the 
wall  by: 


q  =  2irKAT/ln  (b/a) 


(3b) 


where  K  is  the  thermal  conductivity,  eq.  3a  for  AT  can  be 
rearranged  to  yield  an  expression  for  the  maximum  heat 
flux  per  unit  length,  <3max/  to  which  the  cylinder  can  be  subjected 
given  by: 


4tto^(1-v)K 


2a' 


.2  a  2 

b  -a 


In 


n  -1 


(3c) 


Mmax  aE 

For  a  cylinder  of  finite  length,  the  tensile  stresses  exhi¬ 
bit  their  maximum  value  at  the  end  of  the  cylinder  in  the  tangen¬ 
tial  direction.  The  appropriate  values  of  AT  and  q  can  be 

max  ^  max 

obtained  by  dividing  eqs.  3a  and  3c  by  [6]: 


The  f igures-of-merit  for  AT  and  of  eqs.  3a  and  3c 

Ilia.  X  ITlci  X 


are : 


at(l-v)/aE  and  a^.(l-v)K/aE  ( 3e] 

Case  4.  Uniform  internal  heat  generation. 

Thermal  stresses  of  high  magnitude  can  be  caused  by  the 
non-uniform  temperature  distributions  resulting  from  internal 
heat  generation.  This  occurs,  for  instance,  in  such  components 
as  nuclear  fuel  elements,  and  in  materials  subjected  to  micro- 
wave  heating,  nuclear  or  solar  radiation. 

For  steady-state,  spatially  uniform  internal  heat  genera¬ 
tion,  the  role  of  material  properties  can  be  obtained  from  the 
solutions  for  a  solid  circular  cylinder  [71.  The  maximum  rate  of 

internal  heat  generation  per  unit  volume  (H  )  which  can  be 

max 

withstood  is: 

H  =  8a . ( 1-v ) K/aEb2  (4a) 

max  t 

with  the  corresponding  f igure-of-merit : 

a^d-vjK/aE  (4b) 

Eq.  4a  indicates  that  in  order  to  reduce  the  possibility  of 
tensile  thermal  fracture,  the  cylinder  radius  (b)  should  be  kept 
as  small  as  possible. 

Case  5.  Thermal  discontinuities. 

Steady-state  heat  flow  in  materials  with  a  linear  tempera¬ 


ture  distribution  and  absence  of  external  constraints  will  result 


in  zero  thermal  stress.  If,  however,  such  materials  contain 
inclusions  with  heat  conduction  properties  different  from  the 
matrix  materials,  the  temperature  field  near  these  inclusions 
will  be  disturbed  to  be  non-linear.  This,  as  indicated  by  Flo¬ 
rence  and  Goodier  [8]  for  a  spherical  cavity  and  by  Tauchert  [9] 
for  a  general  inclusion,  will  result  in  a  thermal  stress  field  in 
the  immediate  vicinity  around  the  inclusion.  From  the  solution 
for  the  spherical  cavity,  the  maximum  temperature  gradient 

which  can  be  imposed  without  incurring  the  risk  of  tensile  fai¬ 
lure,  can  be  derived: 

V  =  2a.  (l-v)/aEb  (5a) 

max  t  ' 

where  V  is  the  temperature  gradient  in  the  absence 

of  the  cavity  and  b  is  the  cavity  radius. 

The  heat  flux  (q)  can  be  related  to  the  temperature  gradient 

by:  q  =  K  7.  Accordingly,  eq.  5a  can  be  rewritten  in  terms  of 

the  maximum  allowable  heat  flux,  q  : 

'  ^  max 

q  =  2ct .  K ( 1-v ) /aEb  (5b) 

^max  t 

Eqs .  5a  and  5b  yield  the  f igures-of-merit : 

o^(l-v)/aE  and  a^.K(l-v)/aE  (5c) 

Case  6.  Thermal  buckling.  Straight  column  with  uniform 
temperature . 

From  the  solutions  of  Fridman  [10],  and  Burgemeister  and 
Stuep  [11],  the  critical  temperature  difference  (  ATc  )  required 
for  thermoelastic  instability  can  be  expressed  as: 


where  I  is  the  cross-sectional  moment  of  inertia,  L  is  the 
length  of  the  column,  a  is  the  coefficient  of  thermal  expansion 
and  A  is  the  cross-sectional  area. 

The  value  of  the  constant  depends  on  the  boundary  conditions 
and  is  equal  to  unity  for  a  column  for  which  both  ends  are  free 
to  rotate.  For  a  column  where  the  ends  are  not  permitted  to 
rotate  the  value  of  C-^  is  4.0  and  for  a  column  with  one  end  free 
and  the  other  end  restrained  from  rotation  its  value  is  0.25. 

Eq.  6a  reveals  that  the  critical  temperature  difference  is 
an  inverse  function  of  the  coefficient  of  thermal  expansion  only. 

A  thermal  buckling  resistance  f igure-of-merit  can  be  defined  as: 


Case  7.  Post-thermal  buckling.  Straight  column  at  uniform 
temperature . 


When  the  temperature  rise  in  a  straight  column  with  ends 
free  to  rotate  exceeds  the  critical  temperature  AT_.,  it  does  not 
fail  immediately,  but  exhibits  post-thermal  buckling  analyzed  by 
Boley  and  Weiner  [12].  It  will  fail  only  when  the  maximum  bend¬ 
ing  stresses  in  the  column  exceed  the  tensile  strength.  Under 
this  condition,  the  maximum  temperature  difference,  AT  over 
v/hich  a  column  v/ith  square  cross-section  can  be  heated,  is  given 


by  [  13 


2  2  2  2  2 

AT  v  =  AT  +  o.  IT/tt  aE"a 
max  c  t 


Eq.  7a  suggests  the  post-thermal  buckling  f igure-of-merit : 


2,  c-2 

(?£  /aE 
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(7b) 
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B.  Transient  heat  flow 

Case  8.  Convective  Newtonian  heat  transfer.  Instantaneous 
change  in  ambient  temperature. 

An  infinitely  long  solid  circular  cylinder  initially  at  thermal 
equilibrium  is  subjected  to  an  instantaneous  change  AT  in  ambient 
temperature  under  conditions  of  Newtonian  convective  heat  transfer. 
Under  these  conditions  the  maximum  tensile  stress  occurs  at  the 
surface  of  the  cylinder. 

From  the  solutions  of  Jaeger  [14 J  and  Manson  and  Smith  [15]  the 
maximum  change  ATmax  in  ambient  temperature  for  a  range  in  values 
of  the  Biot  number  0  <  8  <  10  to  a  good  approximation  can  be 
expressed : 

1.45a.  (1-v) 

ATmax  =  - st - {1  +  3'41/B)  (8a) 

where  8  =  Rh/K  with  R  being  the  radius,  h  the  heat  transfer  coefficient 

and  K  the  thermal  conductivity. 

Eq.  8a  indicates  that  the  relative  effect  of  the  Biot  number 

on  AT  is  a  function  of  its  magnitude, 
max  ^ 

For  3.41/8  <<  1  eq.  8a  becomes: 

AT  v  s  1.45a.  <l-v)/aE  (8b) 

Ilia  X  L 

Similarly  for  3.41/8  <<  1,  eq.  8a  becomes 

AT  „  “  4.95a.  (l-v)K/ctEbh  (8c) 

Ilia  X  L 

Eqs .  8b  and  8c  suggest  the  f igures-of-meri t : 

o^.  (1-v)  K/aE 


at(l-v) /aE 


( 8d) 


The  relative  influence  of  the  above  figures-of-merit  on  AT 

max 

depends  on  the  magnitude  of  the  Biot  number. 

Case  9.  Controlled  rate  of  heating. 

An  infinite  flat  plate  is  asymmetrically  cooled  by  a 
constant  rate  (T)  of  surface  temperature.  From  the  solutions 
of  Buessum  [16 J  the  maximum  rate  of  cooling  so  that  tensile 
fracture  can  be  avoided  is: 

Tmax  =  3ot(l-v) K/aEb2  (9a) 

where  b  is  the  half-thickness  of  the  plate,  and  k  is  the  thermal 
dif fusivity . 

Eq.  9a  yields  the  f igure-of-merit : 

a^(l-v) K/aE  (9b) 

Radiative  Heat  Transfer 
Case  10.  Opaque  materials. 

A  solid  sphere  initially  at  thermal  equilibrium  at  low 
temperature  is  suddenly  subjected  to  black-body  radiation  at 
temperature  T.  If  thermal  stress  fracture  occurs  for  surface 
temperatures  of  the  sphere  low  enough  that  emitted  radiation  can 
be  neglected,  the  solutions  of  Hasselman  1.17]  show  that  the 
maximum  radiative  heat  flux  qmax  the  sphere  can  be  subjected  to  is: 

qmav  =  5 a.  (l-v)K/uErb  (10a) 

IUclX  t- 

where  b  is  the  radius  of  the  sphere  and  c  =  1-r  where  r  is  the 
reflectivity. 


Ji 


Eq.  10a  can  be  rewritten  in  terms  of  the  maximum  radiation 
temperature,  Tmax  to  which  the  sphere  can  be  subjected. 


4  4 

q  =  P  (T  -  T  ) 

amax  max 


(10b) 


.  »\  v\ 

£  f 

S..N 


where  T  is  the  surface  temperature  of  the  sphere  and  p  is  the 

Stefan-Boltzmann  constant. 

For  T  >  >  T : 

max 


q  =  pT 
nmax  max 


Substitution  of  eq.  10c  for  q  in  eq.  10a  followed  by 

ITlclX 


( 10c) 


rearrangement  yields: 


1/4  1/4 

T  =  (5/pb)  '  {a.  (l-v)K/cxEe} 

IHaX  U 


(lOd) 


Equations  10a  and  lOd  yield  the  f igures-of-merit : 


p t ( 1-v ) K/aE  e 


{ (1-v) K/aE e  } 


1/4 


(lOe) 


tn 


Case  11.  Semi-transparent  materials. 

The  material  is  assumed  to  be  completely  transparent  below  a 

value  of  wavelength  (XQ)  and  totally  opaque  above  this  wavelength  [  18 ] 

Thermal  stresses  will  arise  only  for  radiative  heat  transfer  with 

wavelength  X  >  X  .  Solution  for  q  and  T  can  be  obtained 
^  o  ^max  max 

from  the  earlier  example  of  a  totally  opaque  sphere  initially  at 
low  temperature  suddenly  subjected  to  black-body  radiation: 

q  =  5a  (1-v) K/aEeb (1-F.  )  (11a) 

Mmax  t  Xq 

and 

T  =  (5a.  (l-v)K}1//4/faEcbp  (1-F.  )}1/4  (lib) 

ma  x  t  a 
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where  is  the  fraction  of  energy  in  the  black-body 

o 

spectrum  for  A  £  A  . 

Eqs.  11a  and  11b  lead  to  the  f igures-of-meri t : 

ofc  (1-v)  K/ctEe  (1-  Fa  )  ;  f  afc  (1-v)  K/aEe  ( 1-F  )  l1^4  (11c) 

Ao  o 

Semi-absorbing  materials. 

Radiative  heat  transfer  is  considered  which  penetrates  into 
a  material  with  an  intensity  (q)  as  a  function  of  depth  of  pene¬ 
tration  (x)  given  by: 

q  =  qQe  exp(-yx)  (12a) 

where  qQ  is  the  intensity  of  the  incoming  radiation,  e  =  1-r 
where  r  is  the  reflectivity  and  y  is  the  absorption  coefficient. 
The  resulting  spatially  non-uniform  internal  heat  generation 
gives  rise  to  thermal  stress.  Solutions  for  the  magnitude 
of  these  stresses  were  obtained  in  a  number  of  studies  (19-21) 
specifically  for  an  infinite  flat  plate  originally  at  thermal 
equilibrium  subjected  simultaneously  to  normally  incident 
thermal  radiation  and  cooled  by  Newtonian  convection. 

Temperatures  were  considered  to  remain  low  enough  that 
re-emitted  radiation  over  the  total  thermal  history  was 
negligible  compared  to  the  intensity  of  the  incoming  radiation. 


Appropriate  f igures-of-merit  of  simple  analytical  form  were 
derived  for  the  limiting  conditions  of  the  heat  transfer  coeffi¬ 
cient  h  =  0,  ™  and  optical  thickness,  pa  ->  0,  (where  a  is  the 
plate  half-thickness)  under  conditions  at  which  the  plate  has 
reached  thermal  eqiuilbrium  after  long-term  radiation. 

Case  12.  Asymmetric  heating  and  cooling.  h=  0;  pa  -*  ®. 

The  maximum  heat  flux  (q  _)  to  which  the  plate  can  be  sub- 

^max 

jected  in  order  to  avoid  tensile  failure: 

0 

q  =  12  0  (l-v)K/aEea  (12b) 

max  t 

Case  13.  Asymmetric  heating  and  cooling.  h  =  0;  pa  ->  0 

2 

qmax  =  at(1“y)K/aEeya  (13a) 
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C.  Time-to- failure 

•i 


In  practice,  situations  can  be  encountered  in  which  the  time 
period  required  for  satisfactory  performance  is  of  the  order  of 
the  period  of  time  required  for  the  thermal  stresses  under  tran¬ 
sient  heat  transfer  conditions  to  rise  to  their  maximum  value. 

This  is  the  case,  for  instance,  for  components  of  aerospace 
structures,  with  a  time  of  performance  on  the  order  of  seconds. 

If  thermal  stress  fracture  does  not  occur  over  the  time  period 
over  which  this  component  is  expected  to  provide  satisfactory 
service,  it  is  considered  to  be  safe  even  at  longer  values  of 
time. 

Case  17.  Sphere  subjected  to  an  instantaneous  rise  in  temperature. 

The  time  period  (t*)  required  for  the  tensile  thermal 
stresses  at  the  center  of  the  sphere  to  reach  their  maximum  value 
as  calculated  by  Gruenberg  [22]  is  given  by: 

t*  =  0.0574  b2/<  (17a) 

where  b  is  the  sphere  radius.  Eq.  17a  indicates  that  the 
time  of  maximum  stress  is  an  inverse  function  of  the  thermal 
diffusivity.  This  latter  statement  is  generally  true  for  thermal 
stresses  encountered  in  convective  heat  transfer  involving  an 
instantaneous  change  in  ambient  temperature.  For  these  condi¬ 
tions,  the  magnitude  of  maximum  stress  is  uniquely  defined  by  the 
value  of  the  Biot  number  [12,14].  However,  the  actual  time 
period  required  for  the  stresses  to  reach  a  prescribed  value  is 
an  inverse  function  of  the  thermal  diffusivity.  Accordingly, 
the  f igure-of-merit  appropriate  for  having  large  times-to-failure  is 


3 


N. 
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defined  by: 


■c  (17b) 

D.  Thermoviscoelastic  Stress  relaxation 

Under  conditions  of  mechanical  stress,  materials  can  exhibit 
slow  deformation  by  creep.  Such  creep  deformation  also  is 
expected  to  occur  under  conditions  of  thermal  stress.  Ceramic 
materials  will  exhibit  creep  by  diffusional  processes  at  levels 
of  temperature  at  which  vacancy  concentrations  and  mobility  become 
appreciable.  These  temperatures  correspond  to  about  one  half  to 
two  thirds  of  the  melting  point  of  the  material. 

Creep  deformation  under  conditions  of  thermal  stress  gener¬ 
ally  is  favorable  as  it  leads  to  stress  relaxation.  For  this 
reason,  the  occurrence  of  creep  reduces  the  probability  of  fai¬ 
lure  by  thermal  stresses.  The  magnitude  of  thermal  stresses  due 
to  external  constraints  will  decrease  with  increasing  time. 
Therefore,  thermal  stress  relaxation  by  creep  leads  to  an 
increase  in  the  thermal  condition  (i.e.,  AT  ,  q  ,  etc.)  to 
which  a  material  can  be  subjected  without  thermal  stress  failure, 
to  a  value  greater  than  if  no  stress  relaxation  occurred. 


In 

general 

,  the  higher  the 

creep 

rate , 

the 

more  rapidly 

thermal 

stress 

relaxation  will 

occur 

,  with 

a 

corresponding 

greater 

thermal 

stress  resistance. 

The 

derivation 

of  analytical 

expressions  for  the  magnitude  of  thermal  stress  affected  by  creep 
is  highly  complex,  if  possible  at  all.  In  general,  the  kinetics 
of  creep  is  a  function  of  the  temperature  as  well  as  stress.  For 
this  reason,  in  a  body  subjected  to  non-isothermal  conditions. 


thermal  stress  relaxation  can  occur  by  a  multitude  of  mechanisms. 
Complex  numerical  methods  based  on  finite  element  principles  are 
required  for  such  calculations  (Poritsky  and  Fend  L23]) .  For 
these  reasons,  derivations  of  expressions  for  the  purposes  of 
obtaining  f igures-of-merit  for  thermal  stress  relaxation  in  gen¬ 
eral  cannot  be  obtained  unless  grossly  simplifying  assumptions 
are  introduced.  Even  then,  such  an  analysis  should  be  limited  to 
steady-state  heat  flow. 

In  such  an  analytical  study  [24],  it  was  noted  that  the  rate 
of  creep  in  a  thermally  stressed  structure  corresponding  to  the 
stress  at  the  position  of  lowest  temperature  represents  a  lower 
bound  on  the  rate  of  thermal  stress  relaxation.  This  concept  was 
applied  to  a  hollow  circular  cylinder  undergoing  steady-state 
heat  flow,  treated  earlier.  It  was  pointed  out  that  as  thermal 
stress  relaxation  occurred,  the  values  of  the  maximum  temperature 
difference  or  heat  flux  across  the  wall  could  be  increased,  at  a 
rate  such  that  the  maximum  value  of  tensile  thermal  stress 
remained  constant. 

Case  18.  Hollow  cylinder.  Radial  outward  heat  flow- 

The  maximum  rate  of  rise  of  the  temperature  difference 

(AT)  across  the  cylinder  wall  for  a  ceramic  material  undergoing 
max 

linear  creep  was  derived  to  be  [24 J: 


(AT) 


max 


(1-v) /Can 


(18a) 


where  n  is  the  effective  viscosity  which  relates  the  creep  rate 
(e)  to  the  stress  by  i  =  o/n  and 


C  =  [2*n(|) ]-1[l 

a 


2a 
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b 
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-a 
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In  (|)] 


(18b) 


Similarly,  the  heat  flux  per  unit  length  can  be  increased  at 


a  rate: 


where 


d  =  o . ( 1-v ) K/C ' an 
Mmax  t 


C'  =  [1  -  ~ — 2“  ln  i^/47T 
b  -a 


(18c) 

(18d) 


Eqs .  18a  and  18c  suggest  that  materials  which  exhibit  maxi¬ 

mum  thermal  stress  relaxation  should  be  selected  on  the  basis  of 
the  f igures-of-merit : 

cr.  (1-v)  a,  (l-v)K 

— -  and  — -  (18e) 

an  an 

Important  to  note  is  that  equations  18a  and  18c  represent 
lower  bounds  on  the  rate  of  rise  of  A  ljnax  or  <lmax-  As  time  pro¬ 
gresses  and  the  changes  in  AT  and  q  involve  changes  in  the  abso¬ 
lute  temperature,  these  rates  can  be  modified.  Also,  critical  to 
note  is  that  thermal  stress  relaxation  under  non-isothermal  con¬ 
ditions  can  lead  to  the  generation  of  residual  stresses  of  magni¬ 
tude  comparable  to  the  original  thermal  stress  whenever  condi¬ 
tions  of  thermal  equilibrium  are  achieved  on  cessation  of  the 
flow  of  heat. 

E.  Thermal  Fatigue 

Brittle  ceramic  materials  under  conditions  of  static  or 
cyclic  stress  can  exhibit  the  growth  of  cracks  at  levels  of 
stress  intensity  factor  below  the  critical  stress  intensity  factor 


(i.e.,  for  Kj  <  KIc) .  This  phenomenon  is  referred  to  as  sub- 
critical  crack  growth.  If  for  a  given  stress  value  or  history, 
the  total  extent  of  crack  growth  is  sufficient  for  the  crack 
to  become  critical,  failure  will  ensue.  For  constant  or  cyclic 
stress,  this  mode  of  failure  is  referred  to  as  static  or  cyclic 
fatigue,  respectively. 

Subcritical  crack  growth  can  also  occur  under  conditions  of 
steady-state  or  cyclic  thermal  stress.  The  resulting  failure  is 
referred  to  as  "thermal  fatigue".  Such  fatigue  behavior  can  be 
predicted  from  fracture-mechanical  principles.  Basically,  this 
methodology  requires  a  knowledge  of  the  size  and  geometry  of  the 
propagating  crack,  the  kinetics  of  crack  growth  as  well  as  the 
stress  and  temperature  history  to  which  the  material  is  being 
subjected.  In  principle,  with  this  information,  the  time  or  num¬ 
ber  of  cycles  required  for  the  propagating  crack  to  become  critical 
(i.e.,  the  time  to  failure)  can  be  calculated  for  any  stress  and 
temperature  history. 

Case  19.  Steady-state  heat  flow. 

An  expression  for  the  time  to  failure  under  conditions  of  a 
constant  thermal  stress  was  derived  [25 J  as  follows: 

The  time-to-failure  (tf)  for  a  brittle  material  subjected  to 
a  constant  tensile  stress  (o)  is  given  by  [ 2 6 J . 

K 

tf  =  [ 2/a2Y2  /  Ic [KjdKj/v]  (19a) 

Kli 

where  Y  is  a  geometric  constant  which  relates  the  value  of  tensile 
stress  to  the  mode  I  stress  intensity  factor,  Kj  by: 


where  a  is  the  crack  size.  =  K^(t  =  0)  and  V  is  the  crack 

velocity  which  for  many  brittle  materials  can  be  expressed 
as : 

V  =  AKn  exp  (-Q/RT)  (19c) 

I 

where  A  and  n  are  constants  and  Q  is  the  activation  energy  for 
the  particular  mechanism  responsible  for  the  crack  growth. 

An  expression  for  the  thermal  stresses  can  be  assumed  to  take  the  form 

a  =  CoiEAT/  (1— v )  (19d) 

where  C  is  a  geometric  constant.  For  instance,  for  a  hollow 
circular  cylinder  subjected  to  steady-state  radial  heat  flow, 
the  quantity  C  is  given  by  eq.  3a  presented  in  an  earlier  section 
of  this  report. 

Substitution  of  eq.  19d  into  eq.  19a  and  assuming  that 
KIi  <<  KIc,  where  =  K^.  at  t=0  results  in: 

tf  =  2 ( 1-v ) n  exp  (Q/RT)/[CaEATY]nA(n-2) aL {n-2>/2  (19e) 

where  a^  «  a  at  t  =  0 . 

Eq.  19a  can  be  rewritten  in  terms  of  the  time-to-f ailure  at 
a  given  level  of  heat  flux  (q)  by  using  the  relation  q  =  C-^KAT 
where  is  another  constant  and  K  is  the  thermal  conductivity. 

This  yields: 

tf  =  2(l-v)nC^Kn  exp  (Q/RT)/CaEqY)nA(n-2)  (n_2) /2 


( 19  f ) 


The  appropriate  thermal  fatigue  f igures-of-meri t  can  be 
obtained  from  egs.  19e  and  19f  by  the  deletion  of  the  constant  C 
and  ,  AT  and  q,  respectively. 

One  should  be  careful  to  note  that  the  above  approach  is 
based  on  the  implicit  assumption  that  the  thermal  stress  inten¬ 
sity  factor  can  be  calculated  from  the  value  of  thermal  stress 
obtained  from  thermoelastic  principles  for  a  continuum  and  the 
crack  size  and  geometry  as  given  by  eq.  19b.  This  assumption  is 
expected  to  lead  to  reliable  results  as  long  as  the  crack  size  is 
small  (say  S  1%)  relative  to  the  size  of  the  component  or  struc¬ 
ture  which  is  the  case  for  high-strength  structural  ceramics. 
For  crack  sizes,  however,  which  are  an  appreciable  fraction  of 
the  body  size,  the  above  approach  will  lead  to  an  over-estimate 
of  the  thermal  stress  intensity  factor,  due  to  the  effect  of  the 
crack  on  the  effective  compliance.  As  will  be  shown  later,  this 
effect  is  very  critical  for  the  promotion  of  crack  arrest.  Such 
dependence  can  be  established  by  analytical  or  numerical  techni¬ 
ques  as  those  employed  by  Stahn  et  al  [27,28],  Stern  [29]  and 
Emery  and  co-workers  [30,31].  Since  the  present  approach 
neglects  the  effect  of  the  crack  on  the  compliance,  the  expres¬ 
sions  for  the  time-to-failure,  in  fact,  are  conservative. 

Case  20.  Cyclic  temperature  variation. 

The  prediction  of  thermal  fatigue  behavior  under  conditions 
of  cyclic  temperature  variation  is  far  more  complex  than  for 
steady-state  heat  flow.  Under  cyclic  conditions.  both  the 
stresses  as  well  as  the  temperature  vary  continuously  as  governed 
by  the  mode  of  heat  transfer,  geometry  and  other  variables. 


The  magnitude  of  absolute  temperatures  encountered  is  very 
critical  in  view  of  the  thermally  activated  nature  of  sub-criti¬ 
cal  crack  growth.  For  these  reasons,  the  prediction  of  thermal 
fatigue  under  cyclic  conditions  is  most  conveniently  performed  by 
the  numerical  integration  of  eq.  19a.  This  method  was  followed 
successfully  in  a  number  of  earlier  studies,  as  will  be  indicated 
later  in  the  section  of  experimental  results.  Such  numerical 
integration,  however,  does  not  give  a  direct  feel  to  the  materi¬ 
als  engineer  for  the  relative  role  of  the  material  properties  and 
could  handicap  an  intelligent  choice  of  the  candidate  material 
with  optimum  thermal  fatigue  resistance  for  any  given  applica¬ 
tion. 

A  semi-empirical  equation  based  on  eq.  19a  was  proposed  [32] 
for  the  number  of  cycles  (N)  required  for  failure  of  a  solid  cir¬ 
cular  cylinder  under  conditions  of  cyclic  change  in  ambient  temp¬ 
erature  involving  Newtonian  convective  heating  and  cooling. 
Sub-critical  crack  growth  was  assumed  prominent  during  the  cool¬ 
ing  cycle.  The  cylinder  was  assumed  to  have  attained  thermal 
equilibrium  between  cycles.  The  number  of  cycles  to  failure  (N) 
is : 


N  = 


Ck (exp  0/RTmax) 


A a2  R2Y2(n-2)K  . (n_2) 
max  Ii 


(20a 


where  C  is  a  constant  to  be  evaluated  by  numerical  integration 
of  eq.  19a.  Tmaxis  the  maximum  temperature  encountered  by  the 
material,  which  in  this  case  is  the  higher  of  the  two  temperatures 
between  which  the  cylinder  is  being  cycled.  °max  is  the  maximum 
value  of  tensile  thermal  stress  encountered  in  the  cylinder 


M 


during  the  cooling  part  of  the  cycle.  K  .  is  the  maximum 

Ii 

value  of  mode  I  stress  intensity  factor  during  the  first  cycle: 

K  .  =  a  Ya.1//2  (2 

Ii  max  1 

where  a^  is  the  initial  crack  size. 

Substitution  of  eq.  20b  into  eq.  20a  yields: 

Ck  exp (Q/RT  ) 

M  -  _ inaX _ _ _ _  /  ^ 

Y"(n-2)ailn-;!)/2 

with  omax  for  a  solid  circular  cylinder  (see  case  8)  given 


°max  =  «EAT/1.45(l-v) {1+3.41/0} 


Figures-of-Merit  for  the  Resistance  to 
Crack  Propagation  (K  >  K  ) 

General 


As  mentioned  earlier  in  the  introduction,  thermal  environ¬ 
ments  can  be  encountered  in  practice  which  are  of  such  severity 
that  even  in  the  best  materials,  the  onset  of  thermal  stress 
fracture  cannot  be  avoided.  For  such  conditions,  good  thermal 
stress  resistance  of  brittle  structural  materials  is  based  on 
their  ability  to  retain  their  geometry  and  physical  characteris¬ 
tics  so  that  continued  satisfactory  performance  can  be  rendered. 


in  spite  of  being  fractured.  Experimental  tests  of  thermal 
stress  resistance  based  on  this  criterion  may  measure  the  relative 
loss  in  load-bearing  ability  such  as  the  loss  of  tensile  strength 
due  to  a  single  thermal  cycle,  the  number  of  cycles  required  to 
cause  complete  fracture  or  the  loss  in  weight  due  to  a  single 
thermal  shock  due  to  surface  spalling  and  similar  criteria. 


Clearly,  thermal  shock  resistance  based  on  these  latter  criteria 
is  governed  by  the  extent  of  crack  propagation  following  the  ini¬ 
tiation  of  thermal  stress  fracture.  The  derivation  of  figures- 
of-merit  appropriate  to  such  criteria  requires  the  understanding 
of  the  pertinent  material  properties  which  affect  the  propagation 
of  cracks  in  thermal  stress  fields.  The  literature  studies  con¬ 
cerned  with  this  latter  problem  are  fewer  in  number  than  those 
which  deal  with  the  analysis  of  thermal  stresses  and  the  onset  of 
crack  propagation. 

Case  21.  Flat  plate  with  Griffith  crack.  Steady-state  heat 
flow  perpendicular  to  plane  of  crack. 

From  the  solutions  of  Florence  and  Goodier  [8],  the  maximum 

temperature  gradient  (V  )  permitted  to  avoid  crack  instability 

ms  x 

can  be  derived  to  be: 

Vmax  =  8^/2*a2E£3)  1/2  (21a 

where  y  is  the  total  energy  required  to  create  unit  area  of 
fracture  surface  and  £  is  the  half-length  of  the  crack. 

Eq.  21a  leads  to  the  f igure-of-meri t : 

(Y/a2E)1/2  (21b 

Eq.  21a  can  be  rewritten  in  terms  of  the  maximum  heat  flow 

(Q  )  per  unit  width  of  the  plate: 
max 

0mav  *  8(YK2/27ia2E£3)  1/2  (21c 


with  the  corresponding  figure-of-merit : 

(YK2/a2E)  1/2  (21i 

Case  22:  Three-dimensional  solid  with  penny-shaped  crack. 

Steady-state  heat  flow  perpendicular  to  crack. 

From  the  solutions  of  Nowinski  [33],  the  maximum  temperature 

gradient  (  V_  )  to  avoid  crack  instability  in  terms  of  fracture- 
max' 

mechanical  principles  can  be  written: 


V 

max 


(22i 


where  is  the  mode  II  critical  stress  intensity  factor 

and  b  is  the  radius  of  the  crack. 

The  figure-of-merit  is: 


Clearly,  from  a  general  point  of  view  this  case  is  identical 
to  the  previous  case  21  because  the  critical  stress  intensity 
factor  Kc  can  be  related  to  the  fracture  energy  y  and  Young's 
modulus,  E  by: 

Kc  =  (2yE)  1/2  ( 22< 


This  latter  expression  permits  conversion  of  equations  for 
the  crack  stability  in  thermal  stress  field  expressed  in  terms  of 
energy  principles  to  fracture-mechanical  concepts  and  vice-versa. 


♦For  crack-interaction  in  plane  of  crack  propagation,  the  reader 
is  referred  to  the  study  of  Singh  et  al.  [35]. 


For  the  purpose  of  further  discussion,  it  should  be  noted 
that  if  the  initial  crack  length,  lQ  >  £m  the  crack  will  propa¬ 
gate  in  a  stable  manner  with  the  crack-length  and  aT  related  by 
means  of  eq.  23b.  If,  however,  lQ  <  £m  the  crack  will  propagate 
in  an  unstable  (dynamic)  mode  as  described  below. 

Case  24:  Flat  plate  with  N  parallel  non-interacting  cracks  of 

equal  size  per  unit  area.  Plate  is  cooled  by  temperature 
difference  AT  and  prevented  from  shrinking  by  uniaxial 
constraints  perpendicular  to  plane  of  cracks. 

Unstable  crack  propagation. 

This  case  is  identical  to  Case  23  considered  above,  with  the 

exception  that  the  initial  crack  size  l  <<  £  (eq.  23e). 
c  o  m 

On  crack  instability  at  AT  given  by  eq.  23a  the  final 

1  max  1  ^ 

crack  length  (£f  )  which  results  from  unstable  crack  propagation 


£f  =  (4ttN£q) 


(24a) 


Eq.  24a  indicates  that  is  inversely  proportional  to  the 
initial  crack  length,  £q  and  the  crack  density,  N  but  independent 

from  all  other  material  properties.  This  latter  conclusion 
requires  modification  when  the  fracture  energy  or  stress  inten¬ 
sity  factors  for  the  onset  of  crack  propagation  and  arrest  are 
not  the  same . 

Eq.  24a,  by  means  of  the  Griffith  equation  can  be  rewritten 
in  terms  of  the  original  tensile  strength,  a  ,  which  yields: 


£  f  =  a  t  /8NyE 


(24b) 
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The  appropriate  f igure-of-meri t  is: 

afc2/NyE  (24c 

Again,  by  using  the  Griffith  equation,  eq  (24b)  can  be 
rewritten  in  terms  of  the  tensile  strength  retained  (o^)  follow¬ 
ing  unstable  crack  propagation  as  a  function  of  the  original 
strength  (c  )  : 

a-  =  (yE/o  )  (16N/ir)1/2  (24d 

This  latter  result  indicates  that  the  tensile  load-bearing 
ability  retained  following  dynamic  crack  propagation  is  inversely 
proportional  to  the  corresponding  load-bearing  ability  prior  to 
crack  propagation.  This  conclusion  has  considerable  practical 
implications . 

Figure  2  illustrates  the  crack  stability  and  propagation 
behavior  described  by  eqs.  23  and  24  for  two  values  of  crack  den¬ 
sity.  Although  a  very  simple  mechanical  model  was  used  for 
these  derivations,  the  general  behavior  shown  in  Fig.  2  has  led 
to  the  prediction  of  a  number  of  fracture  phenomena  in  the  ther¬ 
mal  stress  failure  of  brittle  materials.  These  will  be  presented 
in  a  subsequent  section  devoted  to  experimental  behavior. 

The  absence  of  crack-interaction  was  assumed  throughout  in 
the  above  discussion.  Nevertheless  such  interactions  can  be 
taken  into  account  as  shown  by  Singh  and  co-workers  (35J  for  an 
array  of  co-planar  cracks  based  on  the  solutions  of  Yokobori  and 
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Crack  stability  and  propagation  behavior  in 
uniaxially  constrained  flat  plate  with  parallel 
non-interacting  cracks,  cooled  through  a 
temperature  difference  AT  (after  Hasselman  [34]) 
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Ichikawa  {36].  Figure  3  shows  the  dependence  of  4T  on  the  ratio 
of  crack  size  (a)  to  crack  spacing  (d)  as  calculated  by  Singh  et 
al.  for  a  range  of  values  of  the  ratio  of  row-to-row  spacing  (b) 
to  crack  size  (d).  At  low  values  of  a/d  and  high  values  of  b/d, 

AT  depends  on  crack  length  in  a  manner  analogous  to  dependence 
shown  in  Fig.  2.  As  a/d  ->-0.5,  AT  decreases  to  zero  due  to 
crack  coalescence. 

GENERAL  COMMENTS  ON  FIGURES -OF-MERIT 

Presented  above  were  a  total  of  twenty-four  cases  for  the 
thermal  stress  failure  of  brittle  ceramics  for  non- redundant 
modes  of  heat  transfer  and  failure  criteria.  Undoubtedly,  in  the 
near  or  far  future  other  cases  will  be  recognized  and  additional 
f igures-of-merit  will  be  derived.  Regardless  of  the  total  num¬ 
ber,  it  is  clear  that  the  solution  of  the  problem  of  the  high 
susceptibility  of  brittle  structural  ceramics  to  thermal  stress 
failure  is  complex  indeed.  For  any  specific  application,  the 
choice  of  the  material  with  optimum  thermal  stress  resistance 
must  be  preceded  by  a  thorough  analysis  of  the  mechanism(s)  of 
heat  transfer,  the  resulting  magnitudes  and  distribution  of  the 
thermal  stresses  as  well  as  the  most  probable  mode  of  failure, 
all  in  light  of  the  basic  performance  requirements  of  the  compo¬ 
nent  or  structure  in  question.  Because  such  an  analysis  involves 
the  disciplines  of  heat  transfer  and  mechanics  as  well  as  materials 
engineering,  the  problem  of  thermal  stress  fracture  is  most  effi¬ 
ciently  tackled  by  a  multi-disciplinary  approach. 


CRITICAL  TEMPERATURE  DIFFERENCE 


Fortunately,  from  the  point  of  view  of  the  materials  tech¬ 
nologist,  the  material  properties  which  affect  thermal  stress 
fracture,  at  least  for  the  twenty-four  cases  presented  above,  are 
clearly  identified  and  are  uncoupled  from  those  variables  which 
affect  the  transfer  of  heat  to  the  component  or  structure.  Even 
then,  approximately  a  dozen  and  a  half  material  properties  or 
parameters  enter  the  multitude  of  f igures-of-meri t .  A  complete 
characterization  of  any  specific  material  for  the  purpose  of 
predicting  its  thermal  fracture  behavior  is  a  task  of  considera¬ 
ble  magnitude.  Such  a  task  becomes  even  larger  when  it  is  recog¬ 
nized  that  most,  if  not  all,  mechanical,  thermal  or  optical  prop¬ 
erties  or  parameters  which  make  up  the  f igures-of-merit ,  are  a 
function  of  temperature,  wavelength,  environmental  conditions, 
microstructure,  composition  or  impurity  level,  loading  conditions 
and  many  other  variables.  Nominal  property  values  obtained  from 
engineering  handbooks  at  best  are  useful  for  first  approximations 
only . 

In  this  respect,  the  failure  stress  of  brittle  materials 
probably  represents  the  property  subject  to  the  largest  degree  of 
uncertainty.  A  detailed  review  of  the  failure  phenomena  of  brit¬ 
tle  materials  is  beyond  the  scope  of  this  study.  In  general, 
extensive  experimental  data  have  shown  that  the  failure  stress  of 
such  materials  can  be  a  function  of  the  stress  state  (i.e.,  uni¬ 
axial  vs.  multi-axial),  the  stress  distribution  as  well  as  the 
stressed  volume  (i.e.,  statistical  aspects),  loading  rate,  temp¬ 
erature,  past  mechanical  or  thermal  history,  environmental 


(stress-corrosion)  effects  and  numerous  other  variables.  Adding 
to  these  effects  is  the  fact  that  failure  of  brittle  materials 
originates  from  structural  defects  introduced  during  processing 
or  subsequent  mechanical  handling.  Such  defects  usually  are  too 
small  to  be  adequately  characterized  by  non-destructive  techni¬ 
ques  currently  available  for  the  detection  of  flaw  in  metallic 
components.  Because  of  the  uncertainty  in  the  magnitude  of  the 
failure  stress  coupled  with  the  quantitative  uncertainties  of 
many  modes  of  heat  transfer  such  as  those  encountered  in  tran¬ 
sient  convection  and  radiation,  calculation  of  the  anticipated 
thermal  fracture  behavior  of  brittle  structural  ceramics,  with 
the  exception  of  perhaps  the  most  simple  case,  must  be  regarded 
as  order- of- magnitude  estimates  only. 

Nevertheless,  the  equations  and  associated  f igures-of-merit 
presented  above  provide  valuable  guidelines  for  the  design  of 
components  and  selection  of  appropriate  materials.  Minimizing 
the  component  size  or  thickness  in  many  instances  reduces  the 
magnitude  of  thermal  stress,  to  result  in  a  corresponding 
increase  in  the  severity  of  the  heat  flux  (or  other  criteria)  to 
which  the  component  can  be  subjected.  Failure  due  to  thermo¬ 
elastic  instability,  which  can  be  reduced  by  improving  component 
thickness,  probably  represents  the  only  exception  to  the  above 
rule. 

The  optimum  combination  of  material  properties  which  yields 
high  resistance  to  the  initiation  of  thermal  stress  failure 
depends  on  the  role  of  these  properties  in  maximizing  the  figure- 


of-merit  for  a  specific  thermal  environment,  performance  criter¬ 
ion  (AT  ,  q  ,  etc.)  and  failure  mode.  Low  values  of  the  coef- 
ficient  of  thermal  expansion.  Young's  modulus  and  Poisson's  ratio 
are  always  preferred,  in  combination  with  as  high  a  failure  stress 
as  possible.  Whenever  appropriate,  the  value  of  thermal  con¬ 
ductivity  and  thermal  diffusivitv  should  be  as  high  as  possible. 
Regardless  of  which  material  property  is  considered,  no  incom¬ 
patibility  -xists  in  its  role  for  the  f igures-of-merit  presented 
above  for  the  resistance  to  the  initiation  of  thermal  fracture. 

Such  an  incompatibility,  however,  can  be  found  to  exist  in  the 
f igures-of-merit  related  to  crack  stability  and  propagation.  High 
crack  stability,  in  general,  requires  a  low  value  of  the  coefficient 
of  thermal  expansion  and  Young's  modulus,  combined  with  high  values 
of  the  fracture  energy.  These  latter  requirements  of  course  are 
not  unexpected  as  all  f igures-of-merit  involving  the  tensile  fracture 
stress  can  be  converted  into  their  corresponding  fracture-mechanical 
analogue,  by  means  of  the  relation  between  fracture  strength,  Tung's 
modulus  and  fracture  energy,  appropriate  for  a  given  crack  geometry, 
if  known. 

In  contrast,  examination  of  the  f igure-of-merit  appropriate 
to  unstable  crack  propagation  (eq.  24c)  shows  that  rapid  crack 
arrest  (i.e.,  small  values  of  final  crack  length,  £f)  requires  a 
low  value  of  original  tensile  fracture  stress  and  a  high  value  of 


Young's  modulus  of  elasticity.  This  latter  result  derived  in 
terms  of  fracture-mechanical  principles  is  intuitively  understan¬ 
dable,  if  it  is  recognized  that  the  elastic  energy  at  fracture  is 
proportional  to  the  ratio  of  the  square  of  the  fracture  stress  to 
Young's  modulus.  Low  and  high  values  of  these  two  latter  quanti¬ 
ties  will  lead  to  low  values  of  elastic  energy  at  fracture  and 
corresponding  rapid  crack  arrest. 

The  above  incompatibility  in  the  role  of  the  fracture  stress 
and  Young's  modulus  in  the  resistance  to  the  initiation  of  ther¬ 
mal  fracture  on  one  hand  and  rapid  crack  arrest  following  unsta¬ 
ble  crack  propagation  on  the  other  hand,  presents  a  dilemma  to 
the  materials  technologist  in  the  development  and  selection  of 
the  optimum  ceramic. 

Materials  selection  based  on  maximizing  the  resistance  to 
the  initiation  of  fracture  by  increasing  the  fracture  stress  will 
also  lead  to  an  increasing  extent  of  crack  propagation  for  ther¬ 
mal  conditions  of  such  severity  that  failure  cannot  be  avoided. 
If  this  latter  occasion  arises  as  is  frequently  the  case  in  prac¬ 
tice,  material  selection  should  be  based  on  reducing  the  extent 
of  crack  propagation,  which  can  be  achieved  by  selection  of 
materials  with  moderate  to  low  strength.  This  will  result  in  an 
improvement  of  thermal  shock  resistance  in  terms  of  retained  phy¬ 
sical  characteristics  and  geometry  following  thermal  stress  frac¬ 
ture.  However,  this  approach  will  reduce  the  total  load-bearing 
ability,  prior  to  or  following  the  initiation  of  thermal  stress 
fracture,  required  to  fulfill  the  basic  function  of  the  component 


for  any  specific  application.  For  this  reason,  material  selec¬ 
tion  based  on  crack  arrest  appears  most  appropriate  for  materials 
subjected  to  minimal  mechanical  load  during  service.  This,  in 
fact,  is  the  case  for  refractory  liners  for  high- temperature  pro¬ 
cess  vessels.  For  these  types  of  materials  without  exception  the 
criterion  for  high  thermal  shock  resistance  is  based  on  rapid 
crack  arrest  as  determined  by  experimental  tests  which  measure 
the  strength  loss  following  the  thermal  fracture  or  related 
change  in  other  physical  characteristics. 

Again,  the  above  discussion  indicates  the  need  of  a  thorough 
analysis  of  the  thermal  environment,  magnitude  of  stress,  failure 
mode  and  extent  of  crack  propagation  before  a  judicious  choice  of 
the  most  suitable  ceramic  for  any  specific  application  can  be 
made.  For  this  purpose,  especially,  the  qualitative  and  quantita¬ 
tive  aspects  of  crack  propagation  under  conditions  of  thermal 
stress  need  to  be  developed  further  in  considerably  greater 
detai 1 . 

EXPERIMENTAL 

Initiation  of  Thermal  Stress  Fracture 

The  experimental  measurement  of  the  resistance  of  structural 
ceramics  to  thermal  stress  fracture  depends  on  the  information 
desired.  In  view  of  the  qualitative  and  quantitative  uncertain¬ 
ties  in  the  nature  of  heat  transfer  and  associated  thermal 
stresses,  in  practice  the  resistance  of  thermal  stress  fracture 
of  candidate  materials  is  measured  under  conditions  of  heat 


transfer  which  duplicate  the  actual  service  conditions  as  closely 
as  possible.  In  this  manner,  the  critical  temperature  difference 
(AT  )  required  for  failure  in  a  single  quench,  the  number  of 
cycles  required  for  failure  in  thermal  fatigue, or  other  criteria 
of  thermal  stress  resistance  appropriate  for  a  specific  component 
subjected  to  specific  service  conditions  can  be  measured.  In 
industrial  practice,  a  large  number  of  such  tests  were  developed. 
For  practical  purposes,  such  tests  are  entirely  adequate  as  they 
generate  information  essential  for  the  purpose  of  engineering 
design,  performance  prediction  and  failure  analysis.  However, 
since  no  attempts  are  made  to  analyze  the  data  obtained  by  such 
tests  in  terms  of  thermo-elastic  and  failure  theories,  the 
results  obtained  must  be  regarded  as  empirical  only  and  to  be 
appropriate  only  in  the  specific  conditions  of  such  tests.  As  a 
result,  great  care  needs  to  be  taken  in  any  attempt  to  extrapo¬ 
late  the  result  in  other  thermal  conditions  involving  different 
mechanisms  of  heat  transfer,  geometry,  failure  mechanisms  and 
other  relevant  criteria. 

Studies  carried  out  for  the  purpose  of  comparing  predicted 
thermal  fracture  behavior  with  experimental  data  have  met  with 
varying  degrees  of  success.  Coble  and  Kingery  [37]  obtained  rea¬ 
sonable  agreement  between  calculated  and  experimental  data  for 
the  critical  temperature  difference  across  the  wall  thickness  of 
hollow  concentric  circular  cylinders  made  of  porous  aluminum 
oxide  subjected  to  steady-state  heat  flow. 
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Glenny  and  Royston  [38]  as  well  as  Hasselman  [39]  found  fair 
agreement  between  predicted  and  measured  values  for  the  critical 
temperature  difference  required  to  initiate  fracture  in  circular 
rods  of  polycrystalline  aluminum  oxide  subjected  to  a  single 
quench  from  higher  temperature  into  a  water  bath.  Similar  agree¬ 
ment  between  theory  and  experiment  was  found  for  spherical  shapes 
of  polycrystalline  aluminum  oxide  quenched  into  molten  salt  as  a 
quenching  medium  [40].  It  should  be  pointed  out,  however,  that 
the  accuracy  of  predictions  of  results  for  quenching  experiments 
is  handicapped  by  the  considerable  uncertainty  in  the  value  of 
the  heat  transfer  coefficient  for  heated  surfaces  in  fluids  under 
transient  conditions.  In  particular,  the  heat  transfer  coeffi¬ 
cient  for  a  water  bath  is  strongly  affected  by  the  highly  temper¬ 
ature  dependent  effects  of  nucleate  boiling  and  film  formation  as 
indicated  by  the  studies  of  Farber  and  Scorah  [41],  Manson  and 
Smith  [15]  and  Singh,  Tree  and  Hasselman  [42]. 

For  conditions  of  radiative  heat  transfer,  Hasselman  [17] 
found  excellent  agreement  between  calculated  and  measured  values 
of  the  black-body  temperature  required  to  initiate  thermal  stress 
fracture  in  spheres  of  polycrystalline  aluminum  oxide,  initially 
at  room  temperature  and  suddenly  inserted  into  a  furnace  chamber 
at  thermal  equilibrium  at  temperatures  at  which  radiation  repre¬ 
sents  the  primary  mechanism  of  heat  transfer. 

For  ceramic  plates  subjected  to  forced  convective  heat 
transfer  over  a  central  spot,  Faber,  Huang  and  Evans  [43] 
obtained  excellent  agreement  between  calculated  and  predicted 
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values  of  the  critical  temperature  difference  required  to  initi¬ 
ate  failure  at  a  given  rate  of  gas  flow.  The  good  agreement 
obtained  was  facilitated  by  the  calibration  of  the  heat  transfer 
coefficient  by  use  of  a  ceramic  with  accurately  known  failure 
stresses.  These  latter  values  were  established  by  indenting  the 
plates  by  a  diamond  hardness  indentation  followed  by  an  indepen¬ 
dent  strength  test.  By  this  method  the  statistical  aspect  of  the 
brittle  fracture  is  removed.  This  overall  approach  must  be  con¬ 
sidered  unique  in  that  the  magnitude  of  thermal  stress  is  estab¬ 
lished  from  the  fracture  process  itself.  The  magnitude  of  heat 
transfer  coefficient  obtained  from  these  data  can  then  be  used  to 
calculate  the  stresses  in  other  ceramic  plates  to  be  tested. 

The  prediction  of  thermal  fatigue  behavior  of  brittle  ceramics 
has  met  with  a  varying  degree  of  success.  Figure  4  compares  the  ob¬ 
served  and  predicted  number  of  cycles  to  failure  of  circular  rods  of 
a  soda-lime-silica  glass  quenched  into  water  measured  by  Hasselman, 
et  al.,  [44].  However,  it  should  be  noted  that  the  predicted 
thermal  fatigue-life  shown  in  Fig.  4  required  the  use  of  the  Wei- 
bull  theory  in  order  to  take  into  account  that  the  failure  prob¬ 
ability  under  conditions  of  the  multi-axial  thermal  stress  field 
differs  from  the  corresponding  value  under  conditions  of  the  uni¬ 
axial  strength  test.  These  differences  in  failure  probabilities 
lead  to  the  effect  that  the  mean  fracture  stress  under  conditions 
of  the  thermal  quench  is  far  lower  than  the  tensile  fracture 
stress  under  conditions  of  uniaxial  loading.  It  should  be  recog¬ 
nized  that  the  formulation  of  the  Weibull  theory  must  be  consid- 


ered  empirical.  For  this  reason,  the  reasonable  agreement  shown 
in  Fig.  4  to  some  extent  may  well  be  fortuitous.  It  should  also 
be  noted  that  even  slight  uncertainties  in  the  kinetics  of  sub- 
critical  crack  growth  can  introduce  large  relative  errors  in 
predicted  thermal  fatigue-life.  Nevertheless,  the  agreement 


indicated  by  the  observed  and  predicted  results  shown  in  Fig.  4 
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offers  encouragement  towards  future  work  in  this  area. 

The  thermal  fatigue  behavior  of  polycrystalline  hot-pressed 
silicon  nitride  [45]  showed  good  agreement  between  predicted  and 
observed  values  for  ranges  of  crack  sizes  expected  to  occur  in 
this  material.  Nevertheless,  thermal  fatigue-life  was  found  to 
depend  critically  on  the  assumption  of  the  value  of  crack  size  in 
the  original  material  prior  to  thermal  fatigue.  Non-destructive 


tests  for  the  accurate  measurement  of  flaw  size  of  the  order  of 
50-100  pm  have  not  been  developed  at  this  stage  of  technical 
development.  For  this  reason,  the  prediction  of  the  thermal 
fatigue  behavior  of  structural  ceramics  such  as  silicon  nitride 
must  be  regarded  as  still  relatively  uncertain. 


Observations  on  the  Nature  of  Crack  Propagation  During 
Thermal  Stress  Fracture 

In  view  of  its  importance  to  engineering  practice,  crack 
propagation  following  the  onset  of  thermal  stress  fracture  has 
received  quite  a  bit  of  attention.  Literally  explosive  fracture 
of  aluminum  oxide  and  zirconium  carbide  subjected  to  high  radia¬ 
tive  heat  fluxes  was  observed  by  Crandall  and  Ging  [46]  and  Has- 
selman,  Becher  and  Mazdiyasni  [  l  ],  respectively.  In  these  stu- 
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Predicted  and  experimental  cycles-to-failure 
of  soda-lime-silica  glass  rods  with  initial 
flaw  depth  of  8.6  ym  subjected  to  thermal 
fatigue  by  quenching  into  water  bath  at  33  C 
(after  Hasselman  et  al.  [44]). 
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dies,  the  samples  were  fractured  in  a  number  of  separate 
fragments.  Clearly,  this  fracture  mode  is  highly  undesirable  as 
the  crack  sizes  exceeded  the  dimensions  of  the  materials  being 
tested.  This  mode  of  failure  is  quite  commonly  encountered  in 
practice  especially  for  components  of  relatively  large  size. 

Under  other  conditions  of  thermal  stress  fracture  such  as  a 
quench  into  water  or  other  appropriate  fluids,  the  cracks  formed 
will  arrest  within  the  material.  Such  cracks,  because  of  little 
or  no  permanent  crack-opening-displacement,  frequently  are  not 
easily  detected  unless  by  highly  specialized  laboratory  techni¬ 
ques  such  as  electron-optical  method.  Nevertheless,  such  cracks 
can  have  a  significant  effect  on  load-bearing  characteristics. 
For  this  reason  and  for  convenience,  experimental  studies  of  the 
nature  of  crack  propagation  have  concentrated  on  measurements  of 
load-bearing  characteristics,  specifically  tensile  strength  as  a 
function  of  severity  of  the  conditions  of  the  thermal  shock  test. 

The  crack  stability  and  propagation  behavior  shown  in  Figure 
2,  derived  for  a  simple  mechanical  model  with  non-interacting 
cracks,  should  be  qualitatively  valid  also  for  other  stress  con¬ 
ditions  involving  tensile  failure.  As  can  be  inferred  from  Fig. 
2,  Fig.  5  depicts  schematically  the  expected  relative  dependence 
of  crack  size  and  corresponding  tensile  strength  for  unstable  and 
stable  crack  propagation  as  a  function  of  severity  of  thermal 
shock,  expressed  in  such  terms  as  AT  in  a  quenching  experiment, 
rate  of  heat  flow  or  other  appropriate  measure.  For  unstable 
crack  propagation,  crack  size  and  strength  exhibit  a  discontinu- 
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ous  increase  and  decrease,  respectively.  The  crack  size  which 
results  at  AT  is  sub-critical  and  will  not  change  in  size  until 
AT  reaches  the  value  AT'.  For  AT  >  AT'  ,  the  crack  will  propagate 
in  a  stable  manner. 

In  contrast,  for  stable  crack  propagation,  for  A  T  >  AT  , 

c 

crack  and  strength  will  increase  and  decrease,  respectively  in  a 
monotonic  manner  as  shown  in  Figs.  5c  and  5d. 

The  general  strength  behavior  for  unstable  crack  propagation 
shown  in  Fig.  5b  was  verified  experimentally  by  a  number  of 
investigators  for  a  variety  of  structural  ceramics  such  as  alumi¬ 
num  oxide,  silicon  nitride,  berylium  oxide,  silicon  carbide, 
boron  carbide,  zinc  oxide  and  soda-lime  silica  glass  [48-55]. 

Figure  6  shows  typical  data  for  an  industrial  grade  poly¬ 
crystalline  aluminum  oxide  subjected  to  a  quench  into  a  water 
bath  at  room  temperature  followed  by  a  strength  test.  The 
instantaneous  loss  in  strength  at  aTc  as  well  as  the  region  of 

sub-criticality  of  the  cracks  for  AT  <  AT  <  aT'  as  shown  in  Fig. 

c  c 

5b  is  clearly  evident. 

In  the  analysis  of  unstable  crack  propagation  it  was  pred¬ 
icted  that  the  strength  retained  following  the  initiation  of 
thermal  stress  fracture  was  an  inverse  function  of  the  strength 
prior  to  fracture.  In  this  respect,  it  is  of  interest  to  compare 
the  experimental  data  shown  in  Fig.  7  for  another  industrial  alu¬ 
mina  obtained  under  identical  quenching  conditions  as  the  data 
shown  in  Fig.  6. 


TEMPERATURE  DIFFERENCE 


Schematic  representation  of  the  crack  length 
and  corresponding  strength  behavior  as  a 
function  of  severity  of  thermal  shock  for 
unstable  (A,B)  and  stable  (C,D)  crack  propagation 
(after  ref.  [47] ) . 
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Tensile  strength  of  an  industrial  poly¬ 
crystalline  aluminum  oxide  following  a 
quench  from  higher  temperature  into  a 
water  bath  as  a  function  of  quenching 


temperature  difference  (after  Hasselman, 
1970)  (For  further  details  see  ref.  [39]) 


TEMPERATURE  DIFFERENCE 


Tensile  strength  of  a  lower  grade  polycrystalline 
aluminum  oxide  subjected  to  water  quench,  showing 
a  much  lower  relative  strength  loss  at  ATC  than 
for  the  alumina  shown  in  Fig.  6  (For  further 
details  see  ref.  [39]). 


A  comparison  of  the  data  shown  in  Figs.  6  and  7  shows  that 
the  alumina  of  Fig.  6  has  a  higher  initial  strength  in  the 
untested  condition  (i.e.  AT  =  0)  than  the  alumina  of  Fig.  7.  In 
view  of  the  role  of  strength,  it  is  expected  that  a  higher 
quenching  temperature  difference  is  required  to  initiate  failure 
in  the  initially  stronger  alumina  than  in  the  weaker  alumina. 
This  conclusion  is  confirmed  by  the  experimental  data  for  the 
quenching  temperature  difference  required  to  initiate  failure 
which  is  about  300°C  for  the  stronger  alumina  and  around  225°C 
for  the  weaker  alumina. 

From  the  role  of  strength  in  the  extent  of  crack  propaga¬ 
tion,  as  analyzed  earlier,  it  is  predicted  that  the  strength 
retained  following  the  initiation  of  thermal  failure  is  an 
inverse  function  of  the  original  strength.  The  validity  of  this 
conclusion  is  easily  confirmed  by  the  strength  for  aT  >  ATc  given 
in  Figs.  6  and  7.  In  general,  these  data  are  illustrative  of  the 
diametrically  opposite  roles  of  strength  in  the  resistance  of  the 
initiation  of  thermal  stress  fracture  and  the  extent  of  crack 
propagation  and  resulting  degree  of  damage  following  the  onset  of 
crack  propagation.  Considering  the  two  aluminas  of  Figs.  6  and 
7,  clearly  for  thermal  conditions  involving  thermal  stresses  of 
low  magnitude  and  failure  probability,  the  aluminum  oxide  of  Fig. 

6  is  the  preferred  candidate  in  view  of  its  superior  load-bearing 
characteristics  possibly  required  for  its  intended  purpose.  On 
the  other  hand,  for  magnitudes  of  thermal  stresses  for  which  fai¬ 
lure  is  initiated  in  both  aluminas,  the  alumina  of  Fig.  7  is  pre- 


ferred  in  view  of  its  superior  strength  behavior  following  the 
onset  of  thermal  stress  failure. 

Figure  8  shows  the  strength  behavior  typical  of  stable  crack 
propagation  as  found  by  Larson  et  al.  [47]  for  a  high-alumina 
refractory  material  subjected  to  a  water-quench.  Note  that  ori¬ 
ginal  strength  (  AT  =  0)  of  the  material  is  very  low  indeed  com¬ 
pared  to  strengths  commonly  found  in  other  general  structural 
materials.  This  low  strength  value,  in  terms  of  the  simple  ana¬ 
lysis  presented  earlier,  assures  that  the  size  of  failure-initi¬ 
ating  cracks  is  sufficiently  large  that  the  unstable  mode  of 
crack  propagation  is  suppressed  and  that  crack  propagation  will 
occur  only  in  the  stable  mode.  This  latter  mode  of  crack  propa¬ 
gation  is  preferred  because  at  a  given  value  of  AT,  as  easily 
ascertained  from  Fig.  1,  the  extent  of  stable  crack  propagation 
is  less  than  for  unstable  crack  propagation. 

Critical  to  note,  again  as  easily  ascertained  from  Fig.  2, 
is  that  whether  a  crack  will  propagate  in  a  stable  or  unstable 
manner  depends  on  the  total  number  of  cracks  which  participate 
simultaneously  in  the  fracture  process.  This  implies  that  the 
mode  of  crack  propagation  and  associated  strength  loss  for  any 
given  material  will  depend  on  the  nature  of  the  thermal  shock 
environment.  Indeed,  Larson  et  al  [56]  found  that  samples  of  a 
high-alumina  refractory  when  subjected  to  a  water  quench  failed 
by  stable  crack  propagation,  whereas  samples  of  the  same  material 
when  subjected  to  thermal  shock  by  radiation  heating  failed  in 
unstable  mode.  These  differences  in  fracture  mode  can  be  traced 
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8.  Strength  loss  typical  for  stable  crack 
propagation  in  high-alumina  refractory 
subjected  to  quench  into  water  bath 
(after  ref.  [47]) . 
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directly  to  differences  in  crack  density.  During  the  water 
quench  tensile  failure  is  nucleated  at  numerous  surface  sites, 
because  the  total  surface  is  subjected  to  a  high  value  of  tensile 
stress.  In  contrast,  during  radiation  heating,  the  maximum  ten¬ 
sile  stress  is  found  along  a  single  line  or  plane  so  that  fracture 
will  originate  at  a  limited  number  of  sites  only.  Again,  prior  know¬ 
ledge  of  the  mode  of  heat  transfer  must  be  established  before  a 
candidate  material  can  be  selected  on  the  basis  of  the  nature  of 
crack  propagation. 

Figures  6  and  7  clearly  indicate  evidence  for  unstable  crack 
propagation.  It  is  critical  to  point  out  that  such  strength  loss 
behavior,  as  found  by  Ziegler  and  Hasselman  [ 60  ],  can  be  masked 
by  excessive  data  scatter  to  exhibit  an  apparent  stable  failure 
mode  as  shown  in  Figs.  5  and  8. 

The  dependence  of  AT  on  crack  length  as  shown  in  Fig.  2 

also  predicts  that  at  the  transition  from  unstable  to  stable 

crack  propagation  ( at  £  =  £  ),  AT  should  be  independent  of  crack 

o  m  c 

length.  This  conclusion  was  confirmed  by  Gupta  [48]  for  a  number 
of  polycrystalline  aluminum  oxides  with  a  range  of  strength.  For 
all  these  sets  of  data,  AT^  was  found  to  be  independent  of 
strength,  with  the  high  strength  samples  exhibiting  unstable 
crack  propagation  and  the  low  strength  samples  failing  in  a  sta¬ 
ble  mode. 


Finally,  Fig.  2  suggests  that  materials  with  very  high  den¬ 
sity  of  micro-cracks  should  exhibit  exceptional  thermal  shock 
resistance.  Indeed,  this  prediction  is  verified  by  the  dramatic 


results  of  Rossi  [57]  for  a  polycrystalline  magnesium  oxide  in 
which  micro-cracking  was  induced  by  a  dilute  dispersed  phase  of 
spherical  particles  of  tungsten.  In  the  resulting  composite,  the 
micro-cracking  resulted  from  the  large  mismatch  in  the  coeffi¬ 
cients  of  thermal  expansion  of  the  magnesium  oxide  matrix  and 
tungsten  dispersed  phase.  It  may  be  added  that  micro-cracking 
can  cause  an  appreciable  reduction  in  thermal  conductivity  [58]. 
For  this  reason,  micro-cracked  materials  should  prove  to  be  par¬ 
ticularly  suited  to  design  situations  having  the  simultaneous 
requirement  of  high  thermal  shock  resistance  in  combination  with 
high  thermal  insulating  ability  [59]. 

FINAL  COMMENTS 

It  appears  reasonable  to  conclude  that  the  principles  which 
underlie  thermal  stress  fracture  have  been  developed  to  a  consid¬ 
erably  greater  degree  than  the  predictions  of  such  fracture, 
which  in  practice  are  handicapped  by  uncertainties  in  the  quali¬ 
tative  and  quantitative  aspects  of  heat  transfer  as  well  as  un¬ 
certainties  in  the  relevant  material  properties.  The  high  degree 
of  brittleness  of  ceramic  materials  exacerbates  the  problem  even 
further  as  the  safety  features  of  non-linear  deformation  found 
for  many  metals  and  polymers  is  entirely  absent.  Such  properties 
as  the  coefficient  of  thermal  expansion.  Young's  modulus,  Pois¬ 
son's  ratio,  thermal  conductivity  and  optical  properties  which 
are  governed  by  inter-atomic  effects  are  under  control  of  the 
material  specialist  to  a  very  limited  degree  only.  For  this  rea¬ 
son,  improvements  in  the  resistance  of  brittle  ceramics  to  frac- 
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ture  must  rely  primarily  on  improvements  of  strength  and  fracture 
toughness  (energy).  Surface-compression-strengthening  and  the 
formation  of  materials  with  thermal  conductivity  gradients  may 
also  be  used  to  advantage  (61-65].  Alternatively,  for  structural 
ceramics  available  at  present,  design  modifications  may  be 
required  in  order  to  reduce  the  magnitude  of  thermal  stress  to 
acceptable  levels.  Further  understanding  of  the  nature  of  crack 
propagation  in  steady-state  or  transient  thermal  stress  field  can 
be  enhanced  by  additional  theoretical  work.  This  latter  area  is 
judged  to  be  of  particular  practical  importance  as  it  should 
enhance  the  reliability  of  design  of  components  or  structures 
made  from  brittle  ceramics  subject  to  thermal  stress  failure. 
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